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ABSTRACT 

Atmospheric  temperature  structure  fluctuations  quantified  by  CT2 

degrade  the  spatial  and  temporal  coherence  of  electromagnetic  and 
acoustic  waves  propagating  in  the  atmosphere.    A  computer  controlled 
atmospheric  echosounder  developed  at  the  Naval  Postgraduate  School 
measures  a  time  averaged  CT2  profile  of  the  lower  atmosphere.  Assigning 

the  proper  CT2  values  to  the  backscattered  return  signals  depends  on  an 

accurate  calibration  of  this  instrument.    Calibration  involves 
determination  of  the  product  of  the  echosounder"s  transmission  efficiency 
Et  and  reception  efficiency  E 

This  thesis  provides  a  preliminary  investigation  of  a  calibration 
process  using  pulsed  acoustic  energy  backscattered  from  hard  spheres. 
Supporting  software  calculates  the  desired  product  ErEt  based  on  an 

assumption  of  echosounder  efficiency  reciprocity.  Results  of  the 
calibration  process  investigation  indicate  this  assumption  may  be  invalid. 
The  results  also  indicate  the  software  performs  as  intended  and  that  the 
proposed  calibration  method  possesses  sufficient  merit  to  warrant  further 
development. 
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I.    INTRODUCTION 

The  variational  effect  of  wind  shears,  convection,  and  temperature 
gradients  on  atmospheric  density  distributions  characterize  atmospheric 
turbulence.  The  temperature  structure  parameter,  CT2,  and  the  index  of 

refraction  structure  parameter,  C  2,  quantify  two  density  distributions 

that  arise  in  acoustic  and  electromagnetic  propagation.    [Ref.  1] 

The  density  distributions  quantified  by  Cn2  and  CT2  induce  phase 

fluctuations  in  propagating  electromagnetic  waves.  The  density 
distributions  of  atmospheric  turbulence  degrade  the  spatial  and  temporal 
coherences  of  initially  coherent  electromagnetic  radiation.    Similar,  more 
pronounced  effects  occur  in  the  propagation  of  acoustic  waves.    [Ref.  1] 

A  number  of  techniques  for  correcting  these  problems  have  been 
proposed,  particularly  with  respect  to  electromagnetic  propagation.  These 
corrective  techniques  include  the  use  of  modular  mirrors  [Ref.  2],  optical 
phase  conjugation  [Ref's.  2  and  3],  optical  and  digital  signal  processing 
[Ref.  4],  and  simple  avoidance  of  turbulent  air  masses.  The  appropriate 
timing  of  transmissions  or  suitable  location  of  transmission  and  reception 


facilities  accomplishes  avoidance  [Ref's.  5  and  6].    These  techniques  can 
be  applied  alone  or  in  combinations. 

A  knowledge  of  the  history  of  the  density  or  refractive  structure 
along  a  proposed  transmission  path  is  required  to  effectively  employ  a 
number  of  these  compensatory  techniques.    Time-averaged  CT2  or  Cn2 

profiles  charactize  this  history  [Refs.  1  and  71 

Wroblewski  and  Weingartner  developed  and  Moxcey  refined  a  high 
resolution,  computer  controlled  acoustic  sounder,  or  echosounder,  system. 
Their  computer  controlled  echosounder  system  uses  a  transmitted  pulse  of 
acoustic  energy  as  a  probe  of  atmospheric  structure.  The  interaction  of 
the  acoustic  wave  packet  with  atmospheric  turbulence  is  used  to  develop  a 
time-averaged  Cj2  profile  for  a  short  spatial  range  [Ref's.  5,  6  and  81  The 

acoustic  sounder  may  also  be  used  to  verify  certain  long  range  optical 
measurements  of  turbulence  [Ref.  61 

Development  of  the  CT2  profile  uses  a  CT2  expression  developed  from 

the  echosonde  equation  summarized  by  Neff  [Ref.  7]  and  the  empirical 
acoustic  backscatter  cross  section  per  unit  volume  expression  of  Tatarski 
[Ref.  9)  [Refs.  5  and  61  The  accuracy  of  the  absolute  values  of  this  CT2 

profile  depends  on  the  calibration  of  the  acoustic  antenna  parameters 


contained  in  the  CT2  expression     The  antenna  parameters  of  interest  are 

£.3  the  efficiency  of  conversion  of  electrical  to  acoustic  power,  and  E  ,  the 

efficiency  of  conversion  of  acoustic  power  to  electrical  power. 

In  this  thesis  acoustic  energy  backscattered  off  acoustically  hard 
spheres  was  used  to  determine  Er  and  E..    Acoustic  range  equations  were 

software  implemented  to  perform  the  necessary  calculations.  Data  and 
calculational  results  provided  first  order  verification  of  the  calibration 
process. 


II    THEORETICAL    BACKGROUND 


A.      CT2  REVIEW 


Wroblewski  and  Weingartner  applied  the  echosonde  equation, 
Pr-Er[PtEt][e-2«R][^0o(R,f)][AG]_      (.: 


summarized  by  Neff  [Ref.  7]  and  the  acoustic  backscattering  cross 
section  per  unit  volume  expression  given  by  Tatarski  [Ref.  9]  as 

V    c  2 

<Jn(R,f)=   0.0039  k'3   -3-  C2) 

T 
1  o 

to  develop  the  volume  averaged  expression 

CT2=— ! !_   lol  A_LfLR2e2ocR      (3) 

0.0039  ErEt     k/3    cr    AG  Pt 

where 

•  Pr  is  the  electrical  power  of  the  reflected  signal  returned  to  the 

acoustical  antenna  array, 

•  Ep  is  the  efficiency  of  conversion  of  the  returned  acoustic  power  to 

•  Pt  is  the  electrical  power  supplied  to  the  acoustical  array, 

•  Et  is  the  efficiency  of  conversion  of  Pt  to  acoustic  power, 

•  ex  is  the  average  attenuation  per  unit  distance  along  the  transmission 
path, 

•  R  is  the  range  from  the  array  to  the  target  or  backscattering  volume, 


•  c  is  the  average  speed  of  sound  along  the  transmission  path, 

•  t  is  the  length  in  time  of  the  transmitted  acoustic  pulse, 

•  A  is  the  antenna"s  aperture  area, 

•  6  is  the  antenna's  effective  aperture  factor, 

•  k  is  the  wavenumber  of  the  incident  acoustic  energy,  and 

•  TQ  is  the  average  temperature  (in  degrees  Kelvin)  along  the 

transmission  path. 
[Ref's.  5  and  6] 

The  conditions  at  the  time  and  location  of  each  application  determine 
the  value  of  the  terms  T  ,  k,  c,  t,  Pr,  Pt,  R  and  <x  in  the  CT2  expression. 

The  terms  A  and  G  depend  on  the  antenna  design. 

Adverse  field  effects  such  as  dust,  debris  and  vibration  loosened 
connections  cause  the  remaining  two  terms,  Er  and  EL,  to  change  over  the 

life  of  the  system.    This  chapter  develops  the  theoretical  expressions 
supporting  a  determination  of  Ep  and  Et. 

B.      ACOUSTIC  RANGE  EQUATIONS  AND  EFFICIENCIES 

In  the  following  two  sections  the  acoustic  range  equation  of  Neff 
[Ref.  7]  is  re-developed  for  targets  of  small  area  relative  to  the 
cross-section  of  the  ensonifying  field.  Neff  [Ref.  7],  Skolnik  [Ref.  10],  and 
Probert-Jones  [Ref.  1 1]  were  used  for  guidance. 


1.      One  Way  Range  Equation  and  EtGo 

As  defined  previously,  the  electrical  power  supplied  to  the 
acoustic  array  is  Pt.    It  is  converted  at  an  efficiency  Et  to  an  emitted 


acoustic  power, 


P  =  P  E 
ra     rt  ct 


(4) 


P  is  spread  over  a  4tt  solid  angle. 


A  far  field  range  R  is  defined  by  the  condition 


A/2  »  |L|  -  IRI 


(5) 


for  the  geometry  in  Figure  1.    A  is  the  wavelength  of  the  emitted  acoustic 


wave. 


H-Rl 


Acoustic 
Array 


(0,0) 


Tarqet 

of 

Radi  us  a 


Figure  1:  Calibration  (Far  Field)  Geometry 


The  acoustic  intensity, 

la  =  Pa  /  4TTR2,  (6) 

describes  the  emitted  acoustic  power  density  of  an  isotropic  emitter  at  a 
far  field  range  R.  IQ  is  in  units  of  [W/m2]. 

If  the  emitted  power  is  directional,  the  intensity  la  is  modulated 

by  a  geometrical  gain, 

G(O)  =  GolF(O)|2  =  GolF(0,e)|2  (7) 

Note  this  is  not  the  same  G  found  in  equation  ( 1 ).    Now, 

|F(0,0)|2=1,  (8) 

and  G   is  thef ore  the  centerline  gain  of  the  acoustic  array.  [Ref.  11]    The 

centerline  direction  (0,0)  is  taken  to  be  normal  to  the  face  of  the  acoustic 
array.    The  acoustic  intensity  lQ  is  now  written 


R 


a 


'a 


4t)R2 


G(O)  (9) 


In  a  real  propagating  medium  intensity  attenuation  occurs  along 
the  transmission  path.     The  radial  term  e"ocR  accounts  for  this 

attenuation  when  the  atmospheric  path  is  assumed  to  be  reasonably 


7 


homogeneous.  The  acoustic  intensity  is  then 


la   =   -^J5-_r    G(0)e"aR  (io) 

4TTR" 


A  target  of  cross  sectional  area  At  t  intercepts  a  power 


Expanding  this  expression  with  equation  (10)  yields 

For  a  target  of  cross  sectional  area  At .  with  a  sufficiently 

small  solid  angle  At  t  /R2  little  change  in  intensity  occurs  about  the 

direction  Q.    Equation  ( 1 2)  is  now  arranged  as  a  one  way  range 
equation  along  the  centerline  direction  (0,0), 

R  "  L   471  b°    P.  .  J  . 


(13) 


Isolating  EtGo  in  the  above  expression  yields 


4rtR2      <xr    %i_  (M) 

CtU0    "  p  C  A 


8 


Electrical  power  is  expressed  as 


2  w 


p    -       Vrms         _       ^m;  (15) 


where 


•  V       is  the  empirical  root  mean  square  voltage, 

•  V     is  V     2,  and 

ms  rms  ' 

•  Z  is  the  electrical  impedance. 
EjGq  is  now  written 

E  G    =    4TJR2       z        eaR    -&L  <]6> 

«      °  (Vm.<:)te  Atgt      . 

The  ratio  Pt  t/At  t  is  defined  as  the  acoustic  intensity  lr 
received  by  the  target.    Therefore, 

E  g    =    4T1R2      z        eaR    1  <17> 

2.      Two  Way  Range  Equation  and  Er 

For  an  acoustically  hard  target  the  total  power  reflected  is  P.  t 

The  power  reflected  from  the  target  surface  has  a  normalized, 
directional  intensity  distribution  described  by  the  differential  scattering 
cross  section,  <r 


The  direction  from  the  target  to  the  array  is  the  backscattered 
direction.    The  differential  scattering  cross  section  in  the  backscattered 
direction  is  called  the  normalized  backscattered  cross  section,  <rb.  <rb  is 

defined  as  the  power  reflected  toward  the  source  per  unit  solid  angle, 
normalized  by  the  incident  power  density  over  a  4tt  solid  angle. 
The  backscattered  intensity  from  the  target  is 

An  array  of  aperture  area  A  at  range  R  intercepts  an  attenuated,  received 
acoustic  power, 

Pra  =  'tgt  A  e_KR  •  (,9) 

Substituting  equation  ( 1 8)  into  equation  ( 1 9)  yields 

P     =  ^Sl_    -     e"«R  A  (20) 

r°        4«R2      b 

PrQ  is  converted  to  returned  electrical  power,  P„,  at  a  return 

ra  r 

efficiency,  Er; 


P    =  P     E  (21) 

r       ra    r 


0 


Since  electrical  power  is  expressed  as 


a- 


Vrms         _        *  rns 


electrical  power  representing  the  returned  acoustic  signal  from  a 
target  placed  along  an  emitter's  direction  (0,0)  is  written  as 

(Vms)r .  -   (Vms)t  [  Et  GQ    g-OCR  Atqt"|rErffb     e"«R    A    H      (22) 
Zr  Zt      L    47J  R2    JL    4Tf  R2  J    ■ 

For  a  passive,  linear  acoustic  array  circuit  the  concept  of 
impedance  reciprocity  is  invoked  [Ref.  12]  to  yield  a  two  way  range 
equation  for  scattering  from  a  small  target; 


R      r(Vms)i  Ml   e-2«R   A  A      ]4       (23) 

I    (V      }        t  Lr    ,        /2  toil 

This  is  written  to  isolate  E  ■ 


E     _     (Vms)r    (4tJR2)2      e2<xR    j_ 
r  "     (Vrns)t       Et60        A  Afgt    ab     . 


(24) 


A  part  of  the  return  signal's  power  is  unwanted  noise. 
Subtracting  noise  power  from  the  returned  acoustic  power  yields 


r  (Vms)r-(Vms)n    (4T)R2)2       e.2°:R     J_  (25) 


11 


As  return  voltages  are  quite  small  they  are  analyzed  after 
electrical  amplification.    Calling  the  electrical  amplification  a  gain 
Ga  allows  introduction  of  the  final  modification  to  E  : 

e  r 

E      _    (Vms)r-(Vms)n    (4rtR2)2       e2<XR     _L  (26) 

r  "       (Vms)f  Ge2  EtG0        A  Aigt     %     - 

C.      ATTENUATION  AND  RANGE  DETERMINATION 

Employing  equations  (17)  and  (26)  in  the  calibration  process  requires 
calculating  the  attenuation  a  and  the  range  R  from  the  array  to  the 
target. 

1.      Attenuation 

The  calculated  attenuation  a  for  an  atmospheric  propagation 
path  is  assumed  to  be  the  combined  result  of  molecular  and  classical 
absorption  in  the  atmosphere.    The  classical  absorption  a.  is  attributed 

to  heat  conduction  and  viscous  effects  while  the  larger  molecular 
absorption  amoI  is  attributed  to  the  excitation  of  internal  energy  modes  of 

atmospheric  gases  by  the  propagating  sound  energy.    [Ref.  7] 

Information  from  Neff  [Ref.  7],  Businger  [Ref.  13],  Neiburger 
[Ref.  14],  and  Fuller  [Ref.  15]  is  applied  to  make  a  reasonable  determination 
of  a.    The  determination  is  made  using  the  average  measured  temperature 

12 


Tc  in  degrees  Celsius  along  the  propagation  path,  the  measured 

atmospheric  pressure  P  in  millibars,  and  the  measured  relative  humidity 
Rh  in  percent. 

The  molecular  attenuation  coefficient,  a    ,,  is  expressed  by  the 
empirical  relationship 

9  ± 

OC. 


I  r-a+in 


TTI8X 

°'mo1        3048  'ratio 


ocmol  is  in  units  of  [-dB/m]. 


Now, 


amax  =  0.0078  fm  (T*r25  e777(M/T*>  (28) 


is  the  maximum  absorption  at 

(10  +  6600  (lOOf-)  +  44,400  (]OQff)P*       (2g) 
fm  =  (T*)08  . 

f    is  in  units  of  [Hz]. 

m 


Also. 


e_  =       E5Rh  (30) 

P        P  -  Es(  1  +  Rh) 

is  the  mole  ratio  of  water  vapor  in  the  atmosphere.  The  water  vapor 
pressure  e  is  expressed  in  millibars. 


E         10(9.4-2353/To)  (3 ,  j 

expresses  the  saturation  vapor  pressure  E      TQ  is  the  average  absolute 

temperature  along  the  transmission  path,  as  defined  in  Section  A. 

T*  =  (1.8TC+  492)/  519,  (32) 

and 


P*  =  P/  1014.  (33) 


Finally, 


fratio  "  f  I  fm  ■  <34> 

where  f  is  the  operating  frequency  of  the  echosounder  in  [Hz].    [Ref's.  7, 
13,  14and15] 

The  classical  attenuation  coefficient,  occ),  is  approximated  by 

<xcl=  1.74*  1(T10(f)2  .  (35) 

cccl  is  in  units  of  I-dB/ml  [Ref.  7] 

For  a  total  attenuation 

«t  =  acl  +  amol  (36) 

and  the  relationship 
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-aTR  =  10  log  e_ccR  (37) 

the  average  attenuation  per  meter  is  finally  written  as 

«   =    "*  L  "™]    mio  (38) 

Equation  (37)  is  derived  from  unit  definitions  and  equation  one  of  Neff 
[Ref.  7]. 

2.     Range 

The  coordinate  system  for  the  backscattered  cross  section 
determination  uses  the  center  of  the  target  spheres  as  the  coordinate 
system  origin.  Therefore,  the  radius  a  of  a  target  sphere  and  the  time  of 
return  tr  of  an  acoustic  signal  from  the  acoustic  echosounder  are  used  to 

determine  the  range  R  from  the  target  to  the  echosounder; 

R  =  -^+8  .  (39) 

c  is  the  average  speed  of  sound  along  the  transmission  path,  and 

r  "    receive  ~    transmit  ■  *      ' 

For  the  range  of  temperatures  expected  in  normal  echosounder  operation 

the  speed  of  sound  in  meter  per  seconds  in  dry  air  is 

cdry  =  2005^V  <4I> 


For  moist  air  c  becomes 

c   -   cmo,st  "   cdry"  ♦0.14ff/P),  (42) 

where  p/P  is  defined  in  equation  (30).    [Ref.  1] 

D.      BACKSCATTERED  CROSS  SECTION  DETERMINATION 

Employing  equations  (17)  and  (26)  in  the  calibration  process  also 
requires  calculating^,  the  backscattered  cross  section.  Recall  from 

Paragraph  B.2  that  the  backscattered  cross  section  is  the  differential 
scattering  cross  section  in  the  direction  of  the  energy  source. 
1.      Initial  Conditions  and  Assumptions 

In  determining  the  differential  scattering  cross  section  o-  it  is 
assumed  the  incident  wave  state  is  unchanged  after  scattering.  That  is, 
for  a  stationary,  acoustically  hard  target  it  is  assumed 

k   =   incident  :     ^scattered  •  (43) 

and  that  incident  k  has  good  definition,  i.e.  the  incident  wave  can  be 
considered  monochromatic  [Ref.  16]. 

Selection  of  acoustically  hard  spheres  as  targets  further 
simplifies  the  determination  of  the  differential  scattering  cross  section 
solution.    This  target  selection  sets  the  boundary  condition  at  the  target 
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surface  as  the  Neumann  boundary  condition 


iif=° 


(44) 


V  is  the  velocity  potential  defined  by  the  acoustic  propagation  velocity  c 


as 


C  =  7K 


(45) 


n  lies  along  the  direction  normal  to  the  target  surface.    Also, 


p-*# 


(46) 


for  the  pressure  p  and  rest  density  p  of  the  propagating  medium.  [Ref.  16]. 


Yte 


w/71 


Incident  Wave 
Packet  of  N  Waves, 

W.NA»a. 


^ 


Target  of 
Radi  us  a . 


Target  Center 
is  Origin  of  the 
Scattering 
Coordinate 
System. 


Figure  2:  Incident  Wave  and  Target  Description 


The  target  sphere  is  aligned  along  the  centroid  of  the  main 
acoustic  lobe  of  the  acoustic  arrray  echosounder,  that  is,  along  the 
direction  (0,0),  normal  to  the  face  of  the  array.    The  target  sphere  is 
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selected  to  be  sufficiently  small  with  respect  to  incident  beam 
dimensions  so  that  an  essentially  uniform  intensity  is  intercepted  by 
the  sphere  [Ref.  171    Since  the  target  intercepts  a  very  small  proportion 
of  the  centroid  of  the  transmitted  sounding  lobe  no  accommodation  for 
sounding  lobe  geometry  over  the  target's  angular  area  is  made. 
Acoustic  gain  is  simply  taken  as  GQ. 

Target  placement  is  also  at  a  far  field  range.  This  allows 
spherical  incident  waves  to  be  treated  as  planar  incident  waves  and 
places  the  target  outside  the  range  associated  with  the  ring  time  of  the 
echosounder  [Ref.  8].  Placement  of  the  target  is  also  close  enough  to 
the  array  to  support  an  assumption  of  transmission  medium  homogeneity 
along  the  transmission  path. 

Finally,  the  potential  describing  the  target  falls  off  faster  than 
1  /r.    The  describing  potential  for  an  acoustically  hard  sphere  with 
radius  a  is 

r    °°     r  i  a 
0     D  a 


V(r)  = 


(47) 


Thus,  for  r  >  a  V(r)  falls  off  faster  than  1/r. 


Note  that  with  the  exception  of  the  specific  surface  boundary 
condition  for  acoustic  waves  the  above  conditions  and  assumptions  closely 
correspond  to  those  for  plane  electromagnetic  waves  incident  on  a 
perfectly  conducting  sphere,  and  for  perfectly  elastic  nuclear  scattering 
in  a  center  of  mass  coordinate  system.  [Ref's.  16,  18,  19  and  20] 
2.      The  Incident  Wave 

The  asymptotic  form  of  an  incident  plane  wave  approaching  a 
scattering  target  but  beyond  the  influence  of  the  target's  potential  V(r)  is 
the  undisturbed,  or  free,  plane  wave.    In  the  spherical  polar  coordinates 
(r,0,e)  illustrated  in  Figure  3  a  free,  monochromatic  plane  wave  of 
unit  amplitude  propagating  in  the  direction  Oo  given  by  (0o,eo)  is 

f  (r,t)  =  e  fcrtraCe^cosC  0)  +  sin(0o)cos(  8  -  e0)]  e-wrt     (43) 
[Ref.  16] 


Figure  3:  Spherical  Polar  Coordinates 
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For  a  wave  of  unit  amplitude  incident  along  the  z  axis  ,0=0 
and  the  free  wave  f (r,t)  is  presented  as 

f(r,t)  =  e  1kr  cos(0)  e"1a)t  =  e  ikz  e~m  .  (49) 

Note  this  is  a  solution  with  axial  symmetry  for  the  free  wave  equation 
[Ref.  18]. 

The  wave  equation  describing  the  behavior  of  the  undisturbed  or 
free  monochromatic  wave  }(r,t)  is 

(  v2-4^)f(r-,t)  =  0.  (50) 

cz  ar 

c  is  the  speed  of  propagation  of  the  free  wave.    If  the  wave 
function  f (r,t)  has  a  simple  harmonic  time  dependence,  i.e., 

j(r,t)=  e-totJ(D,  (51) 

then  for  a  wave  number 

k  =  o)/c  =  2TT/A  (52) 

the  free  wave  equation  becomes  the  Helmholtz  equation 

[v2  +  k2]f(r,t)  =  0  (53) 

[Ref.  21]. 


20 


The  stationary  solution  j(r)  of  the  free  wave  equation  can  be 
written  as  the  product 

f(r)-  fi(qi)f2(q2)f3(q3)  <54) 

for  orthogonal  curvilinear  coordinates  q.,  q2,  q3    [Ref.  16].  Recall  an 

undisturbed  plane  wave  has  a  cylindrical  symmetry  about  the  axis  of 
propagation  [Ref.  22].    Using  the  spherical  polar  coordinates 

q^r,  q2=0,  q3=  e,  (55) 

the  stationary  solution  describing  an  undisturbed  plane  wave  propagating 
along  the  z  axis  is 

f(r)=  )1(D{2(0)  .  (56) 

Substituting  f(r)  into  the  free  wave  equation  expressed  in 
spherical  coordinates  and  separating  the  coordinate  dependencies  yields 
the  separated  free  wave  equation, 

(kr)2+  f  dFr2iH  =  "  T-J—  T  (s1n  0 f  V-  (57) 

Using  a  separation  constant  of  1(1+ 1 )  the  equation  for  the  angular 
dependency  of  the  free  wave  becomes 

— -! —  (sinjer-^-f.)  +  1(1  +  0  =  O     (58) 

J?  sin  0  d#  dc<   * 
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Letting 

cos(0)  =  x  (59) 

in  the  angular  equation  yields 

[(1  -x2)-^7   -2x-?-  +  l<l  ♦  1)]  U  =  0.     (60) 

L  dxz  dx  * 

The  previous  equation  is  known  as  the  Legendre  differential 
equation  [Ref.  21].    The  most  useful  solutions  to  the  Legendre  differential 
equation  are  the  Legendre  polynomials,  designated  by  P^x).  They  can  be 

generated  by  the  recursion  relationships 

P0(x)=1,  (61a) 

P1(x)  =  x  (61b) 

and 

P,(x)  =  [2-1/1]  x  PM(x)  -  [1-1/1]  P,_2(x)  (61c) 

[Ref.  21].    Therefore, 

f2  ,(0)  =  P,(cos  0)  .  (62) 

From  equation  (57),  the  equation  describing  the  radial 

dependency  of  the  free  wave  is 

,2 


(_Lt   2jL  +  [k2_Ki^)]}f  .        (63) 

"  dr^      r  dr  r^         ' 
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Changing  scale  to  the  variable  s  =  kr  yields  the  differential  equation 

[^^Mi-^lHrO       (64) 

ds        s  ds  s 

[Ref.  171 

This  differential  equation  has  the  spherical  Bessel  function  j^s) 

with  amplitude  A,  as  a  solution; 

flfl(s)  =  A,j,(s)  (65) 

[Ref.  18]. 

The  stationary  solution  of  the  free  wave  equation  in  spherical 
coordinates  is  now  written  as  the  partial  wave 

fyr)  =  P,(cos  0)  [A,  j,(kr)]  .  (66) 

The  most  general  form  of  this  free  wave  equation  solution  is  a  Fourier- 
Bessel  series, 

oo 

J(r\0)=  I  A^CkrJP^cosfl)  (67) 

[Ref's.  17  and  21]. 

A  plane,  free  wave  of  unit  amplitude  propagating  in  the  direction 
0Q=  0  is  then 


e     "   2c  A,  j,(kr)  P,(cos  0)  .  (68) 
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The  norm  of  the  Legendre  polynomials,  the  asymptotic  expression 

sin  (kr-  §1) 


Jifcr) 


2_1  (69) 


i^'  kr 

and  the  process  detailed  in  Elton  [Ref.  18]  are  used  to  determine  A,; 

A,  =  i1  [21  +1]  .  (70) 

Therefore, 

e,kz=    2  i1  [21  +1]  i,(kr)  P,(cos  0)  (71) 

[Ref's.  18  and  21]. 

3.      Detect  ion  Conditions  and  Assumptions 

As  the  dimensions  of  the  incident  acoustic  wave  packet  are 
greater  than  the  size  of  the  target  it  is  expected  that  a  portion  of  the 
incident  wave  packet  continues  to  propagate  past  the  target.    The  total  or 
disturbed  wave  description  is  then  a  superposition  of  the  incident  and 
scattered  wave  descriptions.  [Ref.  17] 

The  acoustic  array  serves  as  both  transmitter  and  receiver,  or 
detector.    This  keeps  the  detector  outside  the  path  of  continued 
propagation  of  the  incident  wave  packet.  [Ref.  17]    Therefore,  only  the 
scattered  wave  description  is  of  interest  in  the  calibration  of  the  array. 
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The  placement  of  the  target  at  a  far  field  range  simultaneously 
places  the  detector  in  the  region  of  asymptotic  behavior  of  the  scattered 
wave. 

4.      The  Disturbed  Wave 

Once  an  incident  wave  described  by  the  velocity  potential  l/(r) 
enters  the  region  of  influence  of  a  scattering  targets  spherical  potential 
V(r)  the  wave  equation  becomes  the  disturbed  wave  equation 

{  vMk2-  V(r)]}   K(r)  =  0  (72) 

[Ref.  1 7].    Comparison  of  the  form  of  the  disturbed  wave  equation  to  the 
free  wave  equation  gives  an  expectation  of  a  solution  of  the  form 

//(r)=  f,(r)^2(0)  .  (73) 

Substituting  K(r)  into  the  disturbed  wave  equation  and 
separating  variables  yields 

! A.  (s1ri0_d_)K=  r2[k2-  V(r)]  +  ^^-r2^-K  .(74) 

V2  sin  jer  60  K  der }  2  L  /J     1/  dr       dr    1 

The  angular  equation  for  the  disturbed  wave  is  identical  in  form 
to  the  angular  equation  for  the  free  wave  and  independent  of  the  Neumann 
boundary  condition  d  Vlbv  =  0.    Therefore, 

l/2](0)  =  P^cos  0) .  (75) 
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Using  the  previous  substitution  of  s  =  kr  yields  the  now  modified 
Bessel  equation  as  the  radial  equation  for  the  disturbed  wave; 

(i;a^+  [1  -  1(ul);V(r)]}i/^o    (76) 

1  ds-       sds  s2  JJ    ] 

Following  the  solution  development  in  Elton  [Ref.  18]  and 
applying  the  Neumann  boundary  condition  at  r  =  a  leads  to  a  solution 
Vy  fl(r)  =  i1  [21  ♦  1  ]  [jjflcr)  -  a,"  h/ ]  >(kr)]  .       (77) 
The  relationship 

.,  j/fca) 

a,  =  i  e~l6i  sin  S.  =  — ~r\ (78) 

'  '  k  1 )  ,     ■ 

h^      (ka) 

is  defined  from  Morse  [Ref.  1 7]  and  Bowman  [Ref.  1 63.    Additionally, 

•  the  j,(kr)  are  the  spherical  Bessel  functions  of  the  first  kind, 

argument  kr,  and  order  1, 

•  h,(1)(kr)  are  the  Hankel  functions  of  argument  kr  and  order  1, 

•  the  Sj  are  the  phase  shifts  of  the  usual  partial  wave  description  for 
spherical  scattering  solutions  [Ref's.  17  and  20], 

•  the  j,'(ka)  are  the  first  derivatives  of  the  j,(kr),  taken  with  respect 
to  the  argument  kr  and  evaluated  at  r  =  a, 

•  h,(1)'(ka)  =  j,'(ka)  ♦  i  y,'(ka),  and  (79) 

•  the  y,'(ka)  are  the  first  derivatives  of  the  y,(kr),  the  spherical  Bessel 

functions  of  the  second  kind  with  argument  kr  and  order  1,  taken  with 
respect  to  the  argument  kr  and  evaluated  at  r  =  a. 
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The  general  series  solution  for  the  total  disturbed  wave  Kr)  is  then 


oo 


V(r)  =    i  i1  [21  ♦  1  ]  P,(cos  0)  [j,(kr)  -  a,*  h/ 1  )(kr)]  .  (80) 

5.      The  Scattered  Wave  and  Its  Amplitude 

The  general  asymptotic  form  for  a  disturbed  wave  Kr,0) 

resulting  from  a  plane  wave  (/.(r,0)  incident  along  00  =  0  and  scattered  by 
a  spherical  potential  V(r)  is 

l/(r,0)  =  eikz  +  s(0)eikr/r  (81) 

[Ref.  18].    This  description  of  the  asymptotic  disturbed  wave  is  a 
superposition  of  the  incident  plane  wave  e,kz  and  a  spherical,  scattered 
wave  s(0)e,kr7r.  The  term  s(0)  is  called  the  scattering  amplitude 
[Refs.  18,  19  and  22]. 

Now,  recalling  the  series  equivalence  for  the  plane  wave  incident 
along  0Q  in  equation  (71 ),  and  comparing  the  general  series  solution  for 

the  disturbed  wave  in  equation  (80)  to  its  asymptotic  form  in  equation 
(81)  gives 


s(0)eikr7r  =  -  2   i]  [21  +  1  ]  P,(cos  0)  a/  h  < '  >(kr) .  (82) 


1  =  0 
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This  is  a  series  expression  for  the  scattered  wave.    Therefore, 


s(O)  =  — —  i  2  [21  +  l]  P,(cos  0)  a.  (83) 

[Ref's.  16  and  18]. 

6.      The  Backscattered  Cross  Section 

The  scattering  amplitude  s(0)  is  related  to  the  differential 
scattering  cross  section  <r(0)  =  v\  t  by 

<r(0)  =  ls(0)|2  (84) 

[Ref's.  18,  19  and  22]. 

Applying  the  series  expression  for  the  scattering  amplitude  in 
equation  (83)  to  the  previous  equation  yields 

<rA.f  =tr(0)=  -2     12    i[21+  I]  a,  P,(cos0)|2.    (85) 

l9l  ^        1=0  '       ' 

For  the  backscattered  direction  0  =  n  this  gives  a  backscattered  cross 
section  of 

'bV-*B-  77   l1i0[2,  +  1lai"1  (86a) 


[Ref's.  16  and  17].   For  the  spherical  target  At  t  =  na2,  so 


vh=  -=-~   I  2  [21+  l]a.  I2.        (86b) 
b   (ka)z  1=0       ' 
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It  should  be  noted  a  plane  wave  of  unit  amplitude  incident  from 
the  positive  z  direction  and  backscattered  along  0  -  0  would  introduce  a 
factor  of  (- 1 )'  into  the  previous  sum  [Ref.  1 6]. 

E.       SPHERICAL  BE55EL  FUNCTIONS  AND  THEIR  DERIVATIVES 

Implementation  of  the  previous  equation  necessitates  calculation  of 

the  a '.    The  recursion  relations 

j0'(ka)  =  H0(ka)  ♦  cos(ka)]/ka  ,  (87a) 

j/Oca)  =  {-[l  +  l][j,(ka)]/ka}  ♦  j,_,(ka) .  (87b) 

y0'(ka)  =  t-y0(ka)  +  sin(ka)]/ka  (87c) 

and 

y/Cka)  =  {-[l  +  ]][y1(ka)]/ka}  +  y1_1(ka)  (87d) 

determine  the  j,'(ka)  and  y,'  (ka)  needed  for  the  calculation  of  a '  in 

equation  (78)  [Ref's.  21  and  23]. 

Determination  of  the  regular  and  irregular  spherical  Bessel  functions 
at  r  =  a  supports  the  use  of  the  recursion  relations  in  equations  (87).    As 
ka  is  real  the  irregular  spherical  Bessel  functions  are  stable  enough  to  use 
the  recursion  relations 
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y0(ka)  =  -[cos(ka)]/ka  .  (88a) 

y  ,(ka)  =  [y0(ka)  -  sin(ka)]/ka  (88b) 

and 

y,(ka)  =  {[21  -  l][y1_,(ka)]/ka}-y1_2(ka)  (88c) 

in  their  determination  [Ref's.  23  and  24]. 

Because  the  size  of  the  argument  of  the  regular  spherical  Bessel 
functions  may  be  very  small  the  method  of  negative  continued  fractions  is 
appropriate  to  the  fast  and  accurate  determination  of  a  converging 
solution  to  these  functions  [Ref.24].    The  regular  spherical  Bessel  function 
reciprocals  are  therefore  determined  by  the  negative  continued  fraction 
calculation 

i    =     i      C-lEb,  ^[b,  2jb,  ,1-Eb,  JEb,^,^  t]...    (89) 

jj(ka)      j0(ka)  [b12] ^\,2^b\,iP\,2^-  * 

Now, 

blm=    2(1  +  m-  l)/ka  (90) 
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and 


=b         ' 


1,™      b,  m   *  -  1 


l,m-1    '  J (91) 

blm-2_J _ 


bl  m-3 


b 


1,2      blJ 

[Ref.  25]. 

While  the  regular  spherical  Bessel  function  reciprocal  calculation 
of  equation  (89)  would  appear  to  involve  an  infinite  number  of  terms,  the 
calculation  can  in  fact  be  terminated  when  the  m1*  numerator  term  equals 
the  m01  denominator  term  to  the  desired  precision  for  the  spherical  Bessel 
function's  solution.    [Ref's.  24  and  25] 
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III.     CALIBRATION    PROCESS    AND  SOFTWARF 

Small  quantities  of  data  were  collected  in  the  anechoic  chamber  at 
the  Naval  Postgraduate  School.  This  data  aided  first  order  verification  of 
the  calibration  process  and  its  supporting  software.    A  brief  description 
of  the  anechoic  chamber  is  provided  in  Appendix  A. 

The  calibration  process  and  data  development  occured  in  two  stages. 
The  first  stage  used  equation  (26)  and  two  way  propagation  path 
measurements  to  determine  <rbErEtG0.    The  second  stage  determined  EtGQ 

using  equation  (17)  and  one  way  propagation  path  measurements. 

E   was  then  determined  from  <r.E  E,Gn,  E,Gn  and  the  calculated  o\  of 

r  brtO'      tO  b 

equation  (86b). 

Because  reciprocity  was  assumed  for  Ep  and  Et 

ErEt  =  Er2  (92) 

provided  the  desired  calibration  information 

Supporting  software  was  written  to  perform  calculations  required  for 

calibration    The  next  section  describes  the  development  of  this  software. 
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A.      CALIBRATION  SOFTWARE  DEVELOPMENT 

Calibration  software  was  written  in  BASIC.  It  was  used  with  a 
BASIC  4.0  compiler  on  a  Hewlett  Packard  310  computer.  Computer 
performance  was  enhanced  by  an  Infotek  floating  point  processor. 

Notation  in  the  programs  sometimes  differs  from  that  in  the  thesis 
text  as  a  result  of  H.P.  keyboard  limitations.  Variable  definitions  are 
located  at  the  beginning  of  subprograms  or  modules.  These  definitions 
are  provided  throughout  a  program  as  an  aid  in  determining  program  to  text 
variable  equivalence. 

The  main  effort  of  software  development  was  calculation  of 
backscattered  cross  sections.  This  was  done  with  equation  (86b). 
Information  available  for  summation  loop  verification  was  differential 
scattering  cross  section  information.  Differential  scattering  section 
software  was  developed  first. 

1.      Differential  Scattering  Cross  Section  Software  Development 
The  differential  scattering  cross  sections  <r  were  calculated 
using  equation  (85)  and  equation  (78). 

To  support  equation  (78)  employment  the  subprograms 
calculating  the  values  of  the  Legendre  polynomials  and  the  values  of  the 
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regular  and  irregular  spherical  Bessel  functions  were  developed  first. 
Each  subprogram  was  developed  and  tested  separately  in  test  programs. 
Each  test  program  was  designed  to  provide  a  tabular  output  for 
verification  against  values  found  in  Abramowitz  and  Stegun  [Ref.  23]. 

The  Legendre  values  were  calculated  using  equations  (61 ).    The 
irregular  spherical  Bessel  functions  were  calculated  using  equations  (88). 
The  regular  spherical  Bessel  functions  were  calculated  using  equations 
(89),  (90)  and  (91). 

Following  verification  these  subprograms  were  assembled  with 
subprograms  calculating  the  derivatives  of  the  regular  and  irregular 
spherical  Bessel  functions,  the  normalized  differential  scattering  cross 
section  <r,  the  scattering  modulus  of  Bowman  [Ref.  16],  and  with  three 
graphics  subprograms.  This  assembly  created  the  differential  scattering 
cross  section  program  found  in  Appendix  B. 

The  values  of  the  regular  and  irregular  spherical  Bessel  function 
derivatives  at  r=a  were  calculated  using  equations  (87). 

Differential  scattering  cross  sections  were  calculated  using 
equation  (85).  The  number  of  orders  1  required  for  the  summation 
calculation  in  equation  (85)  was  determined  by  generating  tabular  outputs 

34 


and  determining  the  number  of  orders  required  for  numerical  stability  to 
occur  to  12  decimal  places.  This  stability  process  was  based  on  the 
concept  of  the  phase  shift  6,  of  each  outgoing  partial  scattered  wave 

becoming  progressively  smaller  with  increasing  1  until  a^  (equation  (78)) 

became  negligible  for  the  precision  desired  [Ref's.  18,  19,  20  and  22], 
Tabular  outputs  covered  a  ka  value  range  of  one  to  thirty  at 

varying  scattering  angles,  including  0  =  tt.  This  ka  value  range  selection 

was  based  on  the  expected  maximum  target  size  to  be  encountered  during 

actual  echosounder  employment. 

Twenty-one  orders  greater  than  the  integer  value  of  the 

argument  ka  were  determined  to  be  required  for  equation  (85)'s 

summation.    The  number  of  orders  required  was  greater  than  the  usually 

cited  relationship 

ka  ~  VK1+  1)  (93) 

gives.  It  should  be  noted  that  the  precision  criteria  of  <r  was  greater  than 

that  normally  encountered.  [Ref's.  18,  19,  20  and  22] 

The  graphics  subprograms  of  the  differential  scattering 

cross  section  program  provided  outputs  for  comparisons  to  Figure  80  in 
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Morse  [Ref.  17]  and  Figure  10. 10  of  Bowman  [Ref.  16].    These  comparisons 
were  used  to  verify  the  functioning  and  output  of  the  differential 
scattering  cross  section  program. 

The  graphics  subprograms  used  to  generate  Figure  4  provide 
linear  polar  plots  of  the  differential  scattering  cross  section  <r(0) 
normalized  by  the  cross  sectional  area  At .  of  the  target  sphere.    These 

linear  polar  plots  show  the  normalized  differential  scattering  cross 
section  of  acoustically  hard  spheres  ensonified  by  plane  acoustic  waves 
incident  from  the  left,  i.e.  from  0  =  tt. 

The  first  view  of  each  sheet  in  Figure  4  shows  the  entire 
normalized  differential  scattering  cross  section  pattern,  including  the 
forward  lobe,  or  "shadow  zone",  along  0  =  0.    The  second  view  of  each 
sheet  is  a  magnified  plot  showing  the  scattering  detail  at  a  scale 
comparable  to  the  normalized  backscattered  cross  section  value.  Note 
the  two  views  are  identical  on  Sheet  1.  The  need  for  magnification  with 
increasing  ka  becomes  apparent  on  subsequent  sheets. 
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Figure  A,  Sheet  1  of  3:    Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  80  of  Morse  [Ref.  1 7]. 
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Figure  A,  Sheet  2  of  3:    Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  80  of  Morse  [Ref.  1 7]. 
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Figure  4,  Sheet  3  of  3:    Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  80  of  Morse  [Ref.  1 7], 
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Figure  4  was  used  for  comparison  to  Figure  80  of  Morse  [Ref.  17]. 
The  radial  scale  of  Morse  was  not  graduated  so  verification  consisted  of 
a  check  for  the  correct  number  and  approximate  angular  location  of 
minimums  and  maximums  in  the  differential  scattering  cross  section 
patterns. 

The  graphics  subprograms  used  to  generate  Figure  5  provide 
semi-logarithmic  plots  of  the  scattering  modulus  5  of  Bowman  [Ref.  16]. 

5  =  k  s(0)  /  V4t7  (94) 

is  the  scattering  modulus  5  of  Bowman  [Ref.  163.    s(g)  is  the  scattering 
amplitude  s(0),  as  defined  by  equation  (83),  with  a  factor  of  (-1)' 
introduced  into  the  sum,  that  is, 


s(0)  =    ^^  i  2    (- 1 )]  (21  *  1 )  P,(cos  0)  a/  .       (95) 
k        i=o  '  I 


s(0)  is  normally  used  to  calculate  <r(0),  per  equation  (84).  The  scattering 
modulus  5  of  Bowman  [Ref.  1 6]  could  also  be  used  to  calculate  <r(0).  Note 
the  introduction  of  (- 1 )'  into  the  summation  corresponds  to  the  case  of  an 
acoustically  hard  sphere  ensonified  by  a  plane  acoustic  wave  incident  from 
0  =  0. 
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Comparison  of  Figure  5  to  Figure  10. 10  of  Bowman  [Ref.  16] 
shows  a  match  in  both  0  and  5.  This  graphical  match  verifies  the 
functioning  and  output  of  the  differential  scattering  cross  section 
program  of  Appendix  B. 
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Figure  5,  Sheet  1  of  3:    Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  1 0. 1 0  of  Bowman  [Ref.  1 6]. 


41 


10 


0"J 


p 


Q 

(- 
U 
01 


60  20  120 

PHI     ( DEGREES) 


150 


80 


C    =    343    M/S. 

F    =   5000   HZ. 

A  =  .0545901454805  M. 

K*A  =  5  . 

NORMALIZED  BCKSCTR  =   1 


0812038G888 


Figure  5,  Sheet  2  of  3:    Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  1 0. 1 0  of  Bowman  [Ref .  1 6], 
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Figure  5,  Sheet  3  of  3:    Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  1 0. 1 0  of  Bowman  [Ref.  1 6]. 
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2.      Backscattered  Cross  Section  Software  Development 

The  differentia]  scattering  cross  section  program  of  Appendix  B 
was  modified  to  yield  the  backscattered  cross  section  program  of 
Appendix  C.    With  the  modification  the  normalized  differential  scattering 
cross  section  is  calculated  for  the  backscattered  direction  only.  Equation 
(86b)  is  used  for  this  calculation. 

Additionally,  the  graphics  subprograms  of  the  differential 
scattering  cross  section  program  were  replaced  by  one  graphics 
subprogram.  The  new  graphics  subprogram  was  designed  to  provide  a 
rectangular  graph  of  the  normalized  backscattered  cross  section,  <rb. 

o"b  is  plotted  as  a  function  of  ka  for  a  specifiable  range  of  values  of  ka. 

Figure  6  was  generated  by  the  backscattered  cross  section 
program.  It  was  used  for  comparison  to  Figure  10. 11  of  Bowman  [Ref.  16]. 
This  comparison  served  as  a  verification  of  the  calculated,  normalized 
backscattered  cross  section  values.  Note  the  ka  axis  is  labeled  K*A. 

Figure  7  was  generated  to  verify  the  normalized  backscattered 
cross  section  profile  was  approaching  a  value  of  one  as  ka  left  the  Mie 
scattering  regime  and  entered  the  "geometric  optics"  region.  Note  the  ka 
axis  is  labeled  K*A. 


AA 


u 

OJ 

Ifi 
(J) 
O 

q: 

(J 

h- 
(J 
en 
*: 
u 

n 

r 
ct: 
o 


10 


K*R 


Figure  6:    Backscattered  Cross  Section  Program  Output 
for  Verification  of  Calculated  Normalized  Backscattered 
Cross  Section  Values  by  Comparison  to  Figure  1 0. 1  I  of 
Bowman  [Ref.  16]. 
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Figure  7:  Backscattered  Cross  Section  Program  Output 
for  Verification  of  the  Normalized  Backscattered  Cross 
Section  Asymptotic  Behavior  with  Increasing  ka. 
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3.      Acoustic  Array  Efficiency  Software  Development 

The  acoustic  array  efficiency  program  in  Appendix  D  determined 
Er  Et.    In  general,  each  major  section  of  Chapter  1 1  of  this  thesis  and  each 

set  of  analyzed  data  is  represented  by  a  subprogram  of  the  acoustic  array 
efficiency  program. 

The  acoustic  array  efficiency  program  of  Appendix  D  includes  the 
non-graphics  subprograms  present  in  the  backscattered  cross  section 
program  of  Appendix  C.  It  also  includes  subprograms  designed  to  furnish 
stored  data,  calculate  attenuation  and  range  using  equations  (27)  through 
(42),  and  use  developed  information  to  calculate  EtGQ  with  equation  (17). 

Equation  (86b)  was  used  to  calculate  <rb. 

The  appropriate  subprogram  and  main  program  calculational 
results  were  then  used  in  equation  (26)  to  determine  E  .  E  was  used  in 

r       r 

equation  (92)  to  determine  Ep  Et. 

B.      DATA  ACQUISITION  AND  SUPPORTING  EQUIPMENT 

Data  collection  occured  28  November  1987,  30  November  1987,  and 
7  December  1987. 
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Data  collected  28  November  and  30  November  provided  information 


needed  for  calculation  of  ffbErEtG0.    Data  collected  7  December  allowed 


calculation  of  E.Gn  and  the  calculation  of  E  's. 

tO  r 


Collected  data  is  presented  in  Appendix  E. 


Two  Way  Propagation  Path  Measurements 


Two  way  propagation  path  data  were  collected  28  and  30 
November.    Measurements  were  used  to  calculate  the  product  trbErEtG0. 


Figure  8  provides  a  connection  schematic  for  the  data  acquisiton 


equipment. 
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Figure  8:    Data  Acquisition  Equipment  Set-up 
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Four  target  spheres  were  available  and  measured  for  size.    They 
included  one  hollow  aluminum  sphere  (sphere  number  one)  and  three  solid 
aluminum  and  brass  spheres.    Table  ( 1 )  lists  the  data  and  the  calculated 


values  of  At  t  =  Tta2,  ka,  and  <rh.    <rh  was  calculated  from  equation  (86b). 


b     "b 


Table  1:  Target  Sphere  Dimensional  Data  and  Normalized 
Backscattered  Cross  Section  Calculational  Results  for 
Frequency  of  5000  [Hz]. 


Target      Sphere  Sphere 

No.  Diameter  [m]       Radius  [m] 


\t  ^ 


ka 


1 

2 
3 

4 


.2546  ±0006        .1273 


1013  ±0001        .05065 


.07634±.00003     .03817 


.06240±.00008     .03120 


.0509 


.004577 


.62 


.003058  2.843 


.9142 


.008060  4.619        1.000 


3.479        1.051 


.6836 


Figure  9  shows  the  target  spheres'  theoretical  normalized 
backscattered  cross  section  distribution  plotted  on  a  graph  generated  by 
the  backscattered  cross  section  program  of  Appendix  C. 
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Figure  9:    Target  Spheres'  Normalized  Backscattered  Cross 
Section  Distribution  as  Calculated  from  Measured  Diameters, 
Using  an  Assumption  of  Perfect  Acoustic  Hardness. 


The  target  spheres  were  suspended  in  turn  in  the  anechoic 
chamber.    Placement  was  made  using  the  criteria  in  Section  D.I  of 
Chapter  II.    The  range  R  from  the  array  to  a  target  was  determined  to  be 


5.30  [m]. 


Temperature  Tc  and  relative  humidity  Rh  along  the  propagation 


path  were  determined  with  a  Weathermeasure  Temperature  and  Relative 
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Humidity  Meter.    Tc  was  usually  about  21  [°c]  with  a  relative  humidity  of 

approximately  50%.    Atmospheric  pressure  P  was  read  from  a  barogram 

supplied  by  the  Naval  Postgraduate  School  meteorological  station. 
Specific  T  ,  Rh,  and  P  data  are  furnished  in  Appendix  E. 

The  unshrouded,  close  packed,  hexagonal  acoustic  array  of 
Moxcey  [Ref.  8]  was  used  to  ensonify  each  target. 

A  =  Ns  tt  rs2  (96) 

determined  the  aperture  area  A  of  the  array.    Ns  =  1 9  was  the  number  of 

speakers  in  the  array  and  rs  =  .0381  [m]  was  the  average  speaker  horn 

radius  [Ref.  8]. 

Moxcey  [Ref.  8]  determined  the  half  width  of  the  array's  main 
lobe  to  be  12°.    The  angular  half  widths  of  targets  two,  three  and  four 
were  all  less  than  5%  of  the  main  lobe's  half  width.  Target  one's  angular 
half  width  was  1 1.5%  of  the  main  lobe's  half  width. 

A  sinusoidal  input  from  a  Hewlett  Packard  3314  A  function 
generator  drove  the  array  via  a  pre-amplifier.    Signal  input  from  the 
HP3314A  function  generator  to  the  pre-amplifier  was  set  for  20  cycles  at 
5000  [Hz].    The  voltage  input  to  the  pre-amplifier  s  power  amplifier 


terminal  (see  Figure  8)  was  3.5  [V]  for  spheres  two,  three  and  four  and  3.0 
[V]  for  sphere  one. 

Twenty  cycles  is  considerably  less  than  the  number  of  cycles 
expected  in  actual  echosounder  application.  It  was  necessary  to  keep  the 
pulse  packet  this  small  because  the  anechoic  chamber's  corner  was  only 
about  1.5  [m]  from  the  target  mounting  location.    Voltage  input  for  sphere 
one  was  reduced  to  3.0  [V]  to  avoid  clipping  of  its  greater  return  signal 
voltage. 

Actual  transmission  voltage  inputs  Vt  to  the  array  from  the 

pre-amplifier  were  measured  with  a  Nicolet  3091  oscilloscope  at  the 
pre-amplifier's  array  terminal.  The  array  was  disconnected  from  the 
pre-amplifier  for  V.  measurements.    Vt  values  are  provided  in  Appendix  E. 

They  are  also  stored  in  the  Trans30  and  Trans35  subprograms  of  the 
acoustic  array  efficiency  program  in  Appendix  D. 

The  Nicolet  oscilloscope  was  connected  to  a  Hewlett  Packard 
7090  A  measurement  plotting  system.  This  connection  provided  a  hard 
copy  plot  of  voltage  vs.  time  traces. 

The  returned  signal  was  amplified  and  filtered. 
Pre-amplifier  gain  Ge  was  determined  asl  1,094  by  Moxcey  [Ref.  8], 
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Filtering  was  done  with  a  Wavetek  Brickwall  filter  model  753  A  to  remove 
the  low  frequency  normal  modes  of  oscillation  of  the  speaker  drivers  from 
the  return  signal.  Filtering  was  accomplished  at  0  [dB]  gain  and  5000  [Hz] 
for  both  low  and  high  pass  settings. 

The  processed  return  signal  voltages  Vr  were  displayed  and 

measured  on  the  Nicolet  3091  oscilloscope.  Vp  measurements  were  made 

at  the  Out  2  terminal  of  the  Wavetek  Brickwall  filter  (see  Figure  8).  They 
were  taken  at  10  [jjlS]  intervals  for  two  full  cycles  along  the  asymptotic 
region  of  the  return  signal  trace. 

The  Nicolet  3091  oscilloscope  was  connected  to  the  Hewlett 
Packard  7090  A  measurement  plotting  system.  Figure  10  shows  the 
Hewlett  Packard  plotter's  traces  for  target  sphere  two's  return  signal 
voltages  as  printed  and  subsequently  labeled.    The  figure  is  representative 
of  the  oscilloscope  traces  for  each  target  return.    Data  for  all  targets  are 
provided  in  Appendix  E. 

The  time  of  return  t  was  determined  as  the  time  difference 

r 

between  the  first  positive  maxima  of  the  transmitted  pulse  packet  and  the 
first  positive  maxima  of  the  target's  return  pulse  packet. 
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Figure  10,  Sheet  1  of  2:    Target  Sphere 
Two's  5000  [Hz]  Return  Signal  Trace 
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Figure  10,  Sheet  2  of  2:    Target  Sphere  Two's 
5000  [Hz]  Return  Signal  Trace,  Scaled  Down 
For  Oscilloscope  Measurement. 
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Noise  voltage  Vn  measurements  were  made  with  the  same 
measuring  equipment  setup  used  for  V  measurements.    The  target  spheres 
were  removed.    V  measurements  were  taken  at  a  time  of  return  t 

n  r 

comparable  to  the  time  of  measurement  of  V  . 

Vn  measurements  were  actually  made  on  occasions  well 

separated  chronologically  from  the  V  measurements.    Therefore,  they 

were  proportionally  corrected  in  line  2530  of  the  acoustic  array 
efficiency  program.    The  corrections  account  for  the  slight  variations  in 
input  voltages  and  propagation  path  attenuations  that  occurred  between 
the  two  measurement  occasions. 

Figure  1 1  shows  a  representative  noise  trace.    Data  from  the 
noise  measurements  is  presented  in  Appendix  E  and  stored  in  the  Noise30 
and  Noise35  subprograms  of  the  acoustic  array  efficiency  program  in 
Appendix  D. 

Unless  otherwise  specified  above  the  remaining  data  from  the 
two  way  propagation  path  measurements  are  stored  in  the  subprograms 
titled  Sphere  1,  5phere2,  5phere3  and  5phere4    All  remaining  data  are  also 
presented  in  Appendix  E. 


55 


V. 


^pwoat*^ 


400  tjiS 


33.80  [mS] 


34  19  [m5] 


Figure  1 1:  Representative  5000  [Hz]  Noise  Trace. 
3.0  [V]  Supplied  by  the  HP3314A  Function  Generator. 


2.      One  Way  Propagation  Path  Measurements 

One  way  propagation  path  data  were  collected  7  December. 
Measurements  were  used  to  calculate  EtGo  with  equation  ( 1 7). 

Data  from  the  one  way  propagation  path  measurements  are 
provided  in  Appendix  E.     Data  are  also  stored  in  the  Gain30  and  Gain35 
subprograms  of  the  acoustic  array  efficiency  program  of  Appendix  D. 
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For  one  way  path  measurements,  i.e.  for  the  determination  of 
EtG0,  an  Ivie  Electronics  Inc.  IE  30  A  Audio  Analyzer  (Serial  *805  A  426) 

was  used  to  measure  received  target  intensity  I   in  [dBJ.    Intensity 

measurements  were  referenced  to  10~12  [W/m2].    The  acoustic  signal 
supplied  for  this  measurement  was  a  continuous  wave  at  5000  [Hz]. 
Measurements  were  taken  at  the  suspension  site  of  the  calibrating  targets. 

A  Fluke  8060  A  True  RMS  Multimeter  was  used  to  measure  the 
root-mean-square  voltage  (V     ).  supplied  to  the  acoustic  array  at  its 

input.  (V     ).  was  squared  to  determine  (V    )t 

The  Weathermeasure  temperature  and  relative  humidity  meter 
was  unavailable  for  Et60  data  collection.    Hygrothermographs  from  the 

Naval  Postgraduate  School's  meteorological  station  were  used  to  estimate 
the  anechoic  chambers  relative  humidity  Rh  and  temperature  Tc  on 

7  December.    Estimates  were  based  on  comparisons  to  the 
hygrothermographs  and  measured  Rh  and  T   of  28  and  30  November. 


57 


C.      RESULTS 

Data  analysis  was  performed  with  the  aid  of  the  acoustic  array 
efficiency  program  of  Appendix  D.  Data  used  in  calculations  is  furnished 
in  Appendix  E. 

The  product  <rbErELGo  was  generated  from  two  way  popagation  path 

data  and  application  of  equation  (26).    The  average  of  calculated  <rb's 
compared  to  the  average  of  the  experimental  product  Q"bErEtG0  gives  a 
normalizing  product  ErEtGo.    The  experimental  product  <rbErEtG0  and  the 
normalizing  product  ErEtGo  are  used  to  calculate  experimental  <rb's. 
These  results  are  supplied  in  Table  2. 

Table  2:  Target  Sphere  Experimental,  Normalized 
Backscattered  Cross  Section  Calculational  Results. 

Target       ka  [m2]      Calc'd  ffbErEtGo       ErEt6o         Exp "' 

No.  <rb  <rb 

1  11.62  .9142  55.90  .947 

2  4.619  1.000  59.88  1.01 

3  3.479  1.051  61.62  1.04 

4  2.843          .6836  37.98  .644 
Avg.                              .912  53.85  59.0 
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The  experimental  <r.  are  plotted  on  a  graph  generated  by  the 


backscattered  cross  section  program  of  Appendix  C.    This  plot  is  Figure 
12.  It  is  similar  to  Figure  9  and  gives  or.  as  a  function  of  ka.  Note  ka  is 


labeled  K*A. 
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Figure  12:    Target  Spheres*  Experimental,  Normalized 
Backscattered  Cross  Section  Distribution  Compared  to 
a  Calculated  Backscattered  Cross  Section  Distribution 
Curve  Generated  by  the  Backscattered  Cross  Section 
Program  of  Appendix  C. 
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As  ka  increases  the  experimental  <rb  distribution  rises  from  just 
below  the  calculated  <rb  curve  to  just  above  it. 

EtGo  was  determined  using  one  way  propagation  path  data  and  equation 

( 1 7).    EtGQj  or  Et,  was  assumed  to  be  the  same  for  spheres  two,  three,  and 

four  because  the  supply  voltage  and  target  locations  were  the  same.    The 
angular  half  widths  of  target  spheres  two,  three  and  four  were  between 
4.6%  and  2.8%  of  the  echosounder*s  main  lobe's  half  width. 

Table  3:  Acoustic  Array  Product  EtGo 

Determined  from  One  Way  Propagation 
Path  Measurements  and  Equation  ( 1 7). 

Target      HP3314A  EtGQ 

No.  Func'n  Gen'r 

Supply  Volt.  [V] 

1  3.0  94.70 

2  3.5  95.11 

3  3.5  95.11 

4  3.5  95.11 

Calculated  <rb,  EtGo  and  two  way  propagation  path  data  were  used  in 
equation  (26)  to  determine  Er.    crb  values  were  calculated  from  the  target 
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spheres"  measured  diameters,  2a.    At  .  was  calculated  as  rra2.    Calculation 
results  are  furnished  in  Table  4. 

Table  4.     E  and  EtEr    Er  Determined  From  Two  Way  Propagation 
Path  Data  and  Equation  (26).  EtEp  Determined  From  Assumption 
of  Efficiency  Reciprocity. 


Target      Al .  [m2]  Calc'd 

No.  ov 


ffb  Atgt  [m2]  Er 


EtEr 


1 


.0509 


.9142 


.04654 


.6456    .4169 


.008060 


.000 


.008060    .6296    .3964 


.004577 


.05 


.004810  .6167  .3803 


.003058 


.6836 


.002090  .5841  .341  1 


The  E  "s  determined  for  Moxcey's  unshrouded  hexagonal  array  are  in 


the  range  of  .58  to  .65.    This  is  above  the  value  range  of  0.5  ±  0. 1 
determined  by  Weingartner  for  his  square  acoustic  array  [Ref.  6]. 

Table  4  results  were  used  to  construct  Figures  1 3  and  1 4.    Figure  1 3 
gives  Er  as  a  function  of  the  calculated,  normalized  backscattered  cross 

section,  g\.  Figure  14  gives  E   as  a  function  of  the  apparent  target  size, 


ffb\gt 
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Figure  13:  Acoustical  to  Electrical  Power  Conversion  Efficiencies, 
E_  as  a  Function  of  the  Normalized  Backscattered  Cross  Sections, 


r. 


<r. ,  Calculated  from  Measured  Sphere  Sizes. 
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Figure  14    Acoustical  to  Electrical  Power  Conversion  Efficiencies, 
Er  as  a  Function  of  the  Apparent  Target  Sizes,  <rbAt .. 
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Recalling  equation  (20), 


A 


Ptgt  =  [|][e-«R3(0)][^] 


R 


(20) 


and  equation  (12), 


P 


tgi 


ra 


4n'R 


-OCR     . 

-    <jh    e  A 

2      »  / 


(12) 


serves  as  a  reminder  that  the  received  acoustic  power,  Pra,  depends  on  the 
apparent  target  size,  tfh\Ql 

Figure  1 4s  trend  and  the  dependency  of  Ppa  on  <rbAt .  were  used  as 
suggestions  for  the  construction  of  Figure  15.  Figure  15  is  a  plot  of  E 
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Figure  15:    Acoustical  to  Electrical  Power  Conversion  Efficiencies, 
E    as  a  Function  of  the  Received  Acoustic  Powers.  P, 


ra' 
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versus  Prn.     Prfl  was  calculated  using  the  definition  (from  equation  (21 )) 


ra'        ra 


Pra=Pr/Er  (97) 

and  the  relationship  (from  equation  ( 1 5)) 

Pr  ■  (Vm5)r  /  Z  .  (98) 

Note  Figure  15"s  trend  is  for  Er  to  increase  with  an  increase  in  PpQ  for 
the  range  of  target  sizes  investigated. 
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IV.     rONCI  IJ510N5    AND  RFr.OMMFNDATIONS 

A.      CONCLUSIONS 

The  larger  magnitudes  of  Er  determined  for  Moxcey's  hexagonal  array 

suggest  it  may  be  more  efficient  than  Weingartner's  square  array. 
However,  the  drivers  for  each  array  are  identical. 

The  0.08  to  0.15  differences  in  efficiencies  are  more  likely  the  result 
of  the  difference  in  calibration  methods.  Weingartner's  measurements 
[Ref.  6]  used  steady  state,  continuous  waves.    The  use  of  steady  state, 
continuous  waves  allows  the  development  of  standing  waves.    This 
calibration  investigation  used  pulsed  waves.  The  use  of  pulsed  waves  did 
not  provide  sufficient  time  for  standing  waves  to  completely  develop  in 
the  backscattered  energy.  [Ref.  26]  This  calibration  method  should  provide 
more  accurate  efficiency  values  since  it  more  nearly  approximates  the 
echosounder  system's  actual  field  operating  configuration. 

Additionally,  steel  braid  and  nylon  lines  were  used  to  suspend  the 
target  spheres.  These  supports  probably  reflected  enough  acoustic  energy 
to  have  a  small  influence  on  results.  Other  fixtures  in  the  anechoic 
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chamber,  such  as  support  clamps  and  lights,  may  have  also  contributed 
small  amounts  of  reflected  energy  to  the  returned  signal. 

Calibration  of  the  computer  controlled  echosounder  required 
determination  of  the  product  ErEt.    The  calibration  process  assumed 

Er  =  EL.    The  calibration  results  in  Tables  3  and  4  for  spheres  two,  three, 

and  four  indicate  Er  varied  with  constant  EtG0,  or  constant  EL.  Such  a 

variation  negates  the  assumption  of  reciprocity  for  Ep  and  Et. 

Figure  15  suggests  the  acoustic  array's  efficiency  E   increased  in 

response  to  an  increase  in  P    and  (V     )  . 

ra  rms  r 

The  acoustic  gain  for  spheres  two,  three  and  four  was  probably  not 
constant,  as  assumed.    However,  G0  is  the  maximum  for  G(O),  and  Er  is 

inversely  proportional  to  GQ,  per  equation  (26).    An  averaging  process 

accounting  for  the  slight  differences  in  gain  encountered  by  spheres  two, 
three  and  four  would  have  increased  the  data  points  spread  observed  in 
Figures  14  and  15. 

Additionally,  the  pre-amplifier's  gain  may  not  have  been  constant  for 
the  range  of  received  voltage  values  encountered.  This  would  have 
introduced  either  an  apparent  increase  in  Er  or  dampened  an  increase  in  Er 
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with  an  increase  in  (V      ).  Such  a  qain  variation  could  have  effectively 

rms  r  a  ' 

negated  the  assumption  of  efficiency  reciprocity.    However,  this 

pre-amplifier  is  approximately  0.1%  linear  for  less  than  10  [V]  peak 

[Ref.  26].  Received  voltages  in  this  investigation  were  all  less  than  10  [VI 

B.      RECOMMENDATIONS 

Because  efficiency  reciprocity  is  now  in  question  it  is  necessary  to 
determine  GQ  by  calculation  or  measurement.  Such  a  determination  would 

allow  the  desired  calibration  quantity  ErEt  to  be  written  as 


EtG0 
Er  Et  =  Er    "if   .  <*» 


One  method  of  calculating  60  would  involve  integrating  the  acoustic 

array's  calculated  or  measured  intensity  pattern  (see  Moxcey  [Ref.  8])  for 
determination  of  an  equivalent  isotropic  emitter.  To  determine  G0  the 

centroid  intensity  of  the  array's  emitted  intensity  pattern  would  be 
compared  to  that  of  the  equivalent  isotropic  emitter.    This  would  be 
repeated  for  each  transmission  power  level  of  interest. 
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Calculation  of  GQ  would  allow  a  determination  of  Et  using  E,G0  .  Such 
a  determination  of  Et  would  allow  a  reciprocity  comparison  to  E  and 
support  a  more  accurate  evaluation  of  ErEt. 

The  accuracy  of  Er  and  EtGQ  could  be  further  refined  as  follows. 

Use  a  calibrating  site  that  allows  the  use  of  a  pulse  train  closer  in 
length  to  actual  field  use. 

Make  target  return  and  noise  voltages  and  site  intensity 
measurements  as  chronologically  close  as  possible. 

Measure  the  barometric  pressure  in  the  propagation  path  directly 
along  with  the  relative  humidity  and  temperature. 

The  received  and  transmitted  electrical  powers  for  Er  and  the 

transmitted  electrical  power  for  EtGQ  were  all  measured  at  the 

pre-amplifier  end  of  the  pre-amplifier/array  transmission  line.  This 
included  transmission  line  power  losses  due  to  impedance  as  efficiency 
decreases.    The  present  calibration  process  includes  the  transmission  line 
impedance  as  part  of  the  array  impedance. 

The  amount  this  impedance  power  loss  actually  affects  Er  could  be 

evaluated  by  determining  the  impedances  seen  in  each  direction  at  the 
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pre-3mplifier's  array  terminal.  This  would  be  an  opportune  time  to 
verify  the  validity  of  transmission  line  impedance  reciprocity. 

Use  the  Nicolet  oscilloscope  to  record  a  triggered  readout  from  a 
calibrated  microphone  suspended  at  the  target  site  at  45°  to  the  centroid's 
axis.    This  would  allow  EtG0  to  be  determined  for  a  pulse  packet  rather 

than  for  a  continuous  wave.    Coupled  with  a  proper  calibrating  site 
selection  this  technique  should  minimize  interference  problems  such  as 
anechoic  chamber  corner  reflections. 

A  larger  number  of  apparent  target  sizes  spanning  the  desired  range 

of  apparent  target  sizes  should  be  used  to  redetermine  Er  for  Moxcey's 

array.  The  Er  distribution  should  be  examined  for  a  smooth  and  asymptotic 

behavior.  Comparative  Ep  values  for  Moxcey's  array  could  be  determined 

using  the  method  outlined  by  Weingartner  [Ref.  6]. 

Weingartner's  square  array's  Ep  could  also  be  determined  with  this 

thesis'  calibration  process  and  compared  to  Weingartner's  E  value  [Ref.  6]. 

The  refined  Er  curve  in  the  now  calibrated  echosounder  system  should 

be  used  to  develop  an  absolute  CT2  plot  of  an  air  mass.  An  established 

method  such  as  a  tower  mounted  vertical  array  of  thermocouples  would  be 
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used  simultaneously  to  generate  a  comparative  absolute  CT2  plot  of  the 
same  air  mass. 

C.      SUMMARY 

Realizing  the  potential  of  the  computer  controlled  echosounder  for 
analyzing  lower  atmospheric  turbulence  required  a  calibration  of  the 
echosounder.    This  thesis  described  theory  and  software  for  performing 
calculations  crucial  to  a  backscattered  cross  section  calibration  of  the 
echosounder.    The  backscattered  cross  section  calibration  was  evaluated 
in  this  thesis. 

E  's  values  indicate  the  evaluated  calibration  process  requires  some 

refinement.  They  also  indicate  the  calibration  method  possesses 
sufficient  merit  to  warrant  further  development  and  that  the  software 
performs  as  intended. 
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APPENDIX    A 
ANECHOIC   CHAMBER   DESCRIPTION 

The  anechoic  chamber  is  located  in  room  019  of  building  232  at  the 
Naval  Postgraduate  School.    Acoustic  research  on  sound  sources,  sound 
receivers  and  sound  scatterers  may  be  conducted  in  the  chamber  with  a 
minimum  of  interference  from  wall  reflection  and  external  noise. 

Wall  reflection  is  minimized  by  102  centimeter  deep  wedges  made  of 
P.F.  612  fiberglass.    The  wedges  are  attached  to  the  walls,  ceiling  and 
floor  of  the  chamber.    Approximately  142  cubic  meters  of  fiberglass  is 
used  to  absorb  99%  of  incident  sound  energy  for  frequencies  greater  than 
100  Hertz.    External  noise  isolation  is  accomplished  by  separating  the 
chambers  inner  concrete  block  sides  from  the  outer  12  inch  concrete 
walls  with  a  two  inch  thick  lining  of  fiberglass  and  cork. 

A  usable  region  of  approximately  8.2  meters  by  4.3  meters  by  3.4 
meters  is  available.    This  region  is  floored  by  a  grid  of  225  wire  cables, 
with  each  cable  placed  in  a  tension  of  150  to  200  pounds  per  square  inch.  * 


*This  Appendix  was  compiled  from  the  description  posted  by  the 
entrance  to  the  anechoic  chamber  instrumentation  and  control  room. 
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APPFNDIX    B 


DIFFERENTIAL    SCATTERING    CROSS  SECTION   PROGRAM 


Examples  of  the  output  of  this  program  are  Figures  4  and  5,  located  in 
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I  ^^iji^^^ilii^5l^i^***i**^i**^^  +  +  +'**T'i*T^ii 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

ACOUSTICALLY  HARD  SPHERE  SCATTERING  OF  PLANE  ACOUSTIC  WAVES 


REM 
REM 
REM 

REM      THIS  PROGRAM  CALCULATES  AND  GRAPHS  THE  DIFFERENTIAL  SCATTERING 
REM    CROSS  SECTION  AND/OR  SCATTERING  MODULUS  FOR  AN  ACOUSTICALLY  HARD 
REM    SPHERE  OF  DIAMETER  D  ENSONIFIED  BY  AN  INCIDENT  PLANE  WAVE  OF 
REM    FREQUENCY  F. 
REM 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
| 
I++++  VARIABLE  DECLARATIONS  AND  DEFINITIONS  ++++++++++++++++++++++++++++ 


REAL  C 
REAL  Freq 
REAL  A 
REAL  L 


REAL    Lroax 


REAL    K 
REAL    Ka 


SPEED  OF  SOUND. 

FREQUENCY  OF  THE  EMITTER.' 

RADIUS  OF  THE  SPHERE. 

ORDER  OF  THE  FUNCTION  CALCULATED 
IN  A  SUBROUTINE.   IT  IS  USED 
AS  A  LOOP  INDEX  AND  TO  DESIGNATE 
THE  ARRAY  ELEMENT  FOR  THE  FUNCTION 
OF  CORRESPONDING  ORDER. 

MAX.  SIGNIFICANT  ORDER  OF  THE  SUM 
DETERMINING  THE  DIFFERENTIAL 
SCATTERING  CROSS  SECTION  (DSCS) 
AND  SCATTERING  MODULUS  (SM). 

THE  WAVE  NUMBER,  2»PI«FREQ/C. 

K*A,  THE  ARGUMENT  OF  THE  REGULAR 
AND  IRREGULAR  SPHERICAL  BESSEL 
FUNCTIONS  AND  THEIR  DERIVATIVES. 
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410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 


REAL  J<  5 1  ) 


REAL  Dj (51  ) 


REAL  Y<  51  ) 


I  ARRAY  FOR  THE  REGULAR  SPHERICAL 
BESSEL  FUNCTIONS  (RSB)  OF 
ARGUMENT  KA  AND  ORDER  0  THROUGH 
LMAX. 

ARRAY  FOR  THE  DERIVATIVES  OF  THE 
RSB  (THE  DRSB). 

ARRAY  FOR  THE  IRREGULAR  SPHERICAL 
BESSEL  FUNCTIONS  (  ISB)  OF 
AR6UMENT  KA  AND  ORDER  0  THROUGH 
LMAX. 

ARRAY  FOR  THE  DERIVATIVES  OF  THE 
ISB  (THE  DISB). 

ARRAY  FOR  THE  LEGENDRE 
POLYNOMIALS  OF  ARGUMENT  COS(PHI). 

THE  ARRAY  REPRESENTING  THE 
SCATTERING  MODULUS  (SM). 

THE  ARRAY  REPRESENTING  THE 

DIFFERENTIAL  SCATTERING  CROSS 

SECTION  NORMALIZED  BY  THE  CROSS 

SECTIONAL  AREA  OF  THE  SPHERE 

(DSCS). 

NOTE  DS(PHI  )  =  4»SM(PHI  )'2/KAA2. 

+  +  +  +  NOTE  ++++ 

THE  ARRAY  DIMENSIONS  (30+21)  NEED 

TO  BE  INCREASED  TO  ACCOMODATE 

KA>30. 


++++  MAIN  PROGRAM  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 


REAL  Dy(51  ) 


REAL  P(51  ) 


REAL  Sm<360> 


REAL  Ds(360) 


CALL  IniUC.Freq  ,A,Lmax  ,Ka  ) 


CALL  Dif sctr( Lmax ,Ka .Dj ( *  )  , J( « >  ,Dy( * 


CALL  Outpt(C  ,Freq  ,A  ,Ka  ,Sm( ♦  )  ,0s ( *  )  ) 


INPUT  AND  CALCULATE  REQUIRED 
PARAMETERS. 

,Y(  * )  ,P(  » ) ,Sm( ♦  )  ,Ds< *  )  ) 
CALCULATE  THE  SM( ♦  )  AND  DS<«). 

OUTPUT  THE  SM(*)  AND/OR  DS<«) 


890 
900 


!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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910  !++++  MAIN  PROGRAM  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++ 

920  ! 

930  PRINTER  IS  1 

940  PRINT 

950  PRINT  "END"  !  ADVISES  THE  USER  OF  THE  MAIN 

9B0  PRINT  i  PROGRAM'S  CONCLUSION. 

970  END 

980  ! 

990  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

1000  ! 

1010  !$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

1020  ! 
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1030 
1040 

1050 

1060 

1070 

1080 

1090 

1  100 

1  1  10 

1  120 

1  130 

1  140 

1  150 

1  160 

I  170 

1  180 

1  190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 


!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Init(C ,Freq  ,A  ,Lmax  ,Ka  ) 

REM    +++++++++++++++-+++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  REQUESTS  THE  APPROPRIATE  INPUT  VARIABLES  TO  DETERMINE 

REM    C,  FREQ,  A,  LMAX  ,  AND  KA. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Temp 
REAL  Circ 
REAL  Diam 


AMBIENT  AIR  TEMPERATURE  IN  DEG  C, 
CIRCUMFERENCE  OF  THE  TGT.  SPHERE 
DIAMETER  OF  THE  TARGET  SPHERE. 


!  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +++  +  +++  +  +  +  +  +  +  +  +  +  +  +  +  +++++++++++++++++++++++++  + 

I 

!++++  REQUESTING  AND  CALCULATING  SYSTEM  PARAMETERS  ++++++++++++++++++++++ 


PRINTER  IS  1 

Choice=0 

WHILE  ChoiceOI  AND  Choice<>2 

PRINT 

PRINT  "DO  YOU  DESIRE. TO  INPUT" 

PRINT  "      ( 1  )  THE  SPEED  OF  "i 

PRINT  "SOUND," 

PRINT  "   OR  (2)  THE  AMBIENT  AIR 

PRINT  "TEMPERATURE?" 

PRINT 

INPUT  Choice 
END  WHILE 

IF  Choice=1  THEN 

PRINT  "PLEASE  ENTER  C  IN  M/S." 

INPUT  C 
END  IF 

IF  Choice=2  THEN 

PRINT  "PLEASE  ENTER  TEMP  IN  "; 

PRINT  "DEGREES  CELSIUS. " 

INPUT  Tenp 

Temp=Temp+273 . 15 

C=20.05«SQR(Tenp ) 
END  IF 


INPUTS  ARE  REQUESTED  ON  THE  CRT, 


**»  DETERMINING  C  *** 


DIRECTLY  ENTER  C  , 


! 

!  OR  CALCULATE  C  FROM  TEMP. 


PRINT 

PRINT  "PLEASE  ENTER  THE  ECHOSOUNDER" ; 

PRINT  "  FREQUENCY  IN  HZ. " 

PRINT 

INPUT  Freq 


I  »«*  DETERMINING  FREQ 
!  DIRECTLY  ENTER  FREQ. 
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1580 
1590 

1600 

16  10 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

1  900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 


l 

i  ♦«*  DETERMINING  KA  ••• 
AND  Cholce<>2  AND  Choice<>3 


DO  YOU  DESIRE  TO  INPUT  " 
(  1  )THE  SPHERICAL  "; 

BESSEL  ARGUMENT  K*A," 

(2  )  THE  SPHERE'S  "j 

CIRCUMFERENCE," 

OR  (3  )  THE  SPHERE'S  " i 

DIAMETER?" 


Choice=0 
WHILE  Choices' 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

INPUT  Choice 
END  WHILE 

IF  Choice=1  THEN 

PRINT  "PLEASE  ENTER  K*A. " 

INPUT  Ka 
END  IF 

IF  Choice=2  THEN 

PRINT  "PLEASE  ENTER  THE  TARGET 

PRINT  "CIRCUMFERENCE  IN  CM." 

INPUT  Circ 

Ka=Circ*(  10"( -2  )  )«Freq/C 
END  IF 

IF  Choice=3  THEN 

FRINT  "PLEASE  ENTER  THE  TARGET 

PRINT  "DIAMETER  IN  CM. " 

INPUT  Diam 

Ka=Diam*(  1 0" ( -2  )  )*P I *Freq/C 
END  IF 


K=2*PI*Freq/C 
A=Ka/K 
Lmax=INT(Ka+21  ) 


!  DIRECTLY  ENTER  KA 


!  CALCULATE  KA  FROM  THE 
!  TARGET  CIRCUMFERENCE, 


OR  CALCULATE  KA  FROM  THE 
TAR6ET  DIAMETER. 


***    DETERMINING  REMAINING  •*« 
***  SYSTEM  PARAMETERS     *** 
CALCULATE  THE  WAUE  NUMBER. 
CALCULATE  THE  SPHERE  RADIUS. 
MAX.  SUMMATION  ORDER  CALCULATED, 


I 

!++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

i 
SUBEND 

! 

+++++++++++++4+++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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2090 
2100 
21  10 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 


!$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$ 
SUB  Difsctr(Lmax  ,Ka  ,D j ( *  )  ,  J( «  > ,Dy( « ) . Y( *  > ,P( * ) ,Sm< « ) ,Ds( * ) > 
REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM    THIS  MODULE  CALCULATES  THE  ARRAYS  SM*)  AND  DSC)  REPRESENTING  THE 
REM    SCATTERING  MODULUS  AND  THE  DIFFERENTIAL  SCATTERING  CROSS  SECTION 
REM    NORMALIZED  EY  THE  CROSS  SECTIONAL  AREA  OF  THE  SPHERE,  RESPECTIVELY, 
REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++t+++++++ 
! 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 

I 

REAL  Fac  !  FACTOR  COMMON  TO  RSM  AND  ISM. 

REAL  Rsm  !  THE  REAL  COMPONENT  OF  THE 

SCATTERING  MODULUS,  SM. 

REAL  Ism  !  THE  IMAGINARY  COMPONENT  OF  THE 

SM. 

FACTOR  COMMON  TO  SM(PHI)  8,  DS(PHI), 

REAL  Phi  !  THE  (POLAR)  ANGLE  OFF  THE  AXIS 

OF  PROPAGATION  OF  THE  INCIDENT 
PLANE  WAVE,  WITH  ORIGIN  AT  THE 
SPHERE  CENTER.   COS(PHI)  IS  THE 
AR6UMENT  OF  THE  LEGENDRE 
POLYNOMIALS  P(*)  OF  ORDER  0 
THROUGH  LMAX. 


I 

!++++  CALCULATING  THE  SM  AND  DSCS  ARRAYS  +++++++++++++++++++++++++++++++ 


CALL  DrsbtLmax  ,Ka  ,Dj(*  )  ,J(«  )  ) 

CALL   Disbdmax  ,Ka,Dy(*  ),Y(*  >) 

FOR    Phl=0    TO    360 

CALL  Leg(Lmax  ,Ka  ,Phi  ,P( *  )  ) 


CALCULATE  THE  DERIVATIVES  OF  THE 
REGULAR  SPHERICAL  BESSEL  FUNCTIONS 
OF  ORDER  0  THROUGH  LMAX. 
CALCULATE  THE  DERIVATIVES  OF  THE 
IRREGULAR  SPHERICAL  BESSEL  FUNC'NS 
!  OF  ORDER  0  THROUGH  LMAX. 


CALCULATE  THE  ARRAY  P(* ) 


Rsn  =  0 
lsm=0 

FOR  L=0  TO  Lnax 

Fac  =  (Dj(L)»(2*L+1)*P(L))/(Dj(L)*Dj(L)  +  Dy(L)*Dy(L>) 

Rsm=Rsm+< Fac*Dy( L  )  ) 

l5m=Ism+(Fac«Dj(L  )  ) 
NEXT  L 


CALCULATING  THE  REAL  AND  IMAGINARY 

COMPONENTS  OF  THE  SM . 
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2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 


Fact=Rsrn*R5n+l5M*Isn 
Sm(Phi  )  =  5QR(Fact  ) 
Ds(Phi  )=4«(Fact  )/(Ka*Ka  ) 

NEXT  Phi 


CALCULATING  THE  ARRAY  REPRESENTING 
THE  SM( PHI  ). 

CALCULATING  THE  ARRAY  REPRESENTING 
THE  DIFFERENTIAL  SCATTERING  CROSS 
SECTION  NORMALIZED  BY  THE  CROSS 
SECTIONAL  AREA  OF  THE  SPHERE. 


! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

I 
I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
I 
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2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2980 
2990 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
31  10 
3120 
3130 
3140 
3150 


SUB  Drsbdmax  ,Ka  ,Dj(  »  )  ,J(  *  )  ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM  THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  REGULAR  SPHERICAL 

REM  BESSEL  FUNCTIONS  (DRSB)  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX . 

REM  THE  CALCULATION  OF  THE  DERIVATIVE  INVOLVES  THE  RSB'S  OF  THE  SAME 

REM  AND  PREVIOUS  ORDER. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  CALCULATING  THE  DERIVATIVES  OF  THE  RSB  +++++++++++++++++++++++++++ 


I  CALCULATE  THE  RSB  VALUES  J(«). 
! 

!  ANGLE  UNITS  FOR  DRSB  CALC'NS. 

!  CALCULATE  THE  INITIAL  DRSB  DJ(0) 

i  CALCULATE  THE  REMAINING  DRSB. 


CALL  Rsb(Lmax ,Ka ,J( » ) ) 

RAD 

Dj(0)=(-J(0)+COS(Ka ) )/Ka 

FOR  L=1  TO  Lnax 

Dj(L)  =  -((L+1  >/Ka)*J(L)+J(L-1  ) 
NEXT  L 


I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I+  +  +  +  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

I 

! 
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31G0 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 


I $$$$$$$$S$S$$$$$$$S$$$$$$$$S$$$$$$$$$$$$S$S$$$$$$$$$$$5$$$3S$$$$$$$$$S$ 
SUB  Rsb<Lmax  ,Ka  ,J( *  )  ) 


REM 

REM 

REM 

REM 

REM 

! 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
THIS  MODULE  CALCULATES  THE  REGULAR  SPHERICAL  BESSEL  FUNCTIONS 
(RSB)  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .    THE 
CALCULATION  USES  THE  CONTINUED  FRACTION  APPROACH. 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


REAL  Rl 
REAL  Numr 

REAL  Denr 

REAL  Scexp 

REAL  Nflag 
REAL  Dflag 
REAL  Bnm 

REAL  Bine 


THE  EVOLVING  RATIO  OF  J(0)/J(L). 

THE  EVOLVING  NUMERATOR  FACTOR  USED 
I  IN  CALCULATION  OF  RL. 
! 

!  THE  EVOLVING  DENOMINATOR  FACTOR 
!  USED  IN  CALCULATION  OF  RL. 
I 

!  EVOLVING  SCALING  EXPONENT  FOR  VERY 
I  LARGE  OR  VERY  SMALL  VALUES  OF  RL. 
! 

I  0  FLAG  FOR  NUMR. 

0  FLAG  FOR  DENR. 

EVOLVING  TERM  USED  IN  THE 
CALCULATION  OF  NUMR  AND  DENR. 

INCREMENT  USED  IN  THE  EVOLUTION 
OF  BNM. 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  CALCULATING  THE  RSB  ++++++++++++++++++++++++++++++++++++++++++++++ 

! 


RAD 


J(0>  =  (SIN(Ka  )  )/Ka 


FOR  L-l  TO  Lnax 
Rl  =  1  .0 
Nunr*0 . 
Denr"1© . 
Scexp-0 . 
Nf lao-0. 
Df lag-0. 
Bnn=1 .0/Ka 
Binc-2.0/Ka 


!  ANGLE  UNITS  FOR  CALC'NS. 

! 

!  •«•  CALCULATE  THE  INITIAL  »•* 

I  •*•  RSB,  J(0).  ••• 

I 

!  »»•  CALCULATING  THE  RATIO  »»* 

!  «*»  J<0)/J(L).  •  «• 

!  INITIALIZING  RL . 

I  INITIALIZING  NUMR. 

!  INITIALIZING  DENR. 

I  INITIALIZING  SCEXP. 

!  INITIALIZING  NFLAG. 

INITIALIZING  DFLAG. 

INITIALIZING  BNM. 

CALCULATE  BINC. 
I 
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3710 

3720 

FOR  Ordcnt=1  TO  L 

3730 

Bnm=Bnn+B l nc 

3740 

IF  Nflag=0  THEN 

3750 

Nufir  =  Bnn-Numr 

37G0 

IF  Numr<>0  THEN 

3770 

Rl=Rl*Numr 

3780 

Numr= 1 . 0/Numr 

3790 

3800 

IF  ABS(R1  )>10 

3810 

ELSE 

3820 

Nflag=1 

3830 

END  IF 

3840 

ELSE 

3850 

Nflag=0 

3860 

END  IF 

3870 

NEXT  Ordcnt 

3880 

3890 

3900 

3910 

3920 

WHILE  NunrODenr  OR 

3930 

Bnro  =  Bnn  +  Binc 

3940 

3950 

IF  Nflag=0  THEN 

3960 

Numr=Bnm-Numr 

3970 

IF  Nunr<>0  THEN 

3980 

Rl=Rl»Numr 

3990 

Numr= 1 . 0/Nunr 

4000 

4010 

IF  ABS(R1  )>10 

4020 

ELSE 

4030 

Nf lag=1 

4040 

END  IF 

4050 

ELSE 

4060 

Nf lag=0 

4070 

END  IF 

4080 

4090 

IF  Dflag=0  THEN 

4100 

Denr=Bnm-Denr 

41  10 

IF  Denr<>0  THEN 

4120 

Rl=Rl/Oenr 

4130 

Denr= 1 . 0/Denr 

4140 

4150 

IF  ABS(R1  )>10 

4160 

ELSE 

4170 

Dflag=1 

4  180 

END  IF 

4190 

ELSE 

4200 

Df lag=0 

4210 

END  IF 

4220 

4230 

END  UHILE 

4240 

*  CALCULATE  THE  L  "NAKED"     « 

*  NUMR  FACTORS  OF  J(0)/J<L).  » 
INCREMENT  BNM. 
CHECK  NUMR=0  FLAG  NOT  SET; 
UPDATE  NUMR  IF  NUMRO0. 

UPDATE  RL  IF  UPDATED  NUMRO0. 
PREPARE  NUMR  FOR  NEXT  EVOLUTION. 
CHECK  SCALING  OF  UPDATED  RL. 
'250  THEN  CALL  Seal e< R  1  ,Scexp  ) 


!  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 
!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


I  *  CALCULATE  THE  REMAINING  NUMR  • 

!  *  AND  DENR  TERMS  OF  J(0)/J(L).  « 

!  NOTE  WHEN  NUMR=DENR  THEN 

!  RL  =  J(0  )/J(L  ). 
Nflag=1  OR  Dflag=1 

I  INCREMENT  BNM. 
! 

!  CHECK  NUMR=0  FLAG  NOT  SET; 

!  UPDATE  NUMR  IF  NUMRO0. 

!  UPDATE  RL  IF  UPDATED  NUMRO0. 

I  PREPARE  NUMR  FOR  NEXT  EVOLUTION, 

!  CHECK  SCALING  OF  UPDATED  RL. 
'250  THEN  CALL  Sea le( Rl  .Scexp  ) 


!  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 

!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 

I 

!  CHECK  DENR-0  FLAG  NOT  SET; 

!  UPDATE  DENR  IF  DENRO0. 

I  UPDATE  RL  IF  UPDATED  DENRO0. 

!  PREPARE  DENR  FOR  NEXT  EVOLUTION. 

I  CHECK  SCALING  OF  UPDATED  RL. 
'250  THEN  CALL  Seal e< Rl  .Scexp  ) 


I  IF  DFLAGO0  THE  PREVIOUS  STEPS  ARE 
!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


i 
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!     »♦•    J(L  )    DETERMINED.    »*« 

I 


4250 

4260  J(L  )  =  J(0  )/(Rl* 10"(Scexp  )  ) 

4270 

4280      NEXT    L 

4290  ! 

4300  ! 

43  10       i +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +++  +  +  +  +  +  +  +  +  +  +  + 

4320 

4330 

4340 

4350      SUBEND 

43E0 


++++    MODULE    CONCLUSION    ++++++++++++++++++++++++++++++++++++++++++++++++ 


4370 
4380 
4390 
4400 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

■T  41  Tp  T>  y  4?  j  j  %£  sT  ■P  -P  X  J  J1  T-  -P  ^'  ■— ■  T'  3^  ■T  J^  'P  T1  T1  T  ■P  -P  P  -P  -P  -P  ■P  •*'  -P  T'  T ;  -P  J  JJ  jyTjJJijjyjj^I^rjjyjiijX'rJ'rt'r'rJJ 
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4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
45G0 
4570 
4580 
4580 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
4700 
4710 


1$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB  Scale( Rl  ,Scexp  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM    THIS  MODULE  SCALES  RL  TO  MAINTAIN  NUMERICAL  ACCURACY  IN  THE 
REM    ONGOING  CALCULATION  OF  RL. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!+  +  +  +  SCALING  RL  +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

l 
■  ! 
IF  Rl>10"250  THEN 

R1=R1 *  1 0"( -250 ) 

Scexp=Scexp+250 
END  IF 
IF  RK10-M-2S0)  THEN 

R1=R1*10A250 

Scexp=Scexp-250 
END  IF 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

I 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 


!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$£$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB  Disb(Lna>.  ,Ka  ,Dy(  »  >  ,Y(  «  )  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM    THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  IRREGULAR  SPHERICAL 
REM    BESSEL  FUNCTIONS  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .   THE 
REM    CALCULATION  OF  THE  DERIVATIVE  USES  THE  ISB'S  OF  THE  SAME  AND 
REM    PREVIOUS  ORDER. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I++++  CALCULATING  THE  DERIVATIVES  OF  THE  ISB  +++++++++++++++++++++++++++ 

i 


CALL  Isb(Lmax  ,Ka  ,Y( *  )  ) 

RAD 

Dy(0)  =  (-Y(0)  +  SIN(Ka))/Ka 

FOR    L=1    TO    Lciax 

Dy(L)  =  -((L+1  )/Ka  )*Y(L)  +  Y(L-1  ) 
NEXT    L 


I  CALCULATE  THE  ISB  VALUES. 
I 

!  ANGLE  UNITS  FOR  CALC'NS. 

!  THE  INITIAL  ISB  DERIVATIVE. 

!  THE  REMAINING  ISB  DERIVATIVES 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

! 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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5040 
5050 
5060 
5070 
5080 
5090 
5100 
51  10 
5120 
5130 
5140 
5150 
51G0 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 
52G0 
5270 
5280 
5290 
5300 
5310 
5320 
5330 


I  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB    Isb(Lna.x  ,Ka  ,Y(  *  )  ) 

REM         ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM         THIS    MODULE    CALCULATES    THE    IRREGULAR    SPHERICAL    BESSEL    FUNCTIONS 

REM        OF    ARGUMENT    KA   AND      ORDER    L    THROUGH   LMAX .         THE    CALCULATION    IS 

REM         ACCOMPLISHED    BY    FORWARD    RECURSION. 

REM         ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!+  +  +  +    CALCULATING    THE    ISB    ++++++++++++++++++++++++++++++++++++++++++++++ 


RAD 

Y(0)=-COS(Ka )/Ka 

Y(  1  )  =  (  Y(0  >-SIN(Ka  )  )/Ka 

FOR    L=2    TO   Lmax 

Y(L)  =  ((2*L-1  >/Ka)*Y<L-1  )-Y(L-2) 
NEXT    L 


ANGLE    UNITS    FOR    CALC'NS. 
INITIAL    ISB'S    CALCULATED, 

!    REMAINING    ISB'S   CALC'D. 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++    MODULE    CONCLUSION    ++++++++++++++++++++++++++++++++++++++++++++++++ 

l 
SUBEND 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

;p;pj)lj)lpjpij}lplplpip5)lj)lp  Ip  -P  ip  -P  -I'  i1  i  -P  it  >P  «fe  »fe  *fc  *P  lp  *P  «P  *P  ^^lp^ji^!j^^^^!p^lp^^tp^^^^^i^^pip!pip  3)  3)  Jfc  JB  3)  ffi  X  B  8  5fc  *S 
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5340 
5350 
53G0 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570 
5580 
5590 
5600 
5610 
5620 
5630 
5640 
5650 


SUB  Leg(l_max  ,Ka  ,Phi  ,P(  *  )  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  CALCULATES  THE  LESENDRE  POLYNOMIALS  OF  ARGUMENT 

REM    COS(PHI)  AND  ORDERS  0  THROUGH  LMAX .   CALCULATION  IS  ACCOMPLISHED 

REM    BY  FORWARD  RECURSION. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  CALCULATING  THE  LEGENDRE  POLYNOMIALS  +++++++++++++++++++++++++++++ 


DEG 
X=COS(Phi ) 


ANGLE  UNITS  FOR  CALC'NS. 
ARGUMENT  OF  THE  LEG.  POL'S, 

INITIAL  LEG.  POL'S.  CALC'D, 


P(0)=l 
P<  1  )  =  X 
FOR  L=2  TO  Lmax  !  REMAINING  LEG.  POL'S.  CALC'D. 

P(L)=<2-<  1/L))*X*P(L-1  )-(  1-(  1/L)  )*P(L-2  ) 
NEXT  L 

! 

I 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

! 
! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


I 
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5G60 
5670 

5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5780 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

5940 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 

6100 

6110 

6120 

6130 

6140 

615C 

6160 

6170 

6180 

6190 


SUB  OutpUC  ,Freq  .A  ,Ka  .Sm<  *  )  ,Ds<  •  >  ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM  THIS  MODULE  CONTROLS  THE  OUTPUT  OF  THE  CALCULATED  DSCS  AND  SM  IN 

REM  GRAPHICAL  PRESENTATIONS. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Desire 
REAL  Rpt 


INTERACTIVE  LOOP  CONTROL. 

POLARPLT  SCALING  CONTROL)  A 
REPETITION  INDEX. 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  DETERMINING  THE  DESIRED  OUTPUT  +++++++++++++++++++++++++++++++++++ 

I 
Des  ire°=0 
WHILE  DesireOI  AND  Desire<>2  AND  Desire<>3 

PRINT  "DO  YOU  DESIRE" 

PRINT  "      (1)  A  PLOT  OF  THE  DSCS," 

PRINT  "      (2)  A  PLOT  OF  THE  SM , " 

PRINT  "   OR  (3)  BOTH?" 

INPUT  Desire 
END  WHILE 
PRINTER  IS  701 


WHILE  Desire<4 

IF  Desire=1  OR  Desire=3  THEN 
FOR  Rpt-1  TO  2 

CALL  Polarplt<Ds( *  )  ,Rpt  ) 
PRINT 
NEXT  Rpt 
END  IF 


i  PLOTTING  THE  DSCS 


IF  Desire=2  THEN 
CALL  Recplt(Sm( *  )  ) 
PRINT 

END  IF 

PRINT  "  " ,"C  =" ;C;"M/S. " 

PRINT  "  " ,"F  =" ;Freq;"HZ 
PRINT  "  "  ,"A  ="iAi"M." 
PRINT  "  "  ,  "K«A  =" ikai " . " 
PRINT  "  "  ."NORMALIZED  "; 
PRINT  "BCKSCTR  =  "  fDs<  1  8i 
PRINT  CHR$( 12 ) 


!  PLOTTING  THE  SM. 


ECHO  SIGNIFICANT  CALCULATED  VALUES 

AND  INPUT  PARAMETERS. 


PRINTER  FORM  FEED, 


IF  Desire=1  OR  Desire=2  THEN  Desire=4 
IF  Desire=3  THEN  Desire=2 
END  WHILE 

I 
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G200 

6210 

G220 

6230 

6240 

6250   SUBEND 

6260 

6270 

6280 

6290 

6300 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

! 
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6310 
6320 
6330 
6340 
6350 
6360 
6370 
6380 
6390 
6400 
6410 
6420 
6430 
6440 
6450 
6460 
6470 
6480 
6490 
6500 
6510 
6520 
6530 
6540 
6550 
6560 
6570 
6580 
6590 
6600 
6610 
6620 
6630 
6640 
6650 
6660 
6670 
6680 
6690 
6700 
6710 
6720 
6730 
6740 


I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB  Polarpl t (Ds( *  )  ,Rpt  ) 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
THIS  MODULE  MAPS  THE  CALCULATED,  NORMALIZED  DSCS  ARRAY  DS( *  )  ON 
LINEAR  POLAR  PLOTS.  THE  FIRST  PLOT  INCLUDES  THE  ENTIRE  DSCS, 
WHILE  THE  SECOND  PLOT  SHOWS  THE  SCATTERING  DETAIL  AROUND  THE 
NORMALIZED  RANGE.    PLOT  SELECTION  IS  GOVERNED  BY  A  REPETITION 
INDEX,  Rpt,  AND  IS  CONTROLLED  IN  THE  CALLING  MODULE  (OUTPT). 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


REAL  Xgunax 

REAL  Ygumax 

REAL  6umax 

REAL  Dtmx 
REAL  Endpt 

REAL  Ring 
REAL  Tick 


THE  MAXIMUM  ABSCISSA  VALUE  IN 
GRAPHICS  DISPLAY  UNITS. 

THE  MAXIMUM  ORDINATE  VALUE  IN 
GRAPHICS  DISPLAY  UNITS. 

THE  MINIMUM  OF  XGUMAX  AND  YGUMAX i 
USED  TO  DETERMINE  PLOTTING  AREA. 

THE  MAX.  VALUE  OF  THE  ARRAY  DS(») 

THE  MAX.  RADIAL  SCALE  VALUE  IN 
USER  DEFINED  UNITS. 

GRID  RING  SPACING  IN  USER  UNITS. 

RADIAL  AXIS  TICK  MARK  SPACINGS  IN 
USER  DEFINED  UNITS. 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  INITIALIZING  THE  PLOTTER  +++++++++++++++++++++++++++++++++++++++++ 

! 
GINIT 

PLOTTER  IS  CRT,  "INTERNAL- 
GRAPHICS  ON 
GCLEAR 

I 
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6750 
6760 

6770 

6780 

6790 

6800 

6810 

6820 

6830 

6840 

6850 

6860 

6870 

6880 

6890 

6900 

6910 

6920 

6930 

6940 

6950 

6960 

6970 

6980 

6990 

7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

7080 

7090 

7100 

71  10 

7120 

7130 

71  40 

7150 

7160 

7170 

7180 

7190 

7200 

7210 

7220 

7230 

7240 

7250 

7260 

7270 

7280 

7290 


I++++  GRAPH  CONSTRUCTION  AND  LABELLING  +++++++++++++++++++++++++++++++++ 

I 

! 
Xgunax=100*MAX<  1  .RATIO  ) 
Ygunax=100»MAX(  1  ,1/RATIO) 
Gunax=MIN(  Xgunax  .Yguriax  > 


CSIZE  3 

LINE  TYPE  1. 

MOVE  6unax ,6umax/2 

LORG  8 

LABEL  "0" 

MOVE  Gumax/2  .Gumax 

LORG  6 

LABEL  "PI/2" 

MOVE  0, Gunax/2 

LORG  2 

LABEL  "PI" 

MOVE  Gunax/2  ,0 

LORG  4 

LABEL  "3*PI/2" 

VIEWPORT  5 ,6unax-5  ,5  ,Gunax-5 


Dtmax=0 

IF  Rpt=l  THEN 

FOR  1=0  TO  360 

IF  Dtmx<Ds(I)  THEN  Dtmx-Dsf I > 

NEXT  I 
END  IF 

Endpt=INT(Dtnx+1 ) 
IF  Endpt<2  THEN  Endpt=2 
SHOW  -Endpt  .Endpt  ,-Endpt  ,Endpt 

LINE  TYPE  4 

Ring=INT(Endpt/5+l  ) 

Tick=.2»Ring 

AXES  Tick .Tick .0 .0,5.5,2 


I 

!  *«  LABELLING  MAJOR  RADIALS  ** 


»»  DEFINES  THE  TOTAL  PLOTTING  •* 
»«  AREA  IN  GRAPH  UNITS.        ** 

I  **  SCALING  THE  GRAPH  TO  ** 

!  **  USER  UNITS.  ** 

I  DETERMINING  MAX.  RADIUS  VALUE 

!  REQ'D  FOR  THE  PLOT. 


LINE  TYPE  3 

DEG 

LORG  7 

FOR  R=1  TO  Endpt/Ring 

LINE  TYPE  1 

MOVE  0  ,R»Ring 

LABEL  R*Ring 

LINE  TYPE  3 

FOR  Angle=0  TO  360 

PLOT  R-Ring«COS( Angle  )  ,R*Ring»SIN( Angle  ) 

NEXT  Angle 

PENUP 
NEXT  R 


!  ISOTROPICALLY  SCALES  THE  GRAPH, 


**  MAJOR  RADIALS  CONSTRUCTION  *• 


! 

I  **  RANGE  RING  CONSTRUCTION  «* 


!  LABELLING  THE  RANGE, 
!  DRAWING  THE  RING. 
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7300  ! 

7310   FOR  Angle=30  TO  150  STEP  30         I  •«  30  DEGREE  RADIALS  *« 

7320     IF  Angle<>90  THEN  !  **  CONSTRUCTION. 

7330      PLOT  Endpt«COS< Angle  )  ,Endpt*SIN( Angle  )  ,-2 

7340       DRAW  Endpt*COS(  180  +  Angle >  ,Endpt«SIN(  180  +  Angle > 

7350    END  IF 

7360   NEXT  Angle 

7370   PENUP 

7380  ! 

7390  ! 

7400   I ++++++++++++++++++++++++++++++++++++++++++++++++++ +++++++++++++++++++++ 

7410   I 

7420   !+  +++  PLOTTING  THE  CALCULATED  DSCS  +++++++++++++++++++++++++++++++++++++ 

7430  ! 

7440   LINE  TYPE  1 

7450   DEG 

7460   FOR  1=0  TO  360  !  PLOTTING  RECTANGULAR  COORDINATES 

7470    PLOT  Ds< I  )*COS(  I  )  ,Ds( I  >*SIN( I  )     !  IN  POLAR  FORMAT. 

7480   NEXT  I 

7490  ! 

7500   !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7510   ! 

7520   !++++  PRINTIN6  THE  PLOT  ++++++++++++++++++++++++++++++++++++++++++++++++ 

7530  I 

7540   DUMP  GRAPHICS  #701 

7550  ! 

7560   !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7570   ! 

7580   I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

7590  ! 

7600   6RAPHICS  OFF 

7610   SUBEND 

7620  ! 

7630   !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7640   ! 

7G50   1$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

7660   ! 
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7G70 
7680 
7690 
7700 
7710 
7720 
7730 
7740 
7750 
7760 
7770 
7780 
7790 
7800 
7810 
7820 
7830 
7840 
7850 
7860 
7870 
7880 
7890 
7900 
7910 
7920 
7930 
7940 
7950 
7960 
7970 
7980 
7990 
8000 
8010 
8020 
8030 
8040 
8050 
8060 
8070 
8080 
8090 
8100 
8110 
8120 
8130 
8140 
8150 
8160 


!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Recplt ( Sm( •  )  > 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM  THIS  MODULE  MAPS  THE  CALCULATED  SCATTERING  MODULUS  ARRAY  SM( * )  ON 

REM  A  SEMI-LOG  PLOT. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Xgumax 

REAL  Ygunax 

REAL  Dtmx 
REAL  Dtmn 
REAL  Mxy 

REAL  Mny 


THE  MAXIMUM  ABSCISSA  VALUE  IN 
GRAPHICS  DISPLAY  UNITS. 

THE  MAXIMUM  ORDINATE  VALUE  IN 
GRAPHICS  DISPLAY  UNITS. 

THE  MAX.  VALUE  OF  THE  ARRAY  SM<») 

THE  MIN.  VALUE  OF  THE  ARRAY  SM( »  ). 

THE  MAX.  ORDINATE  SCALE  VALUE  IN 
USER  DEFINED  UNITS. 

THE  MIN.  ORDINATE  SCALE  VALUE  IN 
USER  DEFINED  UNITS. 


i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  INITIALIZING  THE  PLOTTER  +++++++++++++++++++++++++++++++++++++++++ 

I 
GINIT 

PLOTTER  IS  CRT,  "INTERNAL" 
GRAPHICS  ON 
GCLEAR 

i 
I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
i 
!+  +  +  +  GRAPH  CONSTRUCTION  AND  LABELLING  +++++++++++++++++++++++++++++++++ 


I 


Xgumax=100«MAX<  1  .RATIO  ) 
Ygunax=100»MAX(  1  ,1 /RATIO  ) 

CSIZE  4 

LORG  4 

MOVE  Xgunax/2  .0 

LABEL  "PHI  (DEGREES)" 


i  *«  LABELLING  THE  ABSCISSA  AXIS  «« 
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8170 
8180 
8190 
8200 
8210 
8220 
8230 
8240 
8250 
82G0 
8270 
8280 
8290 
8300 
8310 
8320 
8330 
8340 
8350 
8360 
8370 
8380 
8390 
8400 
8410 
8420 
8430 
8440 
8450 
8460 
8470 
8480 
8490 
8500 
8510 
8520 
8530 
8540 
8550 
8560 
8570 
8580 
8590 
8600 
8610 
8620 
8630 
8640 
8650 
8660 
8670 
8680 
8690 
8700 
8710 


DE6 

LDIR  90 

LORG  6 

MOVE  0  ,Ygumax/2 

LABEL  "5CTRG  MODULUS" 

LDIR  0 

VIEWPORT  15  ,Xgunax-5  ,  10  .Ygunax-5 


Dtnx=0 

Dtnn=Sm(  180  ) 

FOR  1=0  TO  180 

IF  Dtmx<Sm(I)  THEN  Dtmx=Sm< I  ) 
IF  Dtmn>Sm< I  )  THEN  Dtnn  =  Sm< I  ) 

NEXT  I 

Mxy=INT(LGT(Dtmx  )  )+1 

Mny=INT(L6T(Dtr<in  )  ) 

WINDOW  0  ,  180  ,Mny  ,Mxy 

AXES  30,1  ,0,Mny  ,1,1,1.5 

LINE  TYPE  3 
FOR  I=Mny  TO  Mxy 
FOR  J=l  TO  9 

MOVE  0,I+LGT( J  ) 
PLOT  180,I  +  LGT( J  )  ,-1 
PENUP 
NEXT  J 
NEXT  I 

FOR  1=60  TO  180  STEP  60 

MOVE  I  ,Mny 

PLOT  I  ,Mxy  ,-1 

PENUP 
NEXT  I 

CLIP  OFF 
LINE  TYPE  1 

L0R6  6 

FOR  1=0  TO  180  STEP  30 

MOVE  I  ,Mny 

LABEL  I 
NEXT  I 

L0R6  8 

FOR  I=Mny  TO  Mxy 

MOVE  0,1 

LABEL  10"I 

MOVE  0,I+LGT(5) 

LABEL  ( 10"< 1+1  )  )/2 
NEXT  I 


I 

!  **    LABELLING  THE  ORDINATE  AXIS  »« 


•  •  GRAPH  CONSTRUCTION  »» 
DEFINES  THE  TOTAL  PLOTTING  AREA 
IN  GRAPH  UNITS. 

**  ORDINATE  USER  UNITS  *• 
*»  DETERMINED.         «* 


I 

!  ANISOTROPICALLY  SCALES  THE  GRAPH 

!  TO  USER  UNITS. 

!  PLOTS  THE  AXES. 


PLOTS  THE  HORIZONTAL 
LOGARITHMIC  GRID  LINES  I.A.W, 
THE  VERTICAL  SCALE. 


PLOTS  THE  VERTICAL  GRID  LINES 
I.A.W.  THE  HORIZONTAL  SCALE. 


**  LABELLING  THE  ABSCISSA  SCALE  •* 


I  **    LABELLING  THE  ORDINATE  SCALE  ♦• 

!  POWER  OF  10  LABELS. 

!  HALF  THE  NEXT  POWER  OF  10  LABELS. 


93 


8720 
8730 
8740 
8750 
8760 
8770 
8780 
8790 
8800 
8810 
8820 
8830 
8840 
8850 
8860 
8870 
8880 
8890 
8900 
8910 
8920 
8930 
8940 
8950 
8960 
8970 
8980 
8990 
9000 
9010 


CLIP  ON 


I+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  PLOTTING  THE  CALCULATED  SCATTERING  MODULUS  +++++++++++++++++++++++ 


LINE  TYPE  I 
FOR  1=0  TO  180 

PLOT  I  ,LGT(Sn( 1  +  180) ) 
NEXT  I 


THIS  PLOT  IS  FOR  THE  PLANE  WAVE 
INCIDENT  FROM  THE  RIGHT  (0  DEG'S). 
USE  SM( I  )  IN  PLACE  OF  SM( 1  +  180) 
FOR  A  PLANE  UAUE  INCIDENT  FROM  THE 
LEFT  ( 180  DEG'S). 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  PRINTING  THE  PLOT  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 
DUMP  GRAPHICS  #701 

! 
!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 
I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 
GRAPHICS  OFF 
SUBEND 

I 
i +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +++  +  +  +++  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +++  +  +  +  +  +  +  +  +  +  +  +++  +  + 
! 
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APPENDIX    C 


BACK5CATTERED    CROSS  SECTION   PROGRAM 


Examples  of  the  output  of  this  program  are  Figures  6  and  7,  located  in 


Chapter  III. 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

1  10 

120 

130 

140 

150 

1G0 

170 

180 

190 

200 

210 

220 

230 

240 

250 

2G0 

270 

280 

290 

300 

310 

320 

330 

340 

350 

3G0 

370 


I  tt .  tJ)  T.  8! i  T.  fti  tt .*T, ,  *?* it  *V.  $\  *T.  it i  *£  T^  ^  *Tl  ^  ^Ti  *fc  *t  ^  ^  *t  ^T*.  *T^  *fc  T.  'T.  ^  ^.  ^.  T\  ^.  ^.  ^.  ^,  ^  ^  ^  H^  ^  ^  ^,  ^,  ^  ^.  T,  ^  1 1 1 1 1  ^  t  ^  1 1  *tl  ^  ^1  *Tl  *Tt  t  *t  t  f£  *R  *t 

i 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

ACOUSTICALLY  HARD  SPHERE  BACKSCATTERING  OF  PLANE  ACOUSTIC  WAVES 


REM 

REM 

REM 

REM     THIS  PROGRAM  CALCULATES  AND  GRAPHS  THE  NORMALIZED,  BACKSCATTERED 

REM    CROSS  SECTION  FOR  ACOUSTICALLY  HARD  SPHERES  OF  A  SELECTED  RANGE 

REM    OF  DIAMETERS  D  ENSONIFIED  BY  AN  INCIDENT  PLANE  WAVE . 

REM 

REM 

REM 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  MAIN  PROGRAM  VARIABLE  DECLARATION  AND  DEFINITION  +++++++++++++++++ 


REAL  Nfas(301 ) 

REAL  K 
REAL  A 
REAL  Ka 

REAL  Kanin 
REAL  Kama* 
REAL  Stp 


THE  ARRAY  REPRESENTING  THE 
NORMALIZED,  BACKSCATTERED  CROSS 
SECTION  (NBCS). 

THE  WAVE  NUMBER,  2«PI*FREQ/C. 

RADIUS  OF  THE  SPHERE. 

K*A,  THE  ARGUMENT  OF  THE  REGULAR 
AND  IRREGULAR  SPHERICAL  BESSEL 
FUNCTIONS  AND  THEIR  DERIVATIVES. 

INTEGER  OF  THE  MIN.  VALUE  OF  KA. 

INTEGER  OF  THE  MAX.  VALUE  OF  KA. 

KA  INCREMENT. 


t  +  +  +  +  +  +  +  +  +++  +  +  +  +  +  +  +  +  +++  +  +  +  +  +  +  +++  +  +++  +++  +  +  +++++  +++  +  +++++  +  +  +  +  +  +  +++  +  +  +  +  +  +  +  + 
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380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 


++++  MAIN  PROGRAM  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 


CALL  Bcksctr(Nbs(*),Kanin,Kanax,Stp) 
CALL  Recplt(Nbs( *  )  ,  Kan  in  .Kama* ,Stp ) 


CALCULATE  THE  NBS( * ). 
OUTPUT  THE  NBS(*  ). 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MAIN  PROGRAM  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++ 

i 

PRINTER  IS  701 
PRINT  CHR$(  12  ) 
END 

I 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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560 

570 

580 

590 

G00 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 


I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB  Bck5ctr(Nbs(*),Kamin,Kanax,Stp) 


REM 
REM 
REM 
REM 
REM 
I 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


THIS  MODULE  CALCULATES  THE  ARRAY  NBS( *  )  REPRESENTING  THE 
BACKSCATTERED  CROSS  SECTION  NORMALIZED  BY  THE  CROSS  SECTIONAL  AREA 
OF  THE  SPHERE. 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


REAL  L 


REAL  Lmax 


REAL  J(S1 ) 


REAL  Dj  <  51  ) 


REAL  Y(51  ) 


REAL  Dy(51  ) 


REAL  P( 5 1  > 


REAL  Rsm 


REAL  Ism 


REAL  Fac 


ORDER  OF  THE  FUNCTION  CALCULATED 
IN  A  SUBROUTINE.   IT  IS  USED 
AS  A  LOOP  INDEX  AND  TO  DESIGNATE 
THE  ARRAY  ELEMENT  FOR  THE  FUNCTION 
OF  CORRESPONDING  ORDER. 

MAX.  SIGNIFICANT  ORDER  OF  THE  SUM 
DETERMINING  THE  DIFFERENTIAL 
SCATTERING  CROSS  SECTION  (DSCS) 
AND  SCATTERING  MODULUS  ( SM  ) . 

ARRAY  FOR  THE  REGULAR  SPHERICAL 
BESSEL  FUNCTIONS  (RSB)  OF 
ARGUMENT  KA  AND  ORDER  0  THROUGH 
LMAX. 

ARRAY  FOR  THE  DERIVATIVES  OF  THE 
RSB  (THE  DRSB). 

ARRAY  FOR  THE  IRREGULAR  SPHERICAL 
BESSEL  FUNCTIONS  (  ISB  )  OF 
AR6UMENT  KA  AND  ORDER  0  THROUGH 

LMAX. 

ARRAY  FOR  THE  DERIVATIVES  OF  THE 
ISB  (THE  DISB). 

ARRAY  FOR  THE  LE6ENDRE 
POLYNOMIALS  OF  ARGUMENT  -1 . 

THE  REAL  COMPONENT  OF  THE 
SCATTERING  MODULUS. 

THE  IMAGINARY  COMPONENT  OF  THE 
SCATTERING  MODULUS. 

FACTOR  COMMON  TO  RSM  AND  ISM. 

++++  NOTE  ++++ 

THE  ARRAY  DIMENSIONS  (30  +  21  )  NEED 

TO  BE  INCREASED  TO  ACCOMODATE 

KAMAX>30. 
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1 100 
1 1 10 

1  120 
1  130 
1  140 
1  150 
1  1G0 
1  170 
1  180 
1  190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
13E0 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1450 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 


i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  CALCULATING  THE  NBCS  ARRAY  +++++++++++++++++++++++++++++++++++++++ 


PRINTER  IS  1 
INPUT  Kamin.Kamax 

Kanin=INT( Kami n ) 
Karoax  =  INT(  Kamax  ) 
Stp=. 1 

FOR  Ka=Kanin+Stp  TO  Kanax  STEP  Stp 
Lmax=INT(Ka+21  ) 
CALL  Drsb<Lnax  ,Ka ,Dj( • ) ,J( ♦ ) ) 

CALL  Disb(Lmax  ,Ka  ,Dy( •  )  ,Y<* ) ) 

CALL  Leg180(Lmax  ,Ka  ,P< * ) ) 


•  DETERMINE  DESIRED  RANGE  OF  KA 
!  UALUES. 


CALCULATE  THE  DRSB  OF  ARGUMENT  KA 
AND  ORDERS  0  THROUGH  LMAX. 
CALCULATE  THE  DISB  OF  ARGUMENT  KA 
AND  ORDERS  0  THROUGH  LMAX. 
CALCULATE  P<»)  OF  ARGUMENT  -1  AND 
ORDERS  0  THROUGH  LMAX. 
CALCULATING  THE  REAL  AND  IMAGINARY 
COMPONENTS  OF  THE  SCATTERING 
MODULUS. 


Rspi=0 
lsm<=0 
FOR  L=0  TO  Lroax 

Fac  =  (Dj(L  )*(2«L+1  )*P(L  )  )/(Dj(L)*Dj(L)  +  Dy(L)*Dy(L)  ) 
Rsn=Rsn+(Fac*Dy(L ) ) 
Ism=Ism+(Fac*Dj< L ) ) 
NEXT  L 

!  CALCULATING  THE  ARRAY  REPRESENTING 
!  THE  NBCS. 
Nbs( (Ka-Kamin  )/Stp)  =  4*(R5n*R5n+Ism*Isn)/(Ka*Ka> 

! 

NEXT  Ka 

I 

I 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

! 
i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


i 
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1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 


i $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Drsb(Lma>.  ,Ka  ,Dj(  «  )  ,J(  *  >  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  RE6ULAR  SPHERICAL 

REM    BESSEL  FUNCTIONS  (DRSB)  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX 

REM    THE  CALCULATION  OF  THE  DERIVATIVE  INVOLVES  THE  RSB'S  OF  THE  SAME 

REM    AND  PREVIOUS  ORDER. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  CALCULATING  THE  DERIVATIVES  OF  THE  RSB  +++++++++++++++++++++++++++ 

! 

I 


CALL  Rsb(Lmax  ,Ka  ,J( *  )  > 

RAD 

Dj(0)  =  (-J(0)  +  COS(Ka  )  >/Ka 

FOR  L-1  TO  Lmax 

Dj(L  )  =  -(  (L+1  )/Ka)»J(L)  +  J(L-1  > 
NEXT  L 


CALCULATE  THE  RSB  VALUES  J(*). 

ANGLE  UNITS  FOR  DRSB  CALC'NS. 
CALCULATE  THE  INITIAL  DRSB  DJ(0). 
CALCULATE  THE  REMAINING  DRSB. 


! 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 
SUBEND 

! 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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1870 
1880 
1890 
1  900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
21  10 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 


SUB  Rsb(Lrnax  ,Ka,J(«  )  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM    THIS  MODULE  CALCULATES  THE  REGULAR  SPHERICAL  BESSEL  FUNCTIONS 
REM    (RSB)  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX  .    THE 
REM    CALCULATION  USES  THE  CONTINUED  FRACTION  APPROACH. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
i 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 

l 


REAL  Rl 
REAL  Numr 

REAL  Denr 

REAL  Scexp 

REAL  Nflag 
REAL  Dflag 
REAL  Bnm 

REAL  Bine 


THE  EVOLVING  RATIO  OF  J(0)/J(L). 

THE  EVOLVING  NUMERATOR  FACTOR  USED 
IN  CALCULATION  OF  RL . 

THE  EVOLVING  DENOMINATOR  FACTOR 
USED  IN  CALCULATION  OF  RL. 

EVOLVING  SCALING  EXPONENT  FOR  VERY 
LARGE  OR  VERY  SMALL  VALUES  OF  RL. 

0  FLAG  FOR  NUMR. 

0  FLAG  FOR  DENR. 

EVOLVING  TERM  USED  IN  THE 
CALCULATION  OF  NUMR  AND  DENR. 

INCREMENT  USED  IN  THE  EVOLUTION 
OF  BNM. 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
++++  CALCULATING  THE  RSB  ++++++++++++++++++++++++++++++++++++++++++++++ 


RAD 


J(0  )  =  (SIN(Ka  )  )/Ka 


FOR  L=1  TO  Lnax 
Rl  =  1  .0 
Numr=0 . 
Denr=0 . 
Scexp=0 . 
Nf lag=0. 
Dflag=0. 
Bnm=1 .0/Ka 
Binc=2.0/Ka 


I 
! 

i  ANGLE  UNITS  FOR  CALC'NS. 

I 

***    CALCULATE  THE  INITIAL  **« 
***  RSB  ,  J(0  ).  ♦«* 

***  CALCULATING  THE  RATIO  »»* 
*»*  J(0)/J(L).  **» 

INITIALIZING  RL. 
INITIALIZING  NUMR. 
INITIALIZING  DENR. 
INITIALIZING  SCEXP. 
INITIALIZING  NFLAG. 
INITIALIZING  DFLAG. 
INITIALIZING  BNM. 
CALCULATE  BINC. 


00 


2420 

2430 

FOR  Ordcnt-1  TO  L 

2440 

BnmcBnm+Binc 

2450 

IF  Nflag=0  THEN 

24G0 

Nunr»Bnm-NupiF 

2470 

IF  Nunr<>0  THEN 

2480 

Rl=Rl*Numr 

2490 

Numr  = 1 . 0/Numr 

2500 

2510 

IF  ABS<R1  )>10 

2520 

ELSE 

2530 

Nflag-1 

2540 

END  IF 

2550 

ELSE 

2560 

Nflag=0 

2570 

END  IF 

2580 

NEXT  Ordcnt 

2590 

2600 

2610 

2620 

2630 

WHILE  NunrODenr  OR 

2640 

Bnm=Bnpi+Binc 

2650 

2660 

IF  Nflag=0  THEN 

2670 

Nunr=Bnm-Numr 

2680 

IF  Nunr<>0  THEN 

2690 

Rl«=Rl*Numr 

2700 

Numr=1 . 0/Numr 

2710 

2720 

IF  ABS<R1  )>10 

2730 

ELSE 

2740 

Nflag-1 

2750 

END  IF 

2760 

ELSE 

2770 

Nflag=0 

2780 

END  IF 

2790 

2800 

IF  Dflag=0  THEN 

2810 

Denr~Bnm-'Denr 

2820 

IF  Denr<>0  THEN 

2830 

R1=R1/Denr 

2840 

Denr=1 . 0/Denr 

2850 

2860 

IF  ABS(R1  )>10 

2870 

ELSE 

2880 

Dflag=1 

2890 

END  IF 

2900 

ELSE 

2910 

Df lag=0 

2920 

END  IF 

2930 

2940 

END  WHILE 

•  CALCULATE  THE  L  "NAKED" 

•  NUMR  FACTORS  OF  J(0)/J(L).  • 
INCREMENT  BNM. 
CHECK  NUMR=0  FLAG  NOT  SET; 
UPDATE  NUMR  IF  NUMRO0. 

UPDATE  RL  IF  UPDATED  NUMRO0. 
PREPARE  NUMR  FOR  NEXT  EVOLUTION, 
CHECK  SCALING  OF  UPDATED  RL. 
250  THEN  CALL  Scale( Rl  .Scexp  ) 


!  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 
I  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


!  *  CALCULATE  THE  REMAINING  NUMR  • 

I  *  AND  DENR  TERMS  OF  J(0)/J(L).  « 

!  NOTE  WHEN  NUMR=DENR  THEN 

!  RL  =  J(0  )/J(L). 
Nflag=1  OR  Dflag=1 

!  INCREMENT  BNM. 

I 

!  CHECK  NUMR=0  FLAG  NOT  SET; 

!  UPDATE  NUMR  IF  NUMRO0. 

!  UPDATE  RL  IF  UPDATED  NUMRO0. 

I  PREPARE  NUMR  FOR  NEXT  EVOLUTION, 

I  CHECK  SCALING  OF  UPDATED  RL. 
250  THEN  CALL  Scale( Rl  , Scexp  ) 


I  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 

I  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 

I 

!  CHECK  DENR=0  FLAG  NOT  SET; 

I  UPDATE  DENR  IF  DENRO0. 

!  UPDATE  RL  IF  UPDATED  DENRO0. 

!  PREPARE  DENR  FOR  NEXT  EVOLUTION. 

!  CHECK  SCALIN6  OF  UPDATED  RL . 
250  THEN  CALL  Scale< Rl  , Scexp  ) 


IF  DFLAGO0  THE  PREVIOUS  STEPS  ARE 
DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


0 


♦•*  J(L)  DETERMINED.  *** 


2950  i 

29S0  ! 

2970     J(L W(0  )/(Rl«l0A(Scexp  )  )         I 

2980  ! 

2990  NEXT  L 

3000  ! 

3010  ! 

3020  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3030  ! 

3040  !++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

3050  ! 

30G0  SUBEND 

3070  ! 

3080  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3090  ! 

31 00  ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

31  10  ! 
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3120   ! 

31 30   ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

3140   SUB  Scale(Rl  .Scexp  ) 

31 50   REM   ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3160  REM   THIS  MODULE  SCALES  RL  TO  MAINTAIN  NUMERICAL  ACCURACY  IN  THE 

3170   REM   ONGOING  CALCULATION  OF  RL. 

31 80   REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3190   ! 

3200   I++++  SCALING  RL  +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3210  ! 

3220  ! 

3230   IF  Rl>10"250  THEN 

3240     Rl=Rl*10"<-250) 

3250     Scexp=Scexp+250 

32G0   END  IF 

3270   IF  Rl<10"<-250)  THEN 

3280    R1=R1«10"250 

3290     Scexp-Scexp-250 

3300   END  IF 

3310  ! 

3320  I 

3330 

3340 

3350 

3360  ! 

3370   SUBEND 

3380  ! 

3390   !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3400   ! 

3410   I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

3420   ! 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 
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3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 
3710 
3720 
3730 
3740 


I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB  Disb(Lnax  ,Ka  ,Dy( «  )  ,Y<*  )) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM    THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  IRRE6ULAR  SPHERICAL 
REM    BESSEL  FUNCTIONS  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .   THE 
REM    CALCULATION  OF  THE  DERIVATIVE  USES  THE  ISB'S  OF  THE  SAME  AND 
REM    PREVIOUS  ORDER. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++  CALCULATING  THE  DERIVATIVES  OF  THE  ISB  +++++++++++++++++++++++++++ 

I 


CALL  Isb(Lnax  ,Ka.Y(*  )  ) 

RAD 

Dy<0)=-(-Y(0)  +  SIN(Ka)  )/Ka 

FOR    L=1    TO   Lriax 

Dy(L)«-((L+t  >/Ka  >*Y(L)  +  Y(L-1  ) 
NEXT    L 


!  CALCULATE    THE    ISB   VALUES. 

I 

I  ANGLE  UNITS  FOR  CALC'NS. 

!  THE  INITIAL  ISB  DERIVATIVE. 

I  THE  REMAINING  ISB  DERIVATIVES 


! 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 
SUBEND 

! 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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3750 
37B0 
3770 
3780 
3790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4030 
4040 


! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$ 

SUB  Isbdmax  ,Ka  ,Y(  *  )  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  CALCULATES  THE  IRREGULAR  SPHERICAL  BESSEL  FUNCTIONS 

REM    OF  AR6UMENT  KA  AND   ORDER  L  THROUGH  LMAX .    THE  CALCULATION  IS 

REM    ACCOMPLISHED  BY  FORWARD  RECURSION. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  CALCULATING  THE  ISB  ++++++++++++++++++++++++++++++++++++++++++++++ 


RAD 

Y(0)  =  -COS(Ka)/Ka 

Y<  1  )=(  Y(0)-SIN(Ka  )  )/Ka 

FOR  L=2  TO  Lmax 

Y(L)=((2«L-1  >/Ka)«Y(L-1  )-Y(L-2) 
NEXT  L 


ANGLE  UNITS  FOR  CALC'NS. 
INITIAL  ISB'S  CALCULATED 

!  REMAINING  ISB'S  CALC'D. 


j 

I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

SUBEND 

! 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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4050 
4060 
4070 
4080 
4090 
4100 
4!  10 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 


i  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Leg180(l_piax  ,Ka  ,P<  *  )  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  CALCULATES  THE  LEGENDRE  POLYNOMIALS  OF  ARGUMENT  -1  AND 

REM    ORDER  0  THROUGH  LMAX .   CALCULATION  IS  ACCOMPLISHED  BY  FORWARD 

REM    RECURSION. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

l 

I++++  CALCULATING  THE  LEGENDRE  POLYNOMIALS  +++++++++++++++++++++++++++++ 

l 
X=-1  !  ARGUMENT  OF  THE  LEG.  POL'S. 

I 
P(  0 )=1  !  INITIAL  LEG.  POL'S.  CALC'D. 

P<  1  )=X 
FOR  L=2  TO  Lnax  !  REMAINING  LEG.  POL'S.  CALC'D. 

P(L)-(2-(1/L))*X*P(L-1)-(1-(1/L)  )*P(L-2  ) 
NEXT  L 

! 

I 

!++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 
SUBEND 

I 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
I 
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4340 
4350 

4360 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
4460 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4610 
4620 
4630 
4640 
4650 
4660 
4670 
4680 
4690 
4700 
4710 
4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
4860 
4870 
4880 


I 
!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Recplt(Nbs(«),Kapnin,Kanax,Stp> 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM  THIS  MODULE  MAPS  THE  CALCULATED  SCATTERING  MODULUS  ARRAY  SMC)  ON 

REM  A  SEMI-LOG  PLOT. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Xgumax 
REAL  Ygumax 


THE  MAXIMUM  ABSCISSA  VALUE  IN 
GRAPHICS  DISPLAY  UNITS. 

THE  MAXIMUM  ORDINATE  VALUE  IN 
GRAPHICS  DISPLAY  UNITS. 


+++++++++++++++++T+++++++++++++++++++++++++++++++++++++++++++++++++++++ 

+  +  +  +  INITIALIZING  THE  PLOTTER  +++++++++++++++++++++++++++++++++++++++++ 

i 

GINIT 

PLOTTER  IS  CRT, "INTERNAL" 

GRAPHICS  ON 

GCLEAR 

I 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  GRAPH  CONSTRUCTION  AND  LABELLING  +++++++++++++++++++++++++++++++++ 


Xgunax=100«MAX<  1  .RATIO) 
Ygumax=100*MAX(  1  ,1/RATIO) 

CSIZE  4 
LINE  TYPE  1 


LORG  4 

MOVE  Xgumax/2  ,0 

LABEL  "K«A" 

DEG 

LDIR  90 

LORG  6 

MOVE  0  ,Ygumax/2 

LABEL  "NORM'D  BCKSCTRG  CROSS  SEC'N' 

LDIR  0 

VIEWPORT  10  ,Xgunax-5  ,10  .Ygunax-5 


UINDOLI  Kanin  .Kamax  ,0  ,  1  .2 

LINE  TYPE  4 

GRID  Stp«  10  ,  .  1  .Kamin  ,0  ,61  ,  10  .2 


! 

!  **  LABELLING  THE  ABSCISSA  AXIS  ** 


I 

!  *«  LABELLIN6  THE  ORDINATE  AXIS  ** 


!  *»  GRAPH  CONSTRUCTION  *« 

i  DEFINES  THE  TOTAL  PLOTTING  AREA 

i  IN  GRAPH  UNITS. 
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4890 
4900 
491  0 
4920 
4930 
4940 
4950 
49G0 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
51  10 
5120 
5130 
5140 
5150 
51G0 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 
5260 
5270 
5280 
5290 
5300 
5310 
5320 


! 

!  **  LABELLING  THE  ABSCISSA  SCALE 


«*  LABELLIN6  THE  ORDINATE  SCALE  ** 


CLIP  OFF 
CSIZE  3 
LINE  TYPE  1 

LORG  6 

FOR    I=Kamin    TO    Kamax    STEP    Stp*2< 

MO'JE  I  ,0 

LABEL  I 
NEXT  I 

LORG  8 

FOR  1=0  TO  1  STEP  .5 

MOVE  Kami n,  I 

LABEL  I 
NEXT  I 

CLIP  ON 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

+  +  +  +  PLOTTING  THE  CALCULATED  NBCS  +++++++++++++++++++++++++++++++++++++ 

■  ! 
FOR  I=Kamin+Stp  TO  Kamax  STEP  Stp 

PLOT  I  ,Nbs< ( I-Kamin)/Stp ) 
NEXT  I 

! 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  PRINTING  THE  PLOT  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 
DUMP  GRAPHICS  #701 

! 
I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 
!++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 
GRAPHICS  OFF 
SUBEND 

I 
I +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  -t +  +  +  +++  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 
! 
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APPFND1X    n 
ACOUSTIC   ARRAY  EFFICIENCY    PROGRAM 

This  program  was  used  to  analyze  data  and  determine  the  calibration 
product  Er  Et.    The  output  of  this  program  was  also  used  to  construct 

Figures  13,  14,  and  15. 

The  power  supply  mentioned  in  the  Gain30  and  Gain35  subprograms 
was  the  HP33 1 4A  function  generator.  All  data  contained  in  the  data 
subprograms  are  repeated  in  Appendix  E.    Data  units  are  available  in 
Appendix  E  with  their  respective  data. 


20  ! 

30  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

40  REM  ACOUSTIC  ARRAY  EFFICIENCY 

50  REM 

G0  REM      THIS  PROGRAM  CALCULATES  THE  EFFICIENCY  AND  CENTERLINE  GAIN  OF 

70  REM    AN  ACOUSTIC  ARRAY  BASED  ON  THE  RETURN  FROM  AN  ACOUSTICALLY  HARD 

80  REM    SPHERE  OF  DIAMETER  D  ENSONIFIED  BY  AN  INCIDENT  PLANE  WAVE  OF 

90  REM    FREQUENCY  F. 

100  REM 

1 30  REM   ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

140  ! 

150  !++++  VARIABLE  DECLARATIONS  AND  DEFINITIONS  ++++++++++++++++++++++++++++ 

160 

170 

180  REAL  Vmsr 

190 

200 

210 

220  REAL  VnSn 


230 
240 

250 


MEAN  SQUARE  VOLTAGE  RECEIVED  FROM 
THE  TARGET  SPHERE  AND  MEASURED 
AFTER  PRE-AMPLIFICATION. 


MEAN  SQUARE  VOLTAGE  RECEIVED  IN 
THE  ABSCENCE  OF  A  TARGET  AND 
MEASURED  AFTER  PRE-AMP. 
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260 

REAL 

Vmst 

270 

280 

REAL 

R 

290 

300 

310 

320 

330 

REAL 

Et 

340 

350 

3G0 

370 

REAL 

Go 

380 

390 

400 

REAL 

Etgo 

410 

420 

REAL 

At 

430 

440 

450 

REAL 

C  ,Cmoi st 

460 

470 

REAL 

Efreq 

480 

490 

500 

REAL 

K 

510 

520 

REAL 

A 

530 

540 

REAL 

Ka 

550 

560 

570 

580 

REAL 

Aa 

590 

600 

REAL 

Gpre 

610 

620 

REAL 

Nbs 

630 

640 

650 

660 

REAL 

Atgt 

670 

680 

690 

REAL 

Sasph 

700 

710 

720 

730 

740 

750 

REAL 

Prat 

760 

770 

780 

MEAN  SQUARE  VOLTAGE  TRANSMITTED. 

THE  RANGE  FROM  THE  PLANE  OF  THE 
ARRAY  SPEAKER  DIAPHRAGMS  TO  THE 
CENTER  OF  THE  TARGET  SPHERE  AND  TO 
THE  SPECTRUM  ANALYZER  MICROPHONE. 

THE  EFFICIENCY  OF  CONVERSION 
FROM  TRANSMITTED  ELECTRICAL  POUER 
TO  TRANSMITTED  ACOUSTICAL  POUER. 

THE  CENTERLINE  GAIN  OF  THE  ACOUSTIC 
ARRAY. 

ET*GO 

THE  ATTENUATION  COEFFICIENT, 
"ALPHA" . 

THE  SPEED  OF  SOUND  IN  HUMID  AIR. 

THE  SELECTED  FREQUENCY  OF  THE 
ACOUSTIC  ARRAY. 

THE  WAVE  NUMBER,  2*PI*EFREQ/C . 

THE  RADIUS  OF  THE  TARGET  SPHERE. 

K*A.  THE  ARGUMENT  OF  THE  REGULAR 
AND  IRREGULAR  SPHERICAL  BESSEL 
FUNCTIONS  AND  THEIR  DERIVATIVES. 

THE  APERATURE  AREA  OF  THE  ARRAY. 

THE  PRE-AMF'S  ELECTRICAL  GAIN. 

THE  VARIABLE  REPRESENTING  THE 
NORMALIZED,  BACKSCATTERED  CROSS 
SECTION  (NBCS). 

CROSS  SECTIONAL  AREA  OF  THE 
TARGET  SPHERE. 

SURFACE  AREA  OF  AN  IMAGINARY 
SPHERE  CENTERED  ON  THE  ARRAY 
AND  HAVING  A  RADIUS  EQUAL  TO 
THE  RANGE  TO  THE  CENTER  OF  THE 
TARGET  SPHERE. 

THE  RATIO  OF  RECEIVED  TARGET 
ELECTRICAL  POUER  TO  THE 
TRANSMITTED  ELECTRICAL  POUER. 


0 


790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

10G0 

1070 

1080 

1090 

1  100 

1  1  10 

1  120 

1  130 

1  140 

1  150 

1  1S0 

1  170 

1  180 

1  190 

1200 

1210 

1220 

1230 


REAL  Er 


!  THE  EFFICIENCY  OF  CONVERSION 
i  FROM  RECEIVED  ACOUSTICAL  POWER 
i  TO  RECEIVED  ELECTRICAL  POWER. 


I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++  MAIN  PROGRAM  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

l 

I 
CALL  Imt(Vmsr  .Vmsn.Vnst  .R.Etgo.At  ,A,Ka  ,Aa  ,Gpre  ) 

l  INPUT  AND  CALCULATE  REQUIRED 

I  PARAMETERS. 


CALL  Bcksctr(Ka.Nbs  ) 
PRINT  "NBS=";Nbs 


i 

!  CALCULATE  THE  NBCS 


Atgt=PI»A*A 
PRINT    "ATGT="iAtgt 
Sasph=4»PI»R»R 

Prat  =  ( Vmsr — Vmsn  )/( Gpre*Gpre*Vpis t  ) 

Er  =  (Prat  )♦ ( Sasph»Sasph/Etgo  )»( EXP( 2*At *R  )/( Aa* Atgt ) )/Nbs 

! 

PRINT  "Er=";Er 
PRINT  "EtGo=" iEtgo 
PRINT  "Go=" jEtgo/Er 
PRINT  "EtEr=";Er*Er 


+++++++++++++++++++-r+++f +++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MAIN  PROGRAM  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++ 


PRINT  CHR$( 12  ) 

PRINTER  IS  1 

PRINT 

PRINT  "END" 

PRINT 

END 


!  ADVISES  THE  USER  OF  THE  MAIN 
I    PROGRAM'S  CONCLUSION. 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$s$$$$$$$$$$s$$$$$$$s$$$$$$$$$$$$$$$$$ 


1240 
1250 
I2G0 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
13G0 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
I4G0 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
15G0 
1570 
1580 
1590 
1G00 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 


1$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
SUB  Init(  Vmsr .Unsn  ,Vnst  ,R  .Etgo  ,At ,A  ,Ka  ,Aa  .Gpre  ) 


REM 
REM 
REM 
REM 
I 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
THIS  MODULE  DETERMINES  THE  VALUES  OF  THE  PARAMETERS  REQUIRED  FOR 
THE  CALCULATION  OF  Er. 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


REAL  Ch 

REAL  Vr(60) 

REAL  J  ,Jt  .Jn.Jtn 

REAL  Tc 
REAL  Patn 

REAL  Rh 

REAL  Tgttm 

REAL  Dia 
REAL  Nuspkr 
REAL  Diaspkr 
REAL  Vt<60> 
REAL  Noise(60) 
REAL  Tnoise(60) 

REAL  Ten 
REAL  Rhn 


A  CONTROL  VARIABLE  USED  FOR 
CONTROLLING  SELECTION  OF  THE 
CORRECT  GAIN  DATA  SETS. 

THE  ARRAY  REPRESENTING  THE 
MEASURED  TARGET  RETURN  VOLTAGES. 

THE  NUMBER  OF  ARRAY  ELEMENTS  MINUS 
ONE  AND  A  LOOP  CONTROL  VARIABLE. 

THE  AMBIENT  TEMPERATURE  IN  DEG.  C. 

THE  ATMOSPHERIC  PRESSURE  MEASURED 
IN  MILLIBARS. 

THE  RELATIVE  HUMIDITY  MEASURED  IN 
PERCENT. 

THE  MEASURED  TIME  FOR  THE  ACOUSTIC 
SIGNAL  TO  MAKE  THE  ROUND  TRIP 
BETWEEN  THE  ARRAY  AND  THE  TARGET. 

I  THE  MEASURED  DIAMETER  OF  THE 
TAR6ET  SPHERE. 

THE  NUMBER  OF  SPEAKERS  IN  THE 
ARRAY. 

THE  AVERAGE  DIAMETER  OF  AN  ARRAY 
SPEAKER. 

THE  ARRAY  REPRESENTING  THE 
MEASURED  TRANSMISSION  VOLTAGES. 

THE  ARRAY  REPRESENTING  THE 
MEASURED  NOISE  RETURN  VOLTAGES. 

THE  ARRAY  REPRESENTING  THE 
MEASURED  TRANSMISSION  VOLTAGES  FOR 
NOISE  DETERMINATION. 

Tc  FOR  NOISE  DETERMINATION. 

Rh  FOR  NOISE  DETERMINATION. 
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1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
21  10 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 


REAL  Patnn 
REAL  Efreqn 
REAL  Cn 
REAL  Atn 
REAL  Unstn 
REAL  Rspkr 


Patn  FOR  NOISE  DETERMINATION. 

E  f reqn  =  Ef req . 

C  FOR  NOISE  DETERMINATION. 

At  FOR  NOISE  DETERMINATION. 

Unst  FOR  NOISE  DETERMINATION. 

THE  AVERAGE  RADIUS  OF  AN  ARRAY 
SPEAKER. 


i ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I++++  DETERMINING  SYSTEM  PARAMETERS  +++++++++++++++++++++++++++++++++++++ 

I 
! 
PRINTER  IS  1 

PRINT  "WHICH  TARGET  DATA  SET  "; 
PRINT  "DO  YOU  WISH  ANALYZED?" 
PRINTER  IS  701 
INPUT  Sphere 

PRINT  "TARGET  NO.  "(Sphere 
PRINT 

! 
IF  Sphere=1  THEN 

Ch=3.0 

CALL  Gaino(Ch,Etgo  ) 

CALL  Sphere!  ( \)r(  *  )  ,  J  ,Tc  .Patrn  ,Rh  ,Tgt  tm  ,Dia  .Efreq  .Nuspkr  .Diaspkr  ,Gpre  ) 

CALL  Trans30( Jt  ,Vt(«  )) 

CALL  Noise30(Jn,Jtn,Noise<»),Tnoi5e(*),Tcn,Rhn,Patnn  .Efreqn  ) 
END  IF 

! 
IF  Sphere=2  THEN 

Ch=3.5 

CALL  Gaino(Ch  .Etgo  ) 

CALL  Sphere2(l>r(  *)1J,Tc,Patn,Rh,Tgttn,Dia  .Efreq  .Nuspkr  .Diaspkr  .Gpre  ) 

CALL  Trans35( Jt  .Ut( »  )) 

CALL  Nolse35( Jn  ,Jtn,Noise<*)  ,Tnoise< *  )  .Ten  ,Rhn  ,  Patnn  .Efreqn  ) 
END  IF 

! 
IF  Sphere=3  THEN 

Ch=3.S 

CALL  Gaino(Ch  .Etgo  ) 

CALL  Sphere3( Vr(  *  )  , J ,Tc  .Patn  ,Rh  .Tgttn  ,Dia  .Efreq  .Nuspkr .Diaspkr .Gpre  ) 

CALL  Trans35( Jt ,Vt(»  )  ) 

CALL  Noise35(Jn  ,Jtn,Noise( *   )  ,Tnoise( *  )  .Ten  ,Rhn  , Patnn  .Efreqn  ) 
END  IF 

l 
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2300  IF  Sphere-4  THEN 

2310     Ch=3.5 

2320     CALL  Gaino(Ch.Etgo) 

2330      CALL  Sphere4(<Jr( O.J.Tc.Patn.Rh.Tgttm.Dia  .Efreq  .Nuspkr .Diaspkr  ,Gpre  > 

2340     CALL  Trans35( Jt  ,Vt<*  )  ) 

2350      CALL  Noise35( Jn,Jtn,Noise<*)  ,Tnoise( *  )  ,Tcn  ,Rhn  .Patmn .Efreqn ) 

2360  END  IF 

2370  ! 

2380  CALL  Voltsms<Vr(*  )  ,J ,Unsr> 

2390  PRINT  ,'UMSR="  iUmsr 

2400  CALL  Atten(Tcn,Rhn ,Patmn  ,Cn .Atn .Efreqn ) 

2410  CALL  Uoltsms(Tnoise< ♦  )  .Jtn.Umstn) 

2420  CALL  Voltsms(Noise( *  )  .Jn.Umsn  ) 

2430  CALL  Voltsms('Jt<  *  )  ,Jt  ,Unst  ) 

2440  PRINT  "WtST="fVi>ist 

2450  ! 

2460  A=Dia/2 

2470  PRINT  "A=";A 

2480  CALL  Range(Tgttm,A ,Tc ,Rh,Patn,C  ,At  , Efreq  ,R) 

2490  PRINT  "C=";C 

2500  PRINT  "AT="iAt 

2510  PRINT  "ATTEN=";EXP(At* R) 

2520  PRINT  "R=";R 

2530  Vmsn=Umsn*EXP(2*R*< Atn-At  )  )*fms t /<Jns t n 

2540  PRINT  "VMSN=" ;Vnsn 

2550  Ka=2*PI*Efreq»A/C 

2560  PRINT  "KA=";Ka 

2570  Rspkr=Diaspkr/2 

2580  Aa  =  Nuspkr*PI*Rspkr-*Rspkr 

2590  PRINT  "AA=";Aa 

2600  ! 

2610  ! 

2620  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

2630  ! 

2640  I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

2650  I 

2660  SUEEND 

2670  ! 

2680  I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

2690  I 

2700  i $2$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$S$S$ $$$$$$$$$$ $$$$$$$ 

2710  ! 


2720 

2730 

2740 

2750 

27G0 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

28G0 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

30G0 

3070 

3080 

3090 

3100 

31  10 

3120 

31  30 

3140 

3150 


l $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$£$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Gaino<Ch,Etgo) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  DETERMINES  THE  UALUE  OF  EtGo. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

'++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Urns 


REAL  Ir 


REAL  Za 


ROOT  MEAN  SQUARE  UOLTAGE  MEASURED 
IN  GAIN  DETERMINATION. 

MEASURED  INTENSITY  AT  R; 
Ptgt/Atgt. 

ACOUSTIC  ARRAY  IMPEDANCE. 


! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

I++++  DETERMINATION  OF  EtGo  ++++++++++++++++++++++++++++++++++++++++++++ 

| 
I 
IF  Ch  =  3.0  THEN  CALL  Gain30(  Urms  ,Tc  ,Rh  ,R  ,Patm  ,  In  ,Ef  req  ,Za  ) 
IF  Ch=3.5  THEN  CALL  Gai n35( Urns  ,Tc  ,Rh ,R  ,Patn , Ir ,Ef req ,Za ) 

OBTAINING  UALUES  FROM  APPROPRIATE 
DATA. 

CALL  Atten(Tc ,Rh  ,Patn  ,C . At  .Efreq  )    I  CALCULATING  AT. 

Pt=Urms*Vrms/Za 

Etgo=Ir»4*PI*R»R*EXP(At*R)/Pt        i  CALCULATING  EtGo. 


! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

i 

!++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 


I 


SUBEND 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ s$$$$s$$$$$$s $$$$$$$$$$$ 


31G0 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
32B0 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
33G0 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
35G0 
3570 
3580 
3590 
3G00 
3610 
3620 
3550 
3640 
3G50 
3560 


l$$$S$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$S$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Voltsms(Volt< *  )  ,J  ,Vns  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  DETERMINES  THE  MEAN  SQUARE  VALUE  OF  THE  INPUT  VOLTAGE 

REM   ARRAY. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


IREAL  VOLT<»  ) 

REAL  Volts 
REAL  Vavg 
REAL  V 

IREAL  Vns 


THE  ARRAY  REPRESENTING  THE  INPUT 
VOLTAGE  ARRAY. 

THE  SUM  OF  VOLT(»)'S  ELEMENTS. 

THE  AVERAGE  VALUE  OF  VOLT(»). 

THE  DIFFERENCE  BETWEEN  VOLT( I  ) 
AND  VAVG. 

THE  MEAN  SQUARE  VALUE  OF  VOLT(*). 


I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
i 

1 ++++  DETERMINATION  OF  VMS  +++++++++++++++++++++++++++++++++++++++++++++ 

! 

i 

Volt  5=0 

FOR  1=0  TO  J 

Volts=Volts+Volt(  I  ) 
NEXT  I 

Vavg=Volts/( J+1 ) 
FOR  1=0  TO  J 

V=Volt( I >-Vavg 

Vms=( V*V >+Vns 
NEXT    I 
Vns=Vns/< J+1 ) 


i 

I+  +  +  +  MODULE  CONCLUSION  +  +  +  +  +  +  +  +++++++++++++++++++  +  +  +++  ++++++  +  +  +  +  +  +  ++-'-+  + 

l 

SUBEND 

i 

i 
l$SS$$$S$S$5$$$$$$$$$$S$$$$$S3SS$$S$S$$SSSS$S$$SSS$SSS$S$SSSS$S$$SS$S$$$ 
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3670  i 

3680  i $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$ 

3E90  SUB  Range(  Tgttrn  ,A  ,Tc  ,Rh  .Patm  ,C  .At  .Efreq  ,R  ) 

3700  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3710  REM    THIS  MODULE  DETERMINES  THE  AUERAGE  PATH  ATTENUATION  AND  THE  RAN6E 

3720  REM    BETWEEN  THE  ARRAY  AND  THE  TARGET  CENTER. 

3730  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3740  ! 

3750  I++++  DETERMINATION  OF  R  AND  AT  ++++++++++++++++++++++++++++++++++++++++ 

3760  ! 

3770  ! 

3780  Tgttn=Tgttm/1000 

3790  CALL  Atten(Tc .Rh.Patm.C, At  , Efreq  )    !  AT  AND  C  DETERMINED. 

3800  R=Tgttn*C/2+A  !  R  DETERMINED. 

3810  ! 

3820  ! 

3830  • +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3840  ! 

3850  I++++    MODULE    CONCLUSION    ++++++++++++++++++++++++++++++++++++++++++++++++ 

3860  ! 

3870  SUBEND 

3880  i 

3890  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3900  i 

3920  i 
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5930 
3940 
3950 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4030 
4040 
4050 
4060 
4070 
4080 
4090 
4)00 
41  10 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 


! 

!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$ 

SUB  Atten( Tc  ,Rh .Patm  .Cmoist  .At  .Efreq  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  CALCULATES  THE  ATTENUATION  USING  RELATIVE  HUMIDITY 

REM    IN  PERCENT.  ATMOSPHERIC  PRESSURE  IN  MILLIBARS.  AND  TEMPERATURE 

REM    IN  DEGREES  CENTIGRADE. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

I++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Tk 
REAL  Es 
REAL  Pratio 

REAL  Cdry 
REAL  H 
REAL  Past 
REAL  Tast 
REAL  Fm 
REAL  Amax 
REAL  Frat 
REAL  F2 
REAL  Acl  .Anol 


THE  AMBIENT  TEMPERATURE  IN  DEG.  K. 

THE  SATURATION  VAPOR  PRESSURE. 

THE  RATIO  OF  WATER  PRESSURE  IN 
MILLIBARS  TO  ATMOSPHERIC  PRESSURE 
IN  MILLIBARS. 

THE  SPEED  OF  SOUND  IN  DRY  AIR. 


P  ASTERISK. 

T  ASTERISK. 

MAX.  FREQ. 

MAX  EXPECTED  ATTENUATION. 

EFREQ/FM 

FRAT'FRAT 

ATTENUATION  EXPONENT  COMPONENTS. 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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4340 
4350 
43G0 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
44G0 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
45G0 
4570 
4580 
4590 
4G00 
4G10 
4620 
4630 
4640 
4650 
46G0 
4670 
4680 
4690 
4700 
4710 


I++++  CALCULATING  THE  ATTENUATION  ++++++++++++++++++++++++++++++++++++++ 


Tk=Tc+273. 15 


Rh=Rh/!00 

Es=10"(9.4-2353/Tk ) 
Pratio=Es»Rh/(Patn-Es*(  1+Rh  )  ) 

Cdry=20.05*SQR(Tk ) 
Cmoist=Cdry*(  1  +  .  14*Pratio ) 


CONVERSION  OF  TEMP  IN  CELSIUS  TO 
TEMP  IN  KELVIN. 

CONVERSION  OF  RH  IN  PCT.  TO  PRATI0, 


i  CALCULATION  OF  THE  SPEED  OF  SOUND 


l 


H=Pratio*100 
Past=Patm/1014 
Tast  =  (  1  .8*Tc+492  )/519 

Fm  =  ( 10+6600»H+44400»H*H  >*Past/( Tast" . 8 ) 

! 

Anax  =  .0078*Fm*<Tast/>(-2.5  )  )*EXP<7.77»<  1-1  /Tast  )  ) 

i 

Frat=Ef req/Fm 

F2=Frat*Frat 

I 

Amol=Amax«SQR(  .18*. 1 8*F2  +  < 2*F2/( 1+F2  )  T2)/304.8 
Acl  =  l  .74*10"(-10  )*Efreq*Efreq 
At  =  (Anol+AcI  )*LOG< 10)/10 

I 

I 

I 

i ++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

i 

SUBEND 

i 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

l $$$$$$$$$$$$$$$$$$$$$$$ $$$$$j$$$$$$$$$ $$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$ 


19 


4720 
4730 
4740 
4750 
4760 
4770 
4780 
4790 
4800 
4810 
4820 
4830 
4840 
4850 
48B0 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
51  10 
5120 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 


SUB  Bcksctr(Ka  ,Nbs  ) 


REM 
REM 
REM 
REM 

I 

I++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
THIS  MODULE  CALCULATES  THE  BACKSCATTERED  CROSS  SECTION  NORMALIZED 
BY  THE  CROSS  SECTIONAL  AREA  OF  THE  SPHERE. 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


REAL  L 


REAL  Lmax 

REAL  J (51  ) 

REAL  Dj<  51  ) 
REAL  Y( 51  ) 

REAL  Dy(51  ) 

REAL  P<  51  ) 

REAL  Fac 
REAL  Rsm 

REAL  Ism 


ORDER  OF  THE  FUNCTION  CALCULATED 
IN  A  SUBROUTINE.   IT  IS  USED 
AS  A  LOOP  INDEX  AND  TO  DESIGNATE 
THE  ARRAY  ELEMENT  FOR  THE  FUNCTION 
OF  CORRESPONDING  ORDER. 

MAX.  SIGNIFICANT  ORDER  OF  THE  SUM 
DETERMINING  THE  NORMALIZED 
BACKSCATTERING  CROSS  SECTION 
<NBCS). 

ARRAY  FOR  THE  REGULAR  SPHERICAL 
BESSEL  FUNCTIONS  <RSB)  OF 
ARGUMENT  KA  AND  ORDER  0  THROUGH 
LMAX. 

ARRAY  FOR  THE  DERIVATIVES  Of1  THE 
RSB  (THE  DRSB). 

ARRAY  FOR  THE  IRREGULAR  SPHERICAL 
BESSEL  FUNCTIONS  (ISB)  OF 
ARGUMENT  KA  AND  ORDER  0  THROUGH 
LMAX. 

ARRAY  FOR  THE  DERIVATIVES  OF  THE 
ISB  (THE  DISB). 

ARRAY  FOR  THE  LEGENDRE 
POLYNOMIALS  OF  ARGUMENT  COSU80). 

FACTOR  COMMON  TO  RSM  AND  ISM. 

THE  REAL  COMPONENT  OF  THE 
SCATTERING  MODULUS,  SM. 

THE  IMAGINARY  COMPONENT  OF  THE 
SM. 
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5210 
5220 
5230 
5240 
5250 
52B0 
5270 
5280 
5290 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
5370 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
54G0 
5470 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
55B0 
5570 
5580 
5590 
5600 
5610 
5620 


++++  NOTE  ++++ 

THE  ARRAY  DIMENSIONS  (30+2  1)  NEED 

TO  BE  INCREASED  TO  ACCOMODATE 

KA>30. 


I +++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

!++++  CALCULATING  THE  NBCS  +++++++++++++++++++++++++++++++++++++++++++++ 


I 


Lmax=INT(Ka+21  ) 

CALL  Drsb(Lnax  ,Ka  ,Dj< *  )  ,J< * >  ) 

CALL  Disb(Lmax ,Ka  ,Dy(*  )  ,Y(*  )) 

CALL  Leg180(Lnax  ,Ka  ,P( ♦  )  ) 


CALCULATE  THE  DERIVATIVES  OF  THE 
REGULAR  SPHERICAL  BESSEL  FUNCTIONS 
OF  ORDER  0  THROUGH  LMAX. 
CALCULATE  THE  DERIVATIVES  OF  THE 
IRREGULAR  SPHERICAL  BESSEL  FUNC'NS 
OF  ORDER  0  THROUGH  LMAX. 

CALCULATE  THE  ARRAY  P(*). 


Rsm=0 
lsm  =  0 
FOR  L=0  TO  Lnax 

Fac  =  (Dj(L)*(2*L+1  >*P(L))/(Dj(L  )*Dj  (  L  >+Oy(  L  )*Dy  (  L  )  ) 

Rsm=Rsm+<Fac»Dy(L ) ) 

Isn=Ism+(Fac»Dj(L ) ) 
NEXT  L 


CALCULATING  THE  REAL  AND  IMAGINARY 
COMPONENTS  OF  THE  SM. 


Nbs  =  4*(Rsn*Rsrn+l5rn*Isn)/<Ka*Ka) 


CALCULATING    THE    NBCS 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

+  +  +  +    MODULE    CONCLUSION    ++++++++++++++++++++++++++++++++++++++++++++++++ 


l 


SUBEND 


!+++++++  +  +  +++++.+++++++  +++++++++  +  ++++++++  +++  ++++  +  +++  +  +  +++  +  ++++++  +++  +  +  +++  + 
i 

l$$$$$$$$$$$$S$$$$$$$$$$$$$$S$$$$$$$$$$$$$$$$$S$$$$$$$$$$$$S$S$$$$$$$$$S 
I 
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5630 
5G40 
5650 
5560 
5G70 
5680 
5690 
5700 
5710 
5720 
5730 
5740 
5750 
57G0 
5770 
5780 
5790 
5800 
5810 
5820 
5830 
5840 
5850 
58G0 
5870 
5880 
5890 
5900 
5910 
5920 
5930 
5940 


!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$S$$$$$S$$$S$$$$$$$$S$S$$$$$$$$$$$$ 
SUB  Drsb(Lmax  ,Ka  ,Dj(  «  )  ,J(  *  )  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM    THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  REGULAR  SPHERICAL 
REM    BESSEL  FUNCTIONS  ( DRSB  )  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX . 
REM    THE  CALCULATION  OF  THE  DERIVATIVE  INVOLVES  THE  RSB'S  OF  THE  SAME 
REM    AND  PREVIOUS  ORDER. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I++++  CALCULATING  THE  DERIVATIVES  OF  THE  RSB  +++++++++++++++++++++++++++ 

! 


CALL  RsbCLmax ,Ka ,J<  * )  ) 

RAD 

Dj(0)  =  <-J(0)+COS(Ka  )  >/Ka 

FOR  L=1  TO  Lmax 

Dj(L)  =  -((L+l  )/Ka)*J(L)+J(L-1  > 

NEXT  L 


CALCULATE  THE  RSB  VALUES  JO. 

ANGLE  UNITS  FOR  DRSB  CALC'NS. 
CALCULATE  THE  INITIAL  DRSB  DJ(0) 
CALCULATE  THE  REMAINING  DRSB. 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

+  +  +  +  MODULE  CONCLUSION  +  +++++  ++++++  ++++++  +  +  +  +++-r++++++++++++++-f+++++++  + 

i 

SUBEND 

I 

i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I 
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5950 
59G0 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6210 
6220 
6230 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 
6350 
6360 
6370 
6380 


! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Rsbdnax  ,Ka  ,J(  *  >  ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

REM    THIS  MODULE  CALCULATES  THE  REGULAR  SPHERICAL  BESSEL  FUNCTIONS 

REM    (RSB)  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .    THE 

REM    CALCULATION  USES  THE  CONTINUED  FRACTION  APPROACH. 

REM    ++++++++++++++++++++++++++++++++++++++ +++++++++1-++++++++ ++++++++++ 

I 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Rl 
REAL  Numr 

REAL  Denr 

REAL  Scexp 

REAL  Nflag 
REAL  Dflag 
REAL  Bnn 

REAL  Bine 


THE  EVOLVING  RATIO  OF  J(0)/J(L). 

THE  EVOLVING  NUMERATOR  FACTOR  USED 
IN  CALCULATION  OF  RL. 

THE  EVOLVING  DENOMINATOR  FACTOR 
USED  IN  CALCULATION  OF  RL. 

EVOLVING  SCALING  EXPONENT  FOR  VERY 
LARGE  OR  VERY  SMALL  VALUES  OF  RL. 

0  FLAG  FOR  NUMR. 

0  FLAG  FOR  DENR. 

EVOLVING  TERM  USED  IN  THE 
CALCULATION  OF  NUMR  AND  DENR. 

INCREMENT  USED  IN  THE  EVOLUTION 
OF  BNM. 


+++++++++++++++++++++++++++ +++++++++++-T+++ ++++++++++++++++ +++++++++++++ 

++++  CALCULATING  THE  RSB  ++++++++++++++++++++++++++++++++++++++++++++++ 


RAD 


J(0)  =  (SIN(Ka  )  )/Ka 


ANGLE  UNITS  FOR  CALC'NS. 

♦*»  CALCULATE  THE  INITIAL  *** 
♦*♦  RSB,  J<0). 
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6390 

6400 

FOR  L=l  TO  Lmax 

G4I0 

Rl=l .0 

G420 

Numr=0 . 

6430 

Denr=0. 

G440 

Scexp=0 . 

G450 

Nflag=0. 

6460 

Dflag=0. 

6470 

Bnm=l .0/Ka 

6480 

Binc=2.0/Ka 

6490 

6500 

6510 

FOR  Ordcnt=l  TO  L 

6520 

Bnm=Bnm+Bi  nc 

6530 

IF  Nflag=0  THEN 

6540 

Numr=Bnn-Numr 

6550 

IF  Numr<>0  THEN 

6560 

Rl=Rl»Nunr 

6570 

Nunr= 1 . 0/Nunr 

6580 

6590 

IF  ABS(R1  )>10 

6600 

ELSE 

6610 

Nf  lag=l 

6620 

END  IF 

6630 

ELSE 

6640 

Nflag=0 

6650 

END  IF 

6660 

NEXT  Ordcnt 

6670 

6680 

6690 

6700 

6710 

WHILE  NumrODenr  OR 

6720 

Bnm=Bnm+Binc 

6730 

6740 

IF  Nflag=0  THEN 

6750 

Numr=Bnn-Nunr 

6760 

IF  Nunr<>0  THEN 

6770 

Rl=RI*Numr 

6780 

Numr  = 1 . 0/Nunr 

6790 

6800 

IF  ABS(R1  )>10 

6810 

ELSE 

6820 

Nfiag=1 

6830 

END  IF 

6840 

ELSE 

6850 

Nflag-0 

6860 

END  IF 

6870 

6880 

IF  Dflag=0  THEN 

6890 

Denr=Bnn-Denr 

6900 

IF  Denr<>0  THEN 

6910 

R1=R1/Denr 

6920 

Denr= 1 . 0/Denr 

•**  CALCULATING  THE  RATIO 

*  #  * 

•**    J(0)/J(L). 

#  #  # 

INITIALIZING  RL. 

INITIALIZING  NUMR. 

INITIALIZING  DENR. 

INITIALIZING  SCEXP. 

INITIALIZING  NFLAG. 

INITIALIZING  DFLAG. 

INITIALIZING  BNM. 

CALCULATE  BINC. 

•  CALCULATE  THE  L  "NAKED" 

♦  NUMR  FACTORS  OF  J(0)/J(L) 
I  INCREMENT  BNM. 

I  CHECK  NUMR=0  FLAG  NOT  SET; 
I  UPDATE  NUMR  IF  NUMRO0. 


I  UPDATE  RL  IF  UPDATED  NUMRO0. 
!  PREPARE  NUMR  FOR  NEXT  EVOLUTION. 
!  CHECK  SCALING  OF  UPDATED  RL. 
'250  THEN  CALL  Scale( Rl  .Scexp  ) 


I  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 
I  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


!  •  CALCULATE  THE  REMAINING  NUMR  • 

!  «  AND  DENR  TERMS  OF  J<0)/J(L).  » 

!  NOTE  WHEN  NUMR=DENR  THEN 

!  RL  =  J(0)/J(L). 
Nflag=l  OR  Dflag=l 

!  INCREMENT-  BNM. 

I 

!  CHECK  NUMR=0  FLAG  NOT  SET; 

!  UPDATE  NUMR  IF  NUMRO0. 

UPDATE  RL  IF  UPDATED  NUMRO0. 
PREPARE  NUMR  FOR  NEXT  EVOLUTION. 
CHECK  SCALING  OF  UPDATED  RL. 
'250  THEN  CALL  Scale( Ri  , Scexp  ) 


l  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 
!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


CHECK  DENR=0  FLAG  NOT  SET; 
UPDATE  DENR  IF  DENRO0. 

i  UPDATE  RL  IF  UPDATED  DENRO0. 
I  PREPARE  DENR  FOR  NEXT  EVOLUTION. 
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6930 
6940 
6950 
6960 
6970 
6980 
6990 
7000 
7010 
7020 
7030 
7040 
7050 
7060 
7070 
7080 
7090 
7100 
71  10 
7120 
7130 
7140 
7150 
7160 
7170 
7180 
7190 


i  CHECK  SCALING  OF  UPDATED  RL, 
IF  ABS(R1  )>10A250  THEN  CALL  Scale( Rl  .Scexp  ) 

ELSE 


Dflag=1 
END  IF 

ELSE 

Dflag=0 

END  IF 

END  WHILE 


J(L)  =  J<0)/<Rl»l0A(Scexp  )  ) 
NEXT  L 


i  IF  DFLAGO0  THE  PREVIOUS  STEPS  ARE 
!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


! 

i  ***    J(L)  DETERMINED.  *♦♦ 
I 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

| 

SUBEND 

I 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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7200  ' 

72 10  i$$$$$$$$$$$$$$$$$$$$$$$$$S$$S$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

7220  SUB  Scale(Rl  .Scexp  ) 

7230  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7240  REM    THIS  MODULE  SCALES  RL  TO  MAINTAIN  NUMERICAL  ACCURACY  IN  THE 

7250  REM    ONGOING  CALCULATION  OF  RL . 

72 60  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7270  ! 

7280  I++++  SCALING  RL  +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7290  ! 

7300  ! 

7310  IF  R1>10'2S0  THEN 

7320     Rl=RI*10A(-250) 

7330    Scexp=Scexp+250 

7340  END  IF 

7350  IF  Rl<10-(-250)  THEN 

73G0    R1=R1*10*250 

7370     Scexp=Scexp-250 

7380  END  IF 

7390  ! 

7400  ! 

741 0  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7420  I 

7430  !++++    MODULE    CONCLUSION    ++++++++++++++++++++++++++++++++++++++++++++++++ 

7440  I 

7450  SUBEND 

7460  ! 

7470  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7480  i 

7490  i$$$s$$$$$$$$$S$S$S$$$$$$$$$$S$$$$$$$$$$$$$$$$S$$$$$$$$$$S$$$$$$$$$$$$$$ 

7500  ! 
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7510 
7520 
7530 
7540 
7550 
7560 
7570 
7580 
7590 
7600 
7610 
7620 
7630 
7640 
7650 
7660 
7670 
7680 
7690 
7700 
7710 
7720 
7730 
7740 
7750 
7760 
7770 
7780 
7790 
7800 
7810 
7820 


I 

i$$$$ 

SUB  D 

REM 

REM 

REM 

REM 

REM 

REM 

i 

I  ++++ 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

isb(Lnax ,Ka  ,Dy( *  )  ,Y( *  )  ) 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  IRREGULAR  SPHERICAL 
BESSEL  FUNCTIONS  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .   THE 
CALCULATION  OF  THE  DERIVATIVE  USES  THE  ISB'S  OF  THE  SAME  AND 
PREVIOUS  ORDER. 
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

CALCULATING  THE  DERIVATIVES  OF  THE  ISB  +++++++++++++++++++++++++++ 


CALL  Isb(Lnax  ,Ka  ,Y(*  )  ) 

RAD 

Dy(0)  =  (-Y(0)+SIN(Ka  )  >/Ka 

FOR  L=1  TO  Lmax 

Dy(L)  =  -((L+1  )/Ka  )*Y(L)  +  Y<L-1  ) 
NEXT  L 


CALCULATE  THE  ISB  VALUES. 

ANGLE  UNITS  FOR  CALC'NS. 
THE  INITIAL  ISB  DERIVATIVE. 
THE  REMAINING  ISB  DERIVATIVES. 


i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

SUBEND 

i 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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7830  i 

7840  I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$J$S$$$$$$$$$$$$$$$$$$$$$$$$S$$$$$$$S$$$ 

7850  SUB  Isb(Lmax  ,Ka  ,Y( *  )  ) 

7860  REM    +++++++++++++++++++++++++++++++ +++++++-L+++++++++++++++++++++++++++ 

7870  REM    THIS  MODULE  CALCULATES  THE  IRREGULAR  SPHERICAL  BESSEL  FUNCTIONS 

7880  REM    OF  ARGUMENT  KA  AND   ORDER  L  THROUGH  LMAX .    THE  CALCULATION  IS 

7890  REM    ACCOMPLISHED  BY  FORWARD  RECURSION. 

7900  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

7910  i 

7920  I++++  CALCULATING  THE  ISB  ++++++++++++++++++++++++++++++++++++++++++++++ 

7930  ! 

7940  i 

7950  RAD  i  ANGLE  UNITS  FOR  CALC'NS. 

7960  Y(0)=-COS(Ka)/Ka  !  INITIAL  ISB'S  CALCULATED. 

7970  Y(  1  )  =  ( Y(0)-SIN(Ka  )  )/Ka 

7980  FOR  L  =  2  TO  Lmax  !  REMAINING  ISB'S  CALC'D. 

7990    Y(L)  =  ((2*L-1  >/Ka)*Y(L-1  >-Y(L-2) 

8000  NEXT  L 

8010  i 

8020  ! 

8030  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8040  ! 

8050  >++++    MODULE    CONCLUSION    ++++++++++++++++++++++++++++++++++++++++++++++++ 

80G0  ! 

8070  SUBEND 

8080  ! 

8090  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8100  i 

8110  I $$$$$$$$ S$ $$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$SS$$$SS$$$$$$$$$S$$ $$$$$$$$$$$$ 

8120  ! 
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8130  ! 

8140  ! $$$$$$$$$$$$$$$S$$$$$$$$$$$$$$$$$$S$$$$$$$$$$$$$$$S$$$$$$$$$$$$$$$$$$$$ 

8150  SUB  Legl80(Lnax  ,Ka  ,P< *  )  ) 

81G0  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8170  REM    THIS  MODULE  CALCULATES  THE  LEGENDRE  POLYNOMIALS  OF  ARGUMENT  -1  AND 

8180  REM    ORDER  0  THROUGH  LMAX .   CALCULATION  IS  ACCOMPLISHED  BY  FORWARD 

8190  REM    RECURSION. 

8200  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8210  ! 

8220  I++++  CALCULATING  THE  LEGENDRE  POLYNOMIALS  +++++++++++++++++++++++++++++ 

8230  ! 

8240  X=-1  !  ARGUMENT  OF  THE  LEG.  POL'S. 

8250  ! 

82G0  P(0)=1  i  INITIAL  LEG.  POL'S.  CALC'D. 

8270  P( 1  )=X 

8280  FOR  L=2  TO  Lnax  !  REMAINING  LEG.  POL'S.  CALC'D. 

8290    P(L)  =  (2-(  1/L  )  >*X*P(L-1  >-(  1  -<  1  /L  )  )*P(  L-2  > 

8300  NEXT  L 

8310  ! 

8320  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8330  i 

8340  I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

8350  ! 

8360  SUBEND 

8370  ! 

8380  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8390  I 

8400  i$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$s$$$$$$$$$$$S$$$$$$ 

8410  I 
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8420  I 

8430  i $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$S$$ 

8440  SUB  Gain30(Urms ,Tc ,Rh  ,R  .Patrn  ,Ir  .Efreq  ,Za  ) 

8450  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

84G0  REM    THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  EtGo  WITH 

8470  REM    3.0V  INPUT  FROM  THE  POWER  SUPPLY. 

8480  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8490  l 

8500  I++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8510  ! 

8520  ! 

8530  RESTORE  630 

8540  READ  Urns ,Tc ,Rh ,R  ,Pat n  ,  Ir  ,Ef req  ,1a 

8550  G30:    DATA  . 426 ,2  I  . 0 ,50 . 0  ,5 . 30 , 1 02  I . 5 ,2 . 95 ! E-3 ,5000 , 1  5 . 8 

8560  ! 

8570  i 

8580  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8590  i 

8600  I++++  MODULE  CONCLUSION  +  +++  +  +  ++++  ++++++++++  +  ++++++  +  ++++  +++++++  +  +  ->•  ++++++ 

86  10  ! 

8620  SUBEND 

8630  ! 

8640  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8650  i 

8660  i $$$$$$S$$$$$S$S$$$ $$$$$$$ $$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

8670  ! 
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8G80  ! 

86  90  ! $$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

8700  SUB  Ga:n35(Vrms  ,Tc  ,Rh,R  .Patrn  ,Ir  .Efreq  ,Za  ) 

871 0  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8720  REM    THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  EtGo  WITH 

8730  REM    3.5V  INPUT  FROM  THE  POWER  SUPPLY. 

8740  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8750  i 

8760  i ++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8770  i 

8780  ! 

8790  RESTORE  G35 

8800  READ  Urns  ,Tc  ,Rh  ,R  , Patrn  ,Ir  , Efreq  ,Za 

8810  G35:    DATA  . 523 ,2 1 . 0  ,50 . 0  ,5 .30  , 1  02 1 . 5  ,4 . 467E-3  ,5000  ,  1 5 . 8 

8820  I 

8830  i 

8840 

8850 

8860  I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

8870  I 

8880  SUBEND 

8890  i 

8900  ! +++++++++++++++ +++++++++++++++++++++++++++++++++++++++++++++ +++++++++++ 

8910  i 

8930  ! 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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8940  ! 

8950  !$$$$$$$$$$$$$$$$$$$£$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

89G0  SUB  Sphere!  (Dt<»),J,Tc,Patm,Rh,Tgttpt,Dia .Efreq .Nuspkr.Diaspkr  ,Gpre  ) 

8970  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

8980  REM    THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  Er  UITH 

8990  REM    TARGET  NO.  I  EMPLOYED. 

9000  REM    +++++++++++++++++++++++++++++++ +++++++++++++++++++++ ++++++++++++++ 

9010  ! 

9020  •++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

9030  ! 

9040  I 

9050  J=39 

9060  RESTORE  SI 

9070  FOR  1=0  TO  J 

9080  READ  Dt  < I  ) 

9090  NEXT  I 

9100  READ  Tc ,Patm ,Rh  .Tgttm  ,Dia  , Efreq  .Nuspkr .Diaspkr ,Gpre 

9110  SI:   DATA  9. 54  ,  9 . 4  ,8 . 32  ,6 . 42  ,3 . 88  ,  . 96 ,-2 . 07 ,-4 . 9 ,-7. 24 ,-8 . 87 

9120  DATA  -9.64,-9.45,-8.32,-6.41  ,-3 . 83  ,- . 89 ,2 . 1 4 ,4 . 94 ,7 . 26 ,8 . 88 

9130  DATA  9.6 ,9.38.8.22  ,5.26  ,3. 67, .73, -2.3 ,-5. 1 1  ,-7.4  ,-8.96 

9140  DATA  -9.54  ,-9.39,-8.2  ,-6.21  ,-3.61  ,-.65,2.36  ,5. 12 ,7.38,8.91 

9150  DATA  20.2 ,1016.8,52.8,29.91  ,.2546  ,5000,  IS, .0762,1 1094 

9160  i 

9170  ! 

91 80  ! ++4-++++++++++ ++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++ 

9190  ! 

9200  I++++    MODULE    CONCLUSION    ++++++++++++++++++++++++++++++++++++++++++++++++ 

9210  ! 

9220  SUBEND 

9230  I 

9240  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

9250  ' 

9260  i $$$$$$$$$$$S$$$S$SS$$$$$$$$S$$$$$$$$$$$$S$$$$$$$$$$$$S$$$$$$$f $$$$$$$$$ 

9270  ! 
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9280   ! 

92  90   ! $$$$$$$$$$$$$$$$$$$ $$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$ 

9300   SUB  Sphere2(Dt( O.J.Tc.Patm.Rh.Tgttm.Dia  .Efreq  .Nuspkr  .Diaspkr ,Gpre  > 

93 10   REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

9320   REM    THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  Er  WITH 

9330   REM    TARGET  NO.  2  EMPLOYED. 

9340   REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

9350   ! 

9360   I++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

9370  ! 

9380  I 

9390   J=39 

9400   RESTORE  S2 

9410   FOR  1=0  TO  J 

9420     READ  Dt( I  ) 

9430   NEXT  I 

9440   READ  Tc  ,Patm  ,Rh  ,Tgt  trn  ,Dia  ,  Efreq  ,  Nuspkr  ,  Diaspkr  ,Gpre 

9450  S2:   DATA  4 . 72 ,4 . BE  .4 . I 6  .3 . 23  ,2 . 0  I  ,  . 560  ,- . 940  ,-2 . 35  ,-3 . 52 ,-4 . 35 

9460      DATA  -4.77  ,-4 .70  .-4. 1 8  ,-3.25  ,- J . 99  ,- .560  ,. 960  ,2 .34  ,3.51  ,4.34 

9470      DATA  4 . 74  ,4 . 66  ,4 . I  3  ,3 . 1 9  ,  1  . 95  .  . 500  ,- 1 . 00  ,-2 . 40  ,-3 . 57 ,-4 . 38 

9480      DATA  -4 . 78  ,-4 . 68  ,-4 . 1 5  ,-3 . 20  ,- I . 92  ,- . 480  ,  I  . 02  ,2 . 59  ,3 . 54 ,4 . 34 

9490      DATA  2 1 .0  .  1 009. 4  .49. 8  ,30.55  ,  . 1 01 3  ,5000  ,  1 9  ,  . 0762  . 1 1 094 

9500  ! 

9510  ! 

9520 

9530 

9540   i ++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

9550  ! 

9560   SUBEND 

9570  I 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


9580 
9590 
9G00 
96  10 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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9620   ! 

9630   ! $$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$$$$ 

9G40   SUB  Sphere3(Dt(  *),J,Tc,Patm,Rh,Tgttm,Dia  .Efreq  .Nuspkr.Diaspkr.Gpre) 

9650   REM    +++++++++++++++++++++++++^++++++++++++++++++++++++++++++++++++++++ 

9660   REM    THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  Er  UITH 

9670   REM    TARGET  NO.  3  EMPLOYED. 

9680   REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

9690   I 

9700   i ++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

9710  ! 

9720  ! 

9730   J=39 

9740   RESTORE  S3 

9750   FOR  1=0  TO  J 

9760    READ  Dt( I  ) 

9770   NEXT  I 

9780   READ  Tc .Patn ,Rh ,Tgt tm  ,Dia  .Efreq  .Nuspkr  .Diaspkr  ,Gpre 

9790  S3:   DATA  3 . 6  I  ,3 . 40  ,2 . 86  ,2 . 04  ,  1 . 0 1  ,- . I  20  ,- !  . 25 ,-2 . 25  ,-3 . 04  ,-3 . 54 

9800      DATA  -3.68  ,-3.46  .-2.9)  ,-2.07,-1 .04 ,. 1 10,1  .25,2.24,3.00,3.48 

9810      DATA  3. 63, 3. 39, 2. 84, 2. 02,. 970,-. 180,-1 .30,-2.31  ,-3.08,-3.56 

9820      DATA  -3 . 69 ,-3 . 46  ,-2 . 89  ,-2 . 04  ,- . 990 , . 1 60 , 1 . 29 ,2 . 2G  ,3 . 05  ,5. 50 

9830      DATA  2  1  . 3 , 1 009 . 6 ,49 . 8 ,30 . 59  , . 07634 ,5000 , 1 9 . . 0762  ,  1 1 094 

9840  ! 

9850  ! 

9860       !  +++++++++++++  +  +  +++  +  +++  ++++++++ -"-++++++++++++++++  +  +  +++++  ++++++++  +++++++++ 

9870       ! 

9880   I++++  MODULE  CONCLUSION  +++++  ++++++++++++++++  +  +++  +  +++-»-  +++  ++++++++++++++  + 

9890  ! 

9900   SUBEND 

9910  ! 

9920 

9930 

9940 

9950 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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99G0   ' 

9970   I $$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ J $ 

9980   SUB  Sphere4(Dt( *).J,Tc,Patn,Rh,Tgttm,Dia  .Efreq  .Nuspkr.Diaspkr.Gpre) 

9990   REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10000  REM    THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  Er  WITH 

10010  REM    TARGET  NO.  4  EMPLOYED. 

10020  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10030  ! 

10040  !++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10050  l 

100G0  i 

10070  J=39 

10080  RESTORE  S4 

10090  FOR  1=0  TO  J 

10100    READ  Dtd  ) 

101 10  NEXT  I 

10120  READ  Tc ,Patm ,Rh .Tgttm ,Dia  , Efreq  .Nuspkr  .Diaspkr  ,Gpre 

10130  S4:  DATA  2 . 28  ,2 . 1 0  ,  1  . G8  ,  I  . 1 0  ,  . 420  ,- . 320  ,- 1  . 03  ,- 1 . G4  ,-2 . 07  ,-2 . 3 1 

10140     DATA  -2.34 ,-2. 12  ,-1 .71  ,-1  . 12  ,-.440,.  320  ,1  .02  ,1 .62,2.06,2.28 

10150     DATA  2.29,2. 10,1  .68,1  .09,. 400,-. 340,-1. .04,-1 .66 ,-2. 10,-2.32 

10160     DATA  -2.35,-2. 12  ,-!  .70,-1 . 1 1  ,-.400,.  340,1  .04  J .54,2.08,2.28 

10170     DATA  21  .4,1009.9,49.6  ,30. 54, .06240, 5000, 19, .0762  ,1 1094 

10180  ' 

10190  ! 

1 0200  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10210  ! 

10220  I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

10230  ! 

10240  SUBEND 

10250  ! 

10260  ! ++++++++++++++++++++++++++++++++++++++++++++++ +++++++++++++++++++++++++ 

10270  i 

1 0280  I $$$$$$$$$$$$$$$$$$$$$$$$$$$$£$$$$$$$$$$$$$$$$$£$$$$$$$$$$$$$$$$$$ $$$$$$ 

10290  i 
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10300  ! 

10310  !$$$$$$$$$$$$$$$$$$$$ $$£$$$$$$$$$$$$$$$$$$$$$$$$$S$$$$$$$$$$$$$$$$$$$$$$ 

10320  SUB  Trans30( J  ,Dt< «  )  ) 

1 0330  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10340  REM    THIS  MODULE  CONTAINS  DATA  USED  IN  DETERMINATION  OF  Er  WITH 

10350  REM    3.0V  APPLIED. 

1 0360  REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10370  i 

10380  ! ++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10390  ! 

10400  ! 

10410  J=39 

10420  RESTORE  T30 

10430  FOR  1=0  TO  J 

10440     READ  Dt( I  ) 

10450  NEXT  I 

10460  T30:  DATA  . 630  ,  . 680  .  . 590  , . 440  ,  . 320  , . 1 60  ,- . 020  ,- . 230  ,- . 4 1 0 ,- . 550 

10470      DATA  -.G 10, -.6 10  .-.520, -.360 ,-.240 ,-.080,. 1 1 0 , . 320 , . 490 , . 630 

10480      DATA  .6  90  ,  .680,  .590  ,.430,. 320,. 150 ,-.040 ,- . 240 ,- . 420 ,- . 560 

10490      DATA  -.620  ,-.600  ,-.510  ,-.380  ,- . 240  ,- . 070 , . 120,. 330, .500,. 6 40 

10500  ! 

10510  ! 

1 0520  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10530  i 

10540  I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

10550  ! 

10560  SUBEND 

10570  ! 

10580  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10590  I 

10600  i $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$ 

10610  I 
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10G20    ! 

10630    !$$$£$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

10G40    SUB    Trans35( J  ,Dt( *  )  ) 

1 0G50   REM        ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10660   REM        THIS    MODULE    CONTAINS    DATA   USED    IN   DETERMINATION   OF    Er    WITH 

10G70   REM        3.SV   APPLIED. 

1 0680   REM        ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10690    ! 

10700    !++++   LOADING    DATA    +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10710  ! 

10720  ! 

10730    J=39 

10740   RESTORE    T35 

107S0   FOR    1=0    TO    J 

10760  READ   Dt( I > 

10770    NEXT    I 

10780   T35:      DATA    . 830 , . 770 . . 640  , . 490  , . 290 , . 080 ,- . 1 80 ,- . 400 ,- . 6 1 0 ,-. 730 

10790  DATA    -.770, -.720, -.580, -.430  ,-.240, -.020,. 240,. 460,. 660,. 780 

10800  DATA    .830,. 770  ,.630,. 470,. 290,. 060,-. 1 90 ,- . 420 ,- . 600 ,- . 720 

10810  DATA   -.770, -.7 10, -.570, -.420, -.230, -.010,. 240,. 480,. 670.. 790 

10820  ! 

10830  ! 

1 0840    I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

10850 

10860 

10870 

10880    SUBEND 

10890 


++++    MODULE    CONCLUSION    ++++++++++++++++++++++++++++++++++++++++++++++++ 


10900 
10910 
10920 
10930    i 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

ft  ft  ft  ft  m5  ft  41 4i  ft  ft  45  45  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  45  ft  ft  45  45  45  ft  ft  45  ft  ft  ft  45  41 45  4!  ft  45  45 ft  45  ft  ft  4!  ft  ft  ft  ft  ft  ft  ft  3i  ^  't  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft 
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10940 
10950 
10960 
10970 
10980 
10990 
I  1000 
I  1010 
1  1020 
1  1030 
1  1040 
1  1050 
1  1060 
1  1070 
1  1080 
1  1090 
1  1  100 
11110 
1  1  120 
1  1  130 
1  1140 
1  1150 
1  1  160 
1  1  170 
1  1  180 
1  11  90 
1  1200 
11210 
1  1220 
1  1230 
1  1240 
1  1250 
1  1260 
1  1270 
1  1280 
1  1290 
1  1300 
1  1310 
1  1320 
1  1330 
1  1340 
1  1350 
1  1360 
1  1370 
1  1380 


!$S$$S$$$$$$$$$$$$$$$$$$$$$$$$$$$S$$$$$$$$$$$$$$S$$$$S$$$$$$$$$$$$$$$$$$ 
SUB  Noise30( J , Jt ,Dt< * ) ,Tdt< *  >  .Ten ,Rhn .Patmn .Efreqn) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM    THIS  MODULE  CONTAINS  DATA  USED  IN  DETERMINATION  OF  Er  WITH 
REM   3.0U  APPLIED. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
l 

i ++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

'  ) 
I 
J  =  59 

RESTORE  N30 
FOR  1=0  TO  J 

READ  Dt( I  ) 
NEXT  I 
Jt=39 

RESTORE  Tn30 
FOR  1=0  TO  Jt 

READ  Tdt( I  ) 
NEXT  I 

READ  Ten  ,Rhn  .Patmn  .Efreqn 

N30:  DATA  0 ,0 ,- . 02  ,- .03  ,- . 04  ,- . 06  ,- . 08 ,- . 1  ,-  .  1 1  ,- . 1 2 
DATA  -.12,-. 11  .-.1  ,-.08, -.07 ,-.05, -.04, -.01  ,0.0 
DATA  0, 0,-. 02,-. 03, -.05, -.07,-. 08, -.1  ,-. 1 1  ,-. 12  • 
DATA  -. 12  ,-. 12  ,-. 1  1  ,-.08, -.07, -.04 ,-.03, -.02 ,0,0 
DATA  .01  ,0,0,-. 03, -.04, -.07, -.09, -.11  ,-.12, -.12 
DATA  -.12, -.12, -.11, -.08, -.07, -.04, -.03, 0,0,. 01 
Tn30:  DATA  . 69  ,  . 68  ,  . 59  ,  . 44  ,  .  32  , . 1 6  ,- . 02 ,-. 23 ,- . 4 1  ,- . 55 

DATA  -.61  ,-.61  ,-.52  ,-'.36  ,-.24, -.08,.  1  1  ,.32  ,.49,.  63 
DATA  .69,. 68,. 59,. 43,. 32,. 15, -.04, -.24. -.42, -.56 
DATA  -.62  ,-.6  .-.51  ,-.38, -.24, -.07,. 12 ,.33,. 5,. 64 
DATA  20.8,51 .7,1016.7,5000 

I 
i 
i ++++++++++++++++++++ +++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

i ++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

! 
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

Is  5  xi^55xxx>Sxx$SSS$S*$$$S£xx$S  —  $S$SSSSS5$£$SSxS5S5SxSx5Sx5£x**$5S$SS$5S 
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!  1390 
I  1400 
1  1410 
1  1420 
1  1430 
1  1440 
I  1450 
1  1460 
1  1470 
1  1480 
1  1490 
1  1500 
1  1510 
1  1520 
1  1530 
1  1540 
1  1550 
1  1560 
1  1570 
1  1580 
1  1590 
1  1600 
1  1610 
1  1620 
1  1630 
1  1640 
1  1650 
1  1660 
1  1670 
1  1680 
1  1690 
1  1700 
1  1710 
1  1720 
1  1730 
1  1740 
1  1750 
1  1760 
1  1770 
1  1780 
1  1790 
1  1800 
1  1810 


I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$ $$$$$$$ 
SUB  Noise35( J , Jt ,Dt ( ♦  )  ,Tdt ( ♦  )  ,Tcn ,Rhn .Patmn ,Ef reqn ) 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM    THIS  MODULE  CONTAINS  DATA  USED  IN  DETERMINATION  OF  Er  WITH 
REM    3.5V  APPLIED. 

REM    ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++  LOADING  DATA  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 
I 
J=39 

RESTORE  N35 
FOR  1=0  TO  J 

READ  Dt( I  ) 
NEXT  I 
Jt=39 

RESTORE  Tn3S 
FOR  1=0  TO  Jt 

READ  Tdt( I  ) 
NEXT  I 

READ  Ten  ,Rhn  ,Patnn  .Efreqn 

N35:   DATA  - . 06 ,- . 08 ,- . 08  ,- . 08  ,- . 07  ,- . 06  ,- . 04  ,- . 02  .0  ,  . 02 
DATA  .05,. 06 ,.06 ,.06 , .04 , .03 , .02 ,0 ,-.03 ,-.05 
DATA  -.07, -.08, -.09, -.08, -.08, -.06  ,-.04, -.02  ,.0!  ,.03 
DATA  .04 , .07,. 06  ,  .07  ,  .05,  .04,. 01  ,-.01  ,-.04  ,-.06 
Tn35:  DATA  .74 , . 83 , .82 , . 72 , . 53  , . 38  ,  . 1 8  ,-.06  ,-. 31  ,- .53 
DATA  -.68, -.77, -.75, -.64, -.46  ,-.3,-. 1  ,. 14,. 4,. 6 
DATA  .75,. 83,. 81  ,.7,. 5,. 36,. 16  ,-.08, -.33,-. 54 
DATA  -.6 9, -.76 ,- . 74 ,-. 63  ,- . 46  ,-.28, -.08,. 16  ,.4  1  ,.62 
DATA  21  .3,51 .3,1016.6  ,5000 

I 
I 
i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
i 

'++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

! 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

$$$$$$$$$$$$$$$$$$$ $$$$$s$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ss$$$$$$$$$$$ 
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APPFND1X    F 
ACOUSTIC   ARRAY  CALIBRATION    DATA 

This  appendix  contains  data  used  to  determine  the  calibration 
product  Er  Et.    Data  units  are  contained  in  square  brackets. 

For  listings  of  the  same  type  of  data  the  units  and  uncertainty  are 
explicitly  given  for  the  first  datum.  Remaining  data  in  the  listing  have  the 
same  units  and  uncertainty  as  the  first  datum. 

1 .   TWO  WAY  PROPAGATION  PATH  DATA 

Pre-amplifier  gain  Ge  and  array  impedance  Z  were  previously 

determined  by  Moxcey.  Gewas  determined  to  be  1 1,094  Z  was  determined 

to  be  15.8[0].  [Ref.  8] 

The  remaining  data  are  presented  in  the  following  tables. 

Vsu   ,    is  the  voltage  supplied  by  the  HP3314A  function  generator  to 

the  pre-amplifier  (see  Figure  8). 
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Table  5a:  Two  Way  Propagation  Path  Data  for  Sphere  One, 
Collected  30  November  1987;  Run  *2:  Data  Re-measured 
Because  First  Run  V  (with  3.5  [V]  Supply)  Was  Clipped. 


Diameter  =  2a  =  .2546  +  0.0006  [m] 
Tc  =  20.2  ±  0. 1  [°  c] 

Rh  =  52.8  ±  0. 1  [%] 
P=  1016.8  +  0.2  [mb] 
No.  pulses  =  20 


v 

supply 

-  3.0  [V] 

tr  =  30.31  [mS]- 

400  [u.Sl  =  29.91 

+  0.1  [mS] 

t 

Vt 

t 

Vt 

2.25  + 

O.OHmS] 

690  +  5  [mV] 

2.45 

690 

2.26 

680 

2.46 

680 

2.27 

590 

2.47 

590 

2.28 

440 

2.48 

430 

2.29 

320 

2.49 

320 

2.30 

160 

2.50 

150 

2.31 

-20 

2.51 

-40 

2.32 

-230 

2.52 

-240 

2.33 

-410 

2.53 

-420 

2.34 

-550 

2.54 

-560 

2.35 

-610 

2.55 

-620 

2.36 

-610 

2.56 

-600 

2.37 

-520 

2.57 

-510 

2.38 

-360 

2.58 

-380 

2.39 

-240 

2.59 

-240 

2.40 

-80 

2.60 

-70 

2.41 

1  10 

2.61 

120 

2.42 

320 

2.62 

330 

2.43 

490 

2.63 

500 

2.44 

630 

2.64 

640 

141 


Table  5b:  Two  Way  Propagation  Path  Data  for  Sphere  One  Continued, 
Collected  30  November  1987,  Run  *2:  Data  Re-measured 
Because  First  Run  Vr  (with  3.5  [V]  Supply)  Was  Clipped. 


t  V 


r 


33.68  +0.0 1  [mS]     9.54  ±  0.03  [V] 

33.69  9.40 

33.70  8.32 

33.71  6.42 

33.72  3.88 

33.73  0.960 

33.74  -2.07 

33.75  -4.90 

33.76  -7.24 

33.77  -8.87 

33.78  -9.64 

33.79  -9.45 

33.80  -8.32 

33.81  -6.41 

33.82  -3.83 

33.83  -0.890 

33.84  2.14 

33.85  4.94 

33.86  7.26 

33.87  8.88 


t 

V 

r 

33.88 

9.60 

33.89 

9.38 

33.90 

8.22 

33.91 

6.26 

33.92 

3.67 

33.93 

0.730 

33.94 

-2.30 

33.95 

-5.11 

33.96 

-7.40 

33.97 

-8.96 

33.98 

-9.64 

33.99 

-9.39 

34.00 

-8.20 

34.01 

-6.21 

34.02 

-3.61 

34.03 

-0.650 

34.04 

2.36 

34.05 

5.12 

34.06 

7.38 

34.07 

8.91 
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Table  6a:  Two  Way  Propagation  Path  Data  for  Sphere  Two, 
Collected  28  November  1987. 


Diameter  =  2a  =  1 0. 1 3  ±  0.0 1  [cm] 
Tc  =  21.0±  0.1  [°c] 

Rh  =  49.8  ±  0. 1  [%] 
P=  1009.4+  0.2  [mb] 
No.  pulses  =  20 


v 

supply 

=  3.5  [V] 

tr  =  30.95  [mS]  - 

400 [nS] = 

30.55  t 

0.1  [mS] 

t 

Vt 

t 

Vt 

2.26  ± 

0.0  l[r 

TlS] 

830  i  5  [i 

TlV] 

2.46 

830 

2.27 

770 

2.47 

770 

2.28 

640 

2.48 

630 

2.29 

490 

2.49 

470 

2.30 

290 

2.50 

290 

2.31 

80 

2.51 

60 

2.32 

-180 

2.52 

-190 

2.33 

-400 

2.53 

-420 

2.34 

-610 

2.54 

-600 

2.35 

-730 

2.55 

-720 

2.36 

-770 

2.56 

-770 

2.37 

-720 

2.57 

-710 

2.38 

-580 

2.58 

-570 

2.39 

-430 

2.59 

-420 

2.40 

-240 

2.60 

-230 

2.41 

-  20 

2.61 

-    10 

2.42 

240 

2.62 

240 

2.43 

460 

2.63 

480 

2.44 

660 

2.64 

670 

2.45 

780 

2.65 

790 
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Table  6b:  Two  Way  Propagation  Path  Data  for  Sphere  Two  Continued, 
Collected  28  November  1987. 


t 

V 

r 

t 

V 
r 

3413  ±0.01  [mS] 

4.72  +  0.01  [V] 

34.33 

4.74 

3414 

4.66 

34.34 

4.66 

3415 

4.16 

34.35 

4.13 

3416 

3.23 

34.36 

3.19 

3417 

2.01 

3437 

1.95 

3418 

0.560 

3438 

0.500 

3419 

-0.940 

34.39 

-1.00 

3420 

-2.35 

34.40 

-2.40 

3421 

-3.52 

34.41 

-3.57 

3422 

-4.36 

3442 

-4.38 

3423 

-4.77 

34.43 

-4.78 

3424 

-4.70 

3444 

-4.68 

3425 

-4.18 

3445 

-4.15 

3426 

-3.25 

34.46 

-3.20 

3427 

-1.99 

34.47 

-1.92 

3428 

-0.560 

3448 

-0.480 

3429 

0.960 

34.49 

1.02 

34.30 

2.34 

34.50 

2.39 

34.31 

3.51 

3451 

3.54 

3432 

4.34 

34.52 

4.34 
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Table  7:  Two  Way  Propagation  Path  Data  for  Sphere  Three, 
Collected  28  November  1987. 


Diameter  =  2a  =  7.634  +  0.003  [cm] 
Tc  =  21.3  ±0.1  [°c] 

Rh  =  49.8  i  0. 1  [%] 
P=  1009.6  +  0.2  [mb] 
No.  pulses  =  20 
V      ,    =3.5[V] 

supply 

tr  =  30.99  [mS]  -  400  [uS]  =  30.59  ±  0. 1  [m5] 


Vt  data  is  the  same  as  sphere  two's  (see  Table  6a). 


t 

Vr 

t 

Vr 

33.98  ±0.0  1[mS] 

3.61  ±0.03[V] 

34.18 

3.63 

33.99 

3.40 

34.19 

3.39 

34.00 

2.86 

34.20 

2.84 

34.01 

2.04 

34.21 

2.02 

34.02 

1.01 

34.22 

0.970 

34.03 

-0.120 

34.23 

-0.180 

34.04 

-1.25 

34.24 

-1.30 

34.05 

-2.25 

34.25 

-2.31 

34.06 

-3.04 

34.26 

-3.08 

34.07 

-3.54 

34.27 

-3.56 

34.08 

-3.68 

34.28 

-3.69 

34.09 

-3.46 

34.29 

-3.46 

34.10 

-2.91 

34.30 

-2.89 

34.11 

-2.07 

34.31 

-2.04 

34.12 

-1.04 

34.32 

-0.990 

34.13 

0.110 

34.33 

0.160 

34.14 

1.25 

34.34 

1.29 

34.15 

2.24 

34.35 

2.26 

34.16 

3.00 

34.36 

3.05 

34.17 

3.48 

34.37 

3.50 
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Table  8:  Two  Way  Propagation  Path  Data  for  Sphere  Four, 
Collected  28  November  1987. 


Diameter  =  2a  =  6.240  ±  0.008  [cm] 
Tc  =  21.4±  0.1  [°c] 

Rh  =  49.6  +  0. 1  [%] 
P=  1 009.9  ±  0.2  [mb] 
No.  pulses  =  20 
V      ,    =3.5[V] 

supply 

tr  =  3 1 .04  [mS]  -  400  [jiS]  =  30.64  ±  0. 1  [mS] 

Vt  data  is  the  same  as  sphere  twos  (see  Table  6a). 


t 

Vr 

t 

Vr 

33.83  ±0.0 1[mS] 

2.28  i  0.03  [V] 

34.03 

2.29 

33.84 

2.10 

34.04 

2.10 

33.85 

1.68 

34.05 

1.68 

33.86 

1.10 

34.06 

1.09 

33.87 

0.420 

34.07 

0.400 

33.88 

-0.320 

34.08 

-0.340 

33.89 

-1.03 

34.09 

-1.04 

33.90 

-1.64 

34.10 

-1.66 

33.91 

-2.07 

34.11 

-2.10 

33.92 

-2.31 

34.12 

-2.32 

33.93 

-2.34 

34.13 

-2.35 

33.94 

-2.12 

34.14 

-2.12 

33.95 

-1.71 

34.15 

-1.70 

33.96 

-1.12 

34.16 

-1.11 

33.97 

-0.440 

34.17 

-0.400 

33.98 

0.320 

34.18 

0.340 

33.99 

1.02 

34.19 

1.04 

34.00 

1.62 

34.20 

1.64 

34.01 

2.06 

34.21 

2.08 

34.02 

2.28 

34.22 

2.28 
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Table  9:  Noise  Data  for  3.0  [V]  Supply  Voltage, 
Collected  30  November  1987. 


Tc  =  20.8  ±0.1  [° 

c]      Rh 

=  51.7+  0.1  [%] 

P=  1016.7 

±  0.2  [rr 

V 
supply 

-  3.0  [V] 

No. 

pulses  =  20 

Vt  data 

is  the  same  as  sphere  one's  (see 

Table  5a). 

t 

V 

n 

t 

Vn 

33.73  ± 

O.OllmS] 

00  i 

5[mV] 

34.03 

-120 

33.74 

00 

34.04 

-120 

33.75 

-  20 

34.05 

-110 

33.76 

-  30 

34.06 

-  80 

33.77 

-   40 

34.07 

-  70 

33.78 

-  60 

34.08 

-  40 

33.79 

-  80 

34.09 

-  30 

33.80 

-100 

34.10 

-  20 

33.81 

-1  10 

34.11 

00 

33.82 

-120 

34.12 

00 

33.83 

-120 

34.13 

10 

33.84 

-110 

34.14 

00 

33.85 

-100 

34.15 

00 

33.86 

-  80 

34.16 

-  30 

33.87 

-  70 

34.17 

-  40 

33.88 

-  50 

* 

34.18 

-  70 

33.89 

-   40 

34.19 

-  90 

33.90 

-    10 

34.20 

-1  10 

33.91 

00 

34.21 

-120 

33.92 

00 

34.22 

-120 

33.93 

00 

34.23 

-120 

33.94 

00 

34.24 

-120 

33.95 

-  20 

34.25 

-1  10 

33.96 

-   30 

34.26 

-  80 

33.97 

-   50 

34.27 

-  70 

33.98 

-   70 

34.28 

-  40 

33.99 

-   80 

34.29 

-   30 

34.00 

-100 

34.30 

00 

34.01 

-1  10 

34.31 

00 

34.02 

-120 

34.32 

10 
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Table  10a.  Noise  Data  for  3.5  [V]  Supply  Voltage; 
Collected  30  November  1987. 


Tc  =  21.3  ± 

0.1  [° 

c] 

Rh  =  51.3  i 

0.1  [%] 

P=  1016.6 

i  0.2  [mb] 

No.  pulses  ; 

=  20 

V      ,    =  3.5  [V] 

supply 

t 

vt 

t 

Vt 

2.24  ±0.0 1[mS] 

740  +  5  [mV] 

2.44 

750 

2.25 

830 

2.45 

830 

2.26 

820 

2.46 

810 

2.27 

720 

2.47 

700 

2.28 

530 

2.48 

500 

2.29 

380 

2.49 

360 

2.30 

180 

2.50 

160 

2.31 

-   60 

2.51 

-   80 

2.32 

-310 

2.52 

-330 

2.33 

-530 

2.53 

-540 

2.34 

-680 

2.54 

-690 

2.35 

-770 

2.55 

-760 

2.36 

-750 

2.56 

-740 

2.37 

-640 

2.57 

-630 

2.38 

-460 

2.58 

-460 

2.39 

-300 

2.59 

-280 

2.40 

-100 

2.60 

-  80 

2.41 

140 

2.61 

160 

2.42 

400 

2.62 

410 

2.43 

600 

2.63 

620 
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Table  10b:  Noise  Data  for  3.5  [V]  Supply  Voltage  Continued, 
Collected  30  November  1987. 


t 

Vn 

t 

Vn 

33.80  ±0.0 1  [mS] 

-  60  ±  5  [mV] 

34.00 

-  70 

33.81 

-  80 

34.01 

-  80 

33.82 

-  80 

34.02 

-  90 

33.83 

-  80 

34.03 

-  80 

33.84 

-  70 

34.04 

-  80 

33.85 

-  60 

34.05 

-  60 

33.86 

-  40 

34.06 

-  40 

33.87 

-  20 

34.07 

-  20 

33.88 

00 

34.08 

10 

33.89 

20 

34.09 

30 

33.90 

50 

34.10 

40 

33.91 

60 

34.11 

70 

33.92 

60 

34.12 

60 

33.93 

60 

34.13 

70 

33.94 

40 

34.14 

50 

33.95 

30 

34.15 

40 

33.96 

20 

34.16 

10 

33.97 

00 

34.17 

-  10 

33.98 

-  30 

34.18 

-  40 

33.99 

-  50 

34.19 

-  60 
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2.      ONE  WAY  PROPAGATION  PATH  DATA 

One  way  propagation  path  data  are  presented  in  the  following  tables. 
T.  and  Rh  are  estimated.  The  estimates  were  based  on  hygrothermographic 

comparisons  and  28  and  30  November  data. 


Table  11:  One  Way  Propagation  Path  Data  for  3.0  [V], 
5000  [Hz]  Continuous  Wave  Supply  from  the  HP3314A 
Function  Generator;  Collected  7  December  1987. 


(V      ).  =0.246  +  0.001  [V] 
rms  i 

Tc~21.0[°c] 

Rh  ~  50.0  [%] 
P=  1021.5+  0.2  [mb] 

R  =  5.30  [m]  (calculated  from  two  way  propagation  path  data  for 
spheres  two,  three  and  four) 

I   (referenced  to  10~12[W/m2]) 

r 

947  ±  0.1  [dB] 

945 

947 

95.0 

948 

95.1 

942 
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Table  12:  One  Way  Propagation  Path  Data  for  3.5  [V], 
5000  [Hz]  Continuous  Wave  Supply  from  the  HP33 1 4A 
Function  Generator;  Collected  7  December  1987. 


(Vpms)t  =  0.523i  0.001  [V] 

Tc~21.0[°c] 

Rh  *  50.0  [%] 
P=  1021.5+  0.2  [mb] 

R  =  5.30  [m]  (calculated  from  two  way  propagation  path  data  for 
spheres  two,  three  and  four) 

I   (referenced  to  10"12  [W/m2]) 

r 

97.4  1  0. 1  [dB] 

96.4 

95.1 

97.0 

96.2 

97.1 

96.4 

95.4 

96.3 

97.4 
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